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FOREWORD 


One of the main goals of the Float Zone (FZ) growth project of NASA's 
Materials Processing In Space Program Is to thoroughly understand the molten 
zone/ freezing crystal system and all the mechanisms that govern this system. To 
accomplish this, the melt and Interface properties, the heat and mass flows, and 
the dependencies of these on each other and on growth rate and g levels must be 
studied. 


Since the float zone process Involves two solid-melt Interfaces, 
possible gas Interfaces, heat and mass transfers, various driving forces and 
complex heating sources, an analysis of the entire process would be very 
complex. For an Initial investigation, a more feasible approach is to examine 
each component of the process separately, particularly If mathematical models 
are to be manageable. The three principal components are: (1) the shapes of 
the melt and sclld-melt Interfaces, (2) the heat and mass transfers, and (3) the 
heating and cooling sources. This study combined facets of all three 
components . 

The purpose of this 12-month effort was to study and compute the 
surface boundary conditions required to give flat FZ solid-melt Interfaces. The 
successful completion of this study should provide FZ furnace designers with 
better methods for controlling solid-melt interface shapes and for computing 
thermal profiles and gradients. 

This study was undertaken in two phases. The first phase was to 
investigate the solid zones surface boundary conditions required for flat 
solid-melt Interfaces when given the melt zone surface boundary conditons. The 
second phase complemented the first and was to investigate the melt zone surface 
boundary conditions required for flat solid-melt interfaces if given the solid 
zones surface boundary conditions. Dual integral transform methods were used in 
both phases; In addition, the use of various numerical methods for differential 
equations and linear systems of equations were required. 

Using NASA supplied data, the surface boundary conditions required for 
flat solid-melt Interfaces were studied. In addition, complete documentation 
and a simple user's guide are provided for all the computer software required 
during this study. 
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1.0 INTRODUCTION 


1 . 1 OVERVIEW AND STUDY DEFINITION 

Silicon (Si) Is used in a wide variety of electronic devices including 
high power rectifiers, space solar cells. Infrared detector arrays and high 
density integrated circuits. The three principal industrial methods for growing 
silicon crystal Ingots or boules are the float zone (FZ), Czochralski (Cz), and 
cold crucible methods. Because molten silicon acts as a universal solvent, Cz 
grown Si is plagued with crucible contamination which is Intolerable for high 
performance optical and infrared devices. However, because the FZ process is 
containerless, crucible contaminants are avoided. Other advantages of FZ growth 
include uniformity of axial resistivity (on a macroscale), visibility of the 
growth region, low consumable material costs, and high growth rates. Although 
the cold crucible method combines many of the best features of the FZ and Cz 
techniques, the molten Si must be superheated and volatile dopants such as In, 
Ga, and T1 are unfortunately evaporated. 

Because most industrial advances in the FZ growth technologies have 
come about empirically, detailed analysis of the growth process has not kept 
pace with presently used FZ methods. Theoretical modeling of the melt dynamics 
has led to some understanding of the growth process, but it is very Incomplete. 
The characteristics of the FZ melt must be more accurately modeled if an 
understanding of the heat balance and flow, Isotherm shapes, density (including 
inversion) and surface tension variations is to lead to better methods of 
controlling the growth conditions. Moreover, such studies should contribute to 
the design and execution of FZ experiments in low^gravity (g) environments, in 
addition, knowledge gained by studying silicon FZ methods should be applicable 
to other FZ processed materials. 

As noted by E. Kem (1Q| , the main goal of the FZ growth project of 
NASA's Material Processing in Space program is to thoroughly undr stand the 
molten zone/freezing crystal system and all the mechanisms that govern this 
system. In addition, more optimal crystal growth conditions at g-1 and possible 
improvements made by processing in near zero-g environments need to be 
investigated. To accomplish this, the melt and interface properties, the heat 
and mass flows, and the dependencies of these on each other and on the growth 
rate and g levels must be studied. 

To transform a polycrystalline material into a single crystal, it is 
not always necessary to melt the entire sample or charge before growing the 
desired monocrystal. In some cases, it is possible to melt a small portion of 
the original charge, translate this molten zone through the charge, and 
hopefully leave a monocrystal behind the translating molten zone. The actual 
heating sources for this type zone melting process are varied and include 
induction, resistance, electron beam, and laser beam. The molten zone itself 
can be moved through the charge by either moving the heating source over the 
charge or by moving the charge through the heating source. The actual charge 
may, but need not be, contained in some type of crucible or ampoule. If no 
container is Involved, the technique is called a float zone method and is used 
for reactive or high melting point materials. For most float zone applications. 
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the molten zone Is held Intact by surface tension with the occasional aid of a 
magnetic field [6]. A simple illustration of the float zone technique Is given 
In Figure 1-1. 


In order to reduce nonuniformities In such things as resistivities and 
defect distributions, for example, the entire solid-melt system must be 
characterized. Since the float zone process Involves two solid-melt Interfaces, 
possible gas Interfaces, heat and mass transfers, various driving forces and 
complex heating sources, an analysis of the entire float zone process would be 
very complex. A more feasible approach (at least for an Initial Investigation) 
Is to examine each component of the system separately, particularly If the 
mathematical models are to be manageable. Three principal system components 
are: (1) the shapes of the melt and solid-melt Interfaces, (2) the heat and 

mass transfers, and (3) the heating and cooling sources. This study combines 
facets of all three components. 

vniile many investigators, e.g., R. Brown [2], R. Naumann [l^ , and W. 
Wilcox [20], are making significant progress studying t*^e solid-melt interface 
shapes and the thermal gradients at the solid-melt Interfaces for float zone and 
analogous systems subject to specified surface boundary conditions, the 
principal thrust of this effort was to study and compute the surface boundary 
conditions required to give flat FZ solid-melt Interfaces. 

The completion of this study hopefully results In a better 
of the FZ diffusion and growth mechanisms and should provide FZ 
furnace designers with better methods for controlling solid-melt interface 
shapes and for computing thermal profiles and gradients. Moreover, the methods 
developed in this study should aid In the design of FZ heaters that achieve the 
required melt fluxes with minimal energy expenditures and, hence, perhaps reduce 
the system power requirements (a natural concern for any long-term, low-g FZ 
experiment). In particular, if radio frequency heating is used, the methodology 
developed ' In this study should be useful for computing the performance 
requirements and position of auxiliary heating and Insulation required for the 
proper thermal profiles. In addition, the methodology developed in this effort 
might provide, for future studies, a starting point for the more complex and 
realistic case of a slightly concave solid-melt Interface. 

This study was performed In two phases. The first phase analyzed the 
solid zones' surface boundary conditions required for flat solid-melt Interfaces 
when given (a priori) the melt zone surface boundary conditions. The second 
phase complemented Che first and analyzed the melt zone surface boundary 
conditions required for flat solid-melt Interfaces when given (a priori) the 
surface boundary conditions for the solid zones. Dual Integral transform methods 
were used In both phases; in addition, both phases required the use of various 
numerical methods for boundary value problems. Ithough such a study has 
ap'parently never before been undertaken, analogous studies for the 
Brldgman-Stockbarger method have been completed by L. Foster [7], [gj. 

Mathemay.cal descriptions of the problems posed above are stated In 
Section 1-2. In CTiapter 2, various mathematical tools are developed followed by 
some rather Interesting examples. The methodologies used to compute the surface 
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boundary conditions fo'* Che solid regions and melt zone surface boundary 
condlclons required for flac Interface shapes are developed In Chapters 3 and 4. 
The results of various test cases using NASA supplied data are presented In 
Chapter 5 and recoomendaclons for future efforts are given in Chapter 6. A 
simple user's guide to various computer codes (listed In Appendix C) 
Implementing the methods described In Chapters 2, 3, and 4 Is presented In of 
Appendix A. Appendix B contains the proof of a claim made In Section 2.3. 

1.2 MATHEMATICAL STATEMENT OF THE CONTROL PROBLEMS 

Concise mathematical statements of the problems described in the 
previous section are given next using Che numbered equations In the FZ model 
shown In Figure 1-2. First suppose that Che melt zone surface temperature is 
some a priori known (by design or happenstance) distribution h(x) (Figure I-2. 
Equation (FZ8)). Since the temperature at both of the assumed flat solid-melt 
interfaces is the material melting point (Equations (FZl) and (FZ3)), the 
temperature distribution In Che melt zone is known and may be computed by Che 
method described In Section 2.2. Hence, Che axial thermal gradients in the melt 
zone at both of the solid-melt Interfaces are known (see Section 2.2). Invoking 
Equations (FZ2) and (FZ4), Che solid regions' axial thermal gradients at the 
Interfaces are also known.*^ 

For the moment, consider the lower solid region (x<0 In Figure 1-2) 
and let B(r) denote the known required thermal gradient* in the solid region at 
the interface (x*0), i.e.. 


T^(O.r) - B(r) , 0 < r < 1 (1.2.1) 

The basic idea is to compute a temperature distribution f(x), x<0, (henceforth 
called a surface control function) , to be maintained on the surface of Che lower 
solid region such that Che resulting temperature distribution, T(x,r), for the 
lower solid region satisfies Equation (1.2.1). This is concisely stated in 
Problem Pl-1. 


t Equations (FZ2) and (FZ4) of Figure 1-2 guarantee Che conservation of energy 
at the solid-melt interfaces (x-0 and x-Q) . kj^ and kg are the solid and liquid 
thermal conductivities while A is product of the growth rate, solid density and 
latent heat of fusion [15l • 


* Standard machemaclcal nomanelature La used in this report. Both the operator 
and subscript notation are used for partial derivatives, e.g., ^ and T^ both 
denote the partial derivative of T(x,r) with respect to x. For functions of 
one variable, Che "prime" convention for derivatives Is observed, e.g., h''(x) 
denotes Che second derivative of h(x). The Laplacian operator is denoted by 
A and is, in cylindrical coordinates, 

AT-T +T +-T. 

XX rr r r 


1-4 


ORIGINAL PAGE 18 
OF POOR QUALITY 



Boule through 


ORtGINAL PAGE 19 
OF POOR QUALiry 







Figure 1-2 FZ Model 
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Problem Pl-1 Compuco f(x) such that the solutioa T(x,r) of 

AT-P^|j,x< P,0<r< 

TU.l) - £(x), X < 0, 

T(0,r) • A(r'* , 0 < t < I 

also satlaflm* the boundary condition (1.2.1) 



The const..nt Pg Is the solid Peclet number [2Cj and from a practical 
viewpoint, the function A(r) in Equations (1,2.1) is the material melting point. 
Moreover, since numerical methods will be employed. Condition (1.2.1) will only 
be satisfied approximately in practice. 

Having stated the question for the lower solid region, the 
corresponding question for the upper solid region is analogous. Namely, let 
BCr)*^ be the required thermal gradient in the upper solid region at the upper 
solid-melt interface (x^), l.e.. 


3T 

ax 


(x.r) 


B(r), x-Q,0<r<l 


(1.2.3) 


Then find a surface temperature distribution g(x), x>p (henceforth also called a 
surface control function) , to be maintained such that the resulting temperature 
distribution, T(x,r), for the upper solid region satisfies (1.2.3). This is 
concisely stated in Problem Pl-2. 

Problem Pl-2. Determine g(x) such that the solution T(x,r) of 


AT-P^^,x> Q,0<r<l 

T(x,l) - g(x), X > Q, 

T(Q,r) - A(r'' , 0 < r < I 

also satisfies the boundary condition (1.2.3). 


(1.2.4) 


f To help make various computer codes listed in Appendix C easier to follow, 
the thermal gradients in Problem Pl-1 and Pl-2 are both represented by the same 
symbol, B(r); however, these gradients are not necessarily the' same. For 
generality, a similar remark holds l>>r the symbol A(r). 
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As before, In pracclce A(r) la set to the melting temperature and Equation 
(1.2.3) will only be satisfied approximately due to the numerical solution of 
the problem. 

Problems Pl-1 and PI-2. stated above belong to the class of so called 
Ill-posed or over-under posed problems. Unlike most classical second order 
boundary value problems where each portion of the boundary surface Is assigned a 
boundary condition, Problems Pl-l and Pl-2 have two boundary conditions 
(over-posed) assigned to each of their respective solid-melt interfaces (for 
example, in Problem Pl-l, T(0,r)« (r) and Tjj(0,r)«B(r)) and no boundary 
condition (under-posed) assigned to the lateral surfaces of either of the solid 
regions. Indeed, part of the problem is to determine the proper missing 
boundary condition (for example, T(x,l)-f(x) for Problem Pl-l) so as to relax 
the overposing of boundary conditions at the solid-melt interfaces. The 
solutions of Problem Pi-1 and Pl-2 are the subject of Chapter 3. 

Next suppose that the solid regions' surface temperature distributions 
f(x) and g(x) (see Figure 1-2, Equations (FZ6) and (FZIO)) are fixed (by design 
or happenstance) . Since the temperature at both of the solid-melt Interfaces Is 
assumed to be the melting temperature for FZ applications, the temperature 
distributions In both of the solid regions are computable (see Section 2.3). 
Hence the axial thermal gradients In the solid regions at the solid-melt 
Interfaces are computable. Thus, the axial thermal gradients in the melt zone 
at the solid-melt Interfaces (x*0 and x-Q) are known after invoking Equations 
(FZ2) and (FZ4) of Figure 1-2 and are denoted by 


Tj^(O.r) - A(r), 0 < r < 1 
Tjj(Q.r) - B(r), 0 < r < 1 


(1.2.5) 


The problem Is to determine a surface temperature h(x), 0 <_ x < Q (henceforth 
called the melt zone surface control function), to be maintained on the melt 
zone surface such that the resulting temperature distribution, T(x,r), for the 
melt zone satisfies (1.2.5). This Is concisely stated in Problem Pl-3. 

Problem Pl-3 Determine h(x) such that the solution T(x,r) of 


- Pj^ , 0 < X < Q. 0 
TCx.l) > h(x) , 0 < X < Q 

1(0, r) - C(r), 0 < r < 1 

I(Q,r) • D(r), 0 < r < 1 


< r < 1 



( 1 . 2 . 6 ) 


also satisfies the boundary conditions (1.2.5) 
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The coaetant P» is the liquid Peclec nuaber and from a FZ point of 
vl«w, C(r)and D(r) equal the material melting point. As with Problems Pl-I and 
Pl-Z. the numerical nature of the proposed solution method (the subject of 
Section 4.0} means that Conditions (1.2.5) will only be apprc'^-'iiately satisfied. 
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Before turning to chose moral and 
mental aspects of the matter which 
present the greatest difficulties, 
let the inquirer begin by mastering 
more elementary problems. 

— Sherlock Holmes, , "A Study in 
Scarlet" 

2.1 DESCRIPTION OF THE CLASSICAL PROBLEMS AND CHAPTER ODTLINE 

Before developing methods to compute the melt zone and solid regions' 
surface control functions which will yield the desired flat solid-melt 
interfaces'^, two more elementary problems must be dispatched. These are: 

Problem P2-1: Given a surface temperature distribution foi the melt zone, 

compute the resulting interior temperature distribution of the melt zone. 

Problem P2-2: Given a surface temperature distribution for one of the 

semi-infinite solid regions, compute the resulting interior temperature 

distribution for chat region. 

In addition to solving Problems P2-1 and P2-2, methods for 

approximating the interface gradients are presented in this chapter. The 
techniques developed to solve Problems P2-1 and P2-2 will have three important 
functions in this study. First, they will be used to generate the solid and 
melt zone gradients required at Che interfaces. Second, and - probably most 
important, the solution techniques for Problems P2-’. and P2-2 will introduce the 
essential definitions and dual integral transforms which will be used later to 
compute Che desired surface control functions (Chapters 3 and 4). Third, these 
techniques will be used to study how well (or poorly) the computed melt zone (or 
solid region) surface control function performs. 

Problems P2-1 and P2-2 are resolved in Sections 2.2 and 2.3 

respectively. Some numerical test cases are discussed in Section 2.3 along with 
two examples with correspondingly important remarks. 

2.2 SOLOTION OF PROBLEM P2-1 

Suppose Che melt zone of Figure 1-2 is isolated (and perhaps 
translated) as displayed in Figure 2-1. 


T To reduce the terminology, the solid-melt interfaces will henceforth be 
referred to merely as the interfaces. The axial thermal gradient in a solid 
region (or melt zone) at an interface will be referred to as a solid region (a 
melt zone) Interface gradient. 
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T(Xjj.r) - B(r) 

REMARK; FOR FZ 
PROBLEM, A(r) - 
B(r)-MELT PT. 

T(x^,r) - A(r) 
o 

Figure 2-1 Generalized Melt Zone 



Realistically, the melt zone end temperatures A(r) and B(r) are both 
the material melting temperature; however, for sake of illustration, we require 
only that A(r) and B(r) be sufficiently smooth. Problem P2-1 can then be 


mathematically stated as: 

Problem P2-3 : Determine T(x,r) such that 

AT - PT^ , 0 < r < 1 and < x < Xjj (2.2.1) 

T(x ,r) - A(r) , 0 < r < 1 (2.2.2) 

o 

T(Xjj,r) - B(r) , 0 < r < 1 (2.2.3) 

T(x,l) - h(x) , x^ < X < Xjj (2.2.4) 

and 

T^(x,0) - 0 , x^< * < Xjj (2.2.5) 


where A(r), B(r) and h(x) are sufficiently smooth, A( 1 )»h(Xg)and B(l)«h(xjj), and 
P (the Peclet number with the subscript "1” suppressed for convenience) is a 
positive constant. 


Before solving Problem P2-3, some notation is in order: 
Notation N2-1; 

(i) /^(r) - A(r)-A(l) 
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(ii) y(r) - B(r) - B(l) 
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(Hi) » JofA^r) where < A < A. <... 

is Che increasing sequence or rear roots of 

Che Bessel function . 

0 

(Iv) G(x) - Ph'(x) - h"(x) 


Solution Technique; The basic idea is to assume the solution T(x,r) is the sum 
of the lateral surface temperature h(x) plus some unknown function 6(x,r), i.e., 

T(x,r) - 9(x,r) + h(x) 


Problem P 

2-3 can then be 

recast 

as : 




A9 ■ P9 + G , 

X 

0 

< 

r < I 

and X < X < x_ 
o N 

(2.2.6) 


9(x^,r) 


0 

< r < 

1 

(2.2.7) 


9(x^,r) -5(r) 

> 

0 

< r- < 

1 

(2.2.8) 


9(x,l) - 0 , 

X 

0 

< 

X < X 

N 


(2.2.9) 

and 








9^(x,0) - 0 . 


< 

X < X 
• N 


(2.2.10) 

Although 

Equation (2.2. 

6) 

is 

more 

complex Chan Equation 

(2.2.1), Che 


corresponding boundary conditions are greatly simplified. First, the Dlrlchlet 
condition (2.2.4) is replaced by a simple homogenous boundary condition (2.2.9). 
In addition, because /#:i) - 5(1 ) ■ 0, the boundary conditions (2.2.7) and 
(2.2.8) can be further simplified by various Bessel series expansions. For the 
moment, assume 9(x,r) is expanded as 


e(x.r) - 


n«l 


( 2 . 2 . 11 ) 


Then using Che following well known property of Bessel functions [Ifl , 


/ 


<l'„(r)i|/j^(r)rdr - 


0 If a a 


** -^1*0 m) If a - « 


( 2 . 2 . 12 ) 


the functions of Equation (2.2.11) are computed Co be 
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C^(x) 


I 




e(x,r)4ij^(r)rdr 


0 


(2.2.13) 


If the Integral portion of Equation (2.2.13) Is denoted bydgCx), then Equations 
(2.2.11) and (2.2.13) may be combined to fora a dual Integral transform pair: 


e(x.r) 


2t|.^(t)0^'x) 


I 

and 9^(x) - / 0(x,r)i|;^(r)rdr . 
0 


(2.2.14) 


Unfortunately, the desired 9(x,r) of Equation (2.2.14) involves 6 q(x) which in 
turn requires knowing 9(x,r); fortunately, this rather circular problem may be 
resolved by Invoking Green's theorem. If both sides of the partial differential 
equation (2.2.6) are multiplied by Tj)n(r)rdr and the resulting terms integrated 
from r“0 to r»l, a application of Green's theorem combined with the fact that 




r-l 


r-0 


implies 


- ?r(x) - \ - G^(X) , x^<x<Xjj 


where 


G^(x) 

n 


•J" G(x)i|<^(r)rdr ■ G(x) 




(2.2.15) 


(2.2.16) 


Since /#(1) - 1) *0, the smooth functions r) and ff( r) may be represented by 

the following Bessel expansions: 

/^<r) - Zj4J(X r) (2.2.17) 

ff(r) - isjav) (2.2.18) 

a-l “ ° “ 


The above coef f Iclents and could be computed using Integral 
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represcncaclons 0,8). for example, 

1 

\ //#<rti,(X.r)rdr 

1 • 0 

Hc^«ver, to avoid the eventually required ntimerical integration of such 
If presentations, the coefficients and nay be approximated using a least 
squares method as described at the end of this section. Combining Equations 
( 1 . 2 . 7 ), (2.2.8), (2.2.12), and (2.2. I4)-(2.2. 18) , d^^x) may be uncoupled from 

0<K,r) as the solution of the following two point boundary value problem: 



n o 



6 

n 





G 

n 



(2.2.19) 


2 

Since Xjj > 0 , it la well known [ 5 J that Problem (2.2.19) has a unique 
solution. Although the solution of (2.2.19) could be determined by a variation 

of parameters method [ 1 ] , such a technique inevitably requires numerical 
integration. A more straightforward method is to discretize (2.2.19) in the 
following fashion. Pirst, the interval from Xp to Xjj is partitioned by the grid 
loints: 

4 


tj ■ Xp + jAx , M 

where M»Ax " solve the following finite difference analog of 

the boundary value problem (2.2.19): 




(ix)' 


-L-U- . J-1.--.M-1 



( 2 . 2 . 20 ) 



/ 
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The linear system (2.2.20) is cridlagonal anrf auaranteed to have a solution [9] 
if P*AxS.2. MorMver, the solution vector |y , •••, provides a second order 
approximation of (x), l.e., 


I“j - 


In addition, the boundary derivatives of d^may be accurately approximated [ 3 ] 
by the following unbalanced finite differences: 


'S'ix )a (-3y, + 16m, - 36m, + 48m, - 25m )/12Ax 
n . o ■ 4 3 ^ 1 o 

-*«vi + 


Since 


T^(x.r) 


OO 


h(x) 





9 (X) 
n 


(2.2.22) 


Equations (2.2.21) and (2.2.22) may be combined to approximate the axial 
gradients at x«Xq and x^^ (a very Important requirement in Chapters 3 and 4) . 


To finish this section, a short description is given of how the 
coefficients of Equation (2.2.17) are approximated (the same 

technique applies to Equation (2.2.18)). First, denote r^ » (i-l)/M, i ■ I, 

..., M + I and select N « M (typical^ N » 20 and M ■ 100). Define an(M+l) 
by N array L and (>(fl) dimension vector b by the respective elements: 


^IJ “ ^ ^1 


Let a be the solution of the linear least 


squares problem [l7. 


Chapter S]: 


La ■ b 


(2.2.23) 


Then the first N coefficients, of (2.2.17) are approximated by 




N 
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2.3 SOLUTION OF PROBLEM P2-2 

Analogous to the solution technique of Problem P2-1 In Section 2.2 
suppose the lower solid region of Figure 1-2 Is isolated as shown In Figure 2-2 


T(x,l) 




T(«,r) - A(r) 


STATE EQUATION 


at - P 


3T 

s 9x 


REMARK; FOR FZ PROBLEM 
A(r) - MELT PT. 


Figure 2-2 Generalized Lower Solid Region 
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Realistically, the upper end temperature A(r) of the lower solid ra;'-- the 

material melting temperature; however, for the sake of illustration, it is only 
required that A(r) be sufficiently smooth. In addition, it is assumeu i-u;: the 

lateral surface temperature f(x) is smooth, asymptotically constant as x - <» 
and is such that f' and f" approach zero as x » (loosely, this means f(x) 
resembles a horizontal line as x approaches - ■•) . Mathematically, the lower 
solid region case of Problem P2-2 may be stated as: 

Problem P2-4; Determine T(x,r) such that 

AT»PT, 0<r< 1 and x < 0 
x’ 

T(0,r) - A(r), 0 < r < 1 
T(x,l) - f(x) , X < 0 

Tj.(x,0) - 0 , X < 0 

and 

lim max I f (x)-T(x,r) 1 - 0 
x-*-o> 0<r<l 

The functions A(r) and f(x) are assumed sufficiently smooth, and for 

“compatibility, A(l) - f(0). In addition, ^li^f(x) exists and is finite and 

both f' and f " approach zero as x -*• - The constant P is assumed to be 

positive (the subscript "s” is suppressed for convenience). 

The notation e)stabllshed in Section 2.2 will be retained with the 
exception of G(x) which now represents G(x) ■ Pf'(x) - f (x). The solution 
technique is very similar to that used in Section 2.2 First, T(x,r) is expressed 
as 


(2.3.1) 

(2.3.2) 

(2.3.3) 

(2.3.4) 

(2.3.5) 


page 5a 
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T(x,r) - 0(x,r) + f(x) 

and Equations (2.3.1) - (2.3.5) are recast as: 


Ae-P6 +G,0<r<l and x < 0 

X 

(2.3.6) 

8(0, r) •/4(r) , 0 < r < 1 

(2.3.7) 

9(x,l) • 0 , X < 0 

(2.3.8) 

e^(x.O) - 0 , X < 0 

(2.3.9) 

lim max 1 6 (x,r) 1 *0 

3 C+-ao 0<r<l 

(2.3.10) 
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The boundary value problem (BVP) given by Che Equations (2.3.6) 

(2.3.10) Is solved in a manner similar to the solution of the BVP (2.2.6) - 

(2.2.10) . If/Ar) is represented as In Equation (2.2.17), then the BVP (2.3.6) 

- (2.3.10) may be transformed by the dual Integral transform pair (2.2.14) into 

the following boundary value problem: 

2 


6 


P6 ' - A 9 
n n n 


G , X < 0 
n 


6 .( 0 ) - — 




and 


lim 9 (x) « 0 


X 


where is still defined as in (2.2.16). Using a variation of parameters 
technique, the solution of this BVP is given by 


where 




_ , r* -a c _ a 


■ tS’ 


_ -Be _ 3 X 

G^e ° dtl e 


(2.3.11) 


S - ^ 
n ’ 


+ 4A^„ 


a =« (P + S )/2 
u n 


5 - (P - S )/2 
n n 

, r° ^ -St 

3 - \ G e “ ( 

n J n 


and 


- B + /# 
n n 2 




Since each summand in (2.3.11) is the product o^ an exponentially exploding and 
exponentially decaying term, the proof that 11m 9 (x) * 0 Is rather delicate 

and is reserved for Appendix B. The solution of Problem P2-4 is 
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i> 





— (r) - 

T(x,r) - f(x) + 9^(x) 
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(2.3.12) 


Since the float zone process also Involves the upper solid region of 
Figure 1-2, an upper region analog of Problem P2-4 must be solved. After 
translating the upper interface to 3^0 for convenience, the mathematical 
statement of such a problem is: 


Problem P2-5: Determine T(x,r) such that 


AT • PT^ , 0 < r < 1 , X > 0 

(2.3.13) 

T(0,r) - A(r) , 0 < r < 1 

(2.3.14) 

T(x,l) - g (x) , X > 0 

(2.3.15) 

T^(x,0) - 0 , X > 0 

(2.3.16) 

lim max lg(x) - T(x,r)( - 0 

(2.3.17) 


y+oo 0<r<l 


As before, A(r) and g(x) are assumed sufficiently smooth and, for 
compatibility, A(l) - g(0). In addition, lim g(x) exists and is finite, both 
g' and g" approach zero as jc**, and P>0. 


Without belaboring the details, the solution of Problem P2-5 is given 
by Equation (2.3.12) (g(x) obviously replaces f(x)) where &u(x) is still 
represented by Equation (2.3.11) with the G ,S ,o 6^^ unchanged but with 
new Aq and B^, namely n n n 



dc 


and 


-*n * 


As a computational aside, the numerical method described in Section 
2.2 can be used to approximate the analytically defined solutions of this 
section. For example, in the upper solid region case, if the surface 

temperature g(x) is rather constant for x, say, greater than some L, then the 
solution of Problem P2-5 may be approximated for 0 < x < Xjj by the solution of 

Problem P2-3 with Xjj set to, say 3L, and B(r) - g(xnj) and h(x) - g(x). 

Moreover, the gradient at the translated bottom, x - 0, of‘ the upper solid 

region may be accurately estimated by the approximate gradient generated by 
Equations (2.2.21) and (2.2.22). 
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Numerous test cases to numerically verify the above remarks were 
generated with A(r) ■ 0 (to simulate a solid-melt Interface), 0.1 < P <,1.0, and 
• 2L and 3L (L typically on the order of 10). The approximate temperatures 
(and gradients) so obtained were quite accurate. 

2.4 AN UNSETTLING FACT WITH AN ADDED NICE SURPRISE 

Consider for the moment Problem P2-4. If A(r) « 0 (to simulate a 

solid-melt interface) and two temperature distributions T^ and Tj are generated 
corresponding to two surface temperature conditions (x) and \ then If 
f^ ss fjt ic is reasonable to expect T^^ « Tj . In addition, If fj s fj near x 
« 0, then It Is also reasonable to expect that the corresponding thermal 
gradients of Tj and T^ will be close at x - 0. These Intuitive observations are 
Indeed^ true and may be rigorously proven after such concepts as "close" aie 
precisely defined. All of this, however, might lead to the assumption that If 
fl and fj are not close, then the corresponding thermal gradients ana 
temperature distributions near the simulated solid-melt interface (x » 0) are 
probably not close. This, of course, is not always true, and will be 
illustrated in this section by two examples. In fact, the second example will 
demonstrate the somewhat unsettling fact that It Is quite possible for f^(x) to 
exponentially explode while fjCx) remains nicely bounded with the corresponding 
thermal gradients at x •• 0 virtually Indistinguishable. In light of the 
development presented in Section 1.2, this Implies there might exist many varied 
surface control functions, all of which provide the required (or nearly so) 
thermal gradient at the desired Interface. If this is the case, then the float 
zone furnace designer may have at his disposal many different prospective 
surface control functions to choose from (a nice surprise). For example, the 
designer might select a surface control function that requires a minimum of 
power. 


The two examples in this section clearly demonstrate that small 
changes In the thermal gradient at the end boundary (x ■ 0) can result in a 
rather large change In the resulting surface control function. This, as noted 
before, can provide an entire family of useful surface control functions if the 
FZ designer is willing to permit a slight "misfit" (albeit small) between the 
desired and obtained thermal gradients at the end boundary (x ■ 0 for the 
following examples). Unfortunately, this also means that an attempt to measure 
the sensitivity of the required surface control functions to changes in the 
material or system parameters (which obviously produce changes In the desired 
interface thermal gradient) can be quite misleading and should probably not be 
attempted. 

Example E2-1; In this example, P - 0.1, A(r) - 0 and the lower solid region 
case (x<0) is selected. The nominal surface temperature fj is illustrated in 
Figure 2-3; the surface temperatures f^ , ..., fj (also illustrated in Figure 2-3) 
are perturbations of the nominal f^ . 

Letting T^ denote the thermal distributions corresponding the 
surface temperatures f^, the relative difference (measured in both L^ andL°° 
norma between the nominal gradient of T^ and each of the gradients of , i ■ 
2, .., 5, at the end boundary, x ■ 0, is Illustrated in Figure 2-4, 


For a function h(r), 0<r;^l, the L^ and L“ norms are (respectively) 

and )|h]j - max jh(r) j . 

« 0<r<l 


Ih 


'=•[1 


h (r)dr 
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Relative Difference Between 
Noalnal and Perturbed Gradient 


Figure 2-3 Notninel and Perturbed 
Surface Teoperacurea 



Flguve 2-4 Influence of Perturbations of Che 
Surface Tcnipuracure on the Thermal 
Gradient 



Noce chac even for Che cases where Che nominal surface 
cemperaCure, Is percurbed relaclvely close Co Che end boundary (x - 0), che 
corresponding per(urbaclons of che chermal gradlencs ac che x ■ 0 boundary are 
sclll remarkably close Co ChaC of che nominal gradlenC. 

Example E2-2; In Chls example, P • 0.1, A(r) ■ 0, and Che upper solid region 
(cranslaced Co x > 0) Is selecCed. Suppose Ic Is required chac che chermal 
gradlenC aC x ■ 0 be IdenClcally 1, l.e., T^(0,r) *1. A parclcular surface 
concrol funcclon gj(x) which will give Che desired resulc Is Che exponenclally 
growing surface CemperaCure: 


gj(x) - -10 + 10EXP(x/10), x>0 

In face, che corresponding Chermal dlsCribution T- is IdenClcall; jial Co g.. 
Suppose goCx), ..., g^(x) are surface conCrol funcclons ChaC eq . -(x) on an 
InCerval [0,Zj], 1 - 2,'..., 5 buC are asympCoClcally consCanC as .. ^iows (see 
Figure 2-5.) “ The relaClve and L* differences beeween Che Chermal gradlencs 
ac X * 0 of Che corresponding CemperaCure dlscrlbuclons ?2 , .... and che 
chermal gradlenC of T, Is lllusCraCed In Figure 2-6. Noce chac even when g^ 
(a bounded surface CemperaCure) separaces from g^ (an unbounded surface 
CemperaCure) racher close Co Che end boundary (x*0), che cwo corresponding 
chermal gradlencs ac x > 0 are remarkably close (see Figure 2-6, z > O.S). 
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It is an old maxim of mine that when 
you have excluded the impossible, 
whatever remains, however improbable, 
must be the truth. 


The main thrust of this chapter is to provide solution methods for 
Problems PI— 1 and Pl-2. Beginning with Problem Pl-1 , suppose a temperature 
distribution T(x,r) is required to satisfy two known boundary conditions 

T(0,r) - A(r) (3.0.1) 

T^(0,r) - B(r) (3.0.2) 

at the lower solid region's end boundary (the melt-solid interface in practice) 
as depicted in Figure 3-1. 


—Sherlock Holmes, "The Adventure 
of the Beiryl Coronet" 

I 



AX ■ P^ -|“ (State Equation) 


T(0,r) - A(r) 
|^T(0,r) - B(r) 
T(x,l) - f(x) 



KNOwN 


UNKNOWN 


”£.MARK: f^or FZ problems, 

A(r) ’iMel ting Temperature 


Figure 3-1 


FZ Lower Solid Region Problem 
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Recall from Section 1.2 that the problem is to find some (at this point unknovm) 
cooling control function f(x), x<0, such that the solution of the veil-posed 
boundary value problem: 


AT - , x< 0 . 0 < r < 1 


(3.0.3) 


T(0,r)-A(r) , 0 < r < 1 


(3.0.4) 


and 


T(x»l)-f(x) , X < 0 (3.0.5) 

also satisfies the addition boundary condition (3.0.2). The basic idea of the 
proposed method is to solve the boundary value problem (3.0.3)-(3.0.5) by the 
method described in Section 2.3 and, in the process find a sufficient number of 
conditions to allow the calculation of the desired, but unknown, f(x). First, 
from a practical point of view, any viable control function f(x) should become 
rather constant as the distance from the lower interface increases. Thus 
it is expected that: 


11a f(x) exists and is finite 

X -» - BO 

and 

11a f' (x) - 11a f"(x) - 0 

yx an TC*- « 

Proceeding as in Section 2.3, denote^ 

T(x,r) - 8(x,r) + f(x) 

G(x) - Pf* - f'* 

/#(r) - A(r) - f(0) 

^(r) - B(r) - f' (0) 

For f(x) to be compatible with A(r) and B(r), 

A(l) - fvO) 
and 

B(l) - f’ (0) 

Then Equations (3.0.2)-(3.0.5) reduce to 


(3.0.6) 


(3.0.7) 


(3.0.8) 


(3.0.9) 


+ For Che moment, suppress the solid subscript "s", i.e., 


P -P. 
s 
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1 


A9 • P0^ + Gfx) 
9(0,r) •/4(t) 
9^(0,r) -?(r) 

9 (x,l) - 0 


J 
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(3.0.10) 


Denote i(> (r) " J (X r) where Xj < X, < X^ < ••• are the real roots of the 
a* o n 

zero order Beasel function. As in Chapter 2, let 




ana 


0 (x,r) 


9j,(x) 


~2 

M-1 Xj^) 


I 9(x,r)tjjj^ 

Jq 


(r)rdr 


(3.0.11) 


(3.0.12) 


form a dual Integral transform pair. If/i^r) and ^r) are expanded In the 
Bessel series: 


^ 7 




and denoting 


(3.0.13) 


n-1 


S(r) • J (A r) 

non 


(3.0.14) 


n"l 


■ I 


Gjj(x) - I G(x) vj;^(r)rdr ■ G(x) ^ 


>^1 (Xj^) 


M 


then operating on (3.0.10) by the Integral transform (3.0.12) yields 


(3.0.15) 


'll 2 


(3.0.16) 


(3.0.17) 
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(3.0.18) 


In practice, the Bessel coefficients, and in (3.0.13) and (3.0.14) are 
approximated by least square methods as described In Section 2.2. In addition, 
for FZ applications > 0 because A(r) is constant (the material melting 
temperature). Denoting 




- (P + S„)/2 


. (P - S„)/2. 


the solution of (3.0. l6)-(3. 0.18) is then: 


* V f - 

^ Jn ^ ^ 




(“H 


(3.0.19) 


GnCt) 


Since G|^(x) ■ (Pf’ (x) - f*'(x) ) (a^Y 2 approaches zero as x proceeds toward 
negative infinity (see (3.0.7)) ,* an" argument similar to that found in Appendix B 
will show that the second summand in (3.0.19) approaches zero as x approaches 
negative infinity; since is positive, the first summand of (3.0.19) shares a 
similar fate. In light of (3.0.6), it is reasonable to assume (or require 
depending on the point of view) that 

11m max I T(x,r)-f (x) ( - 0 
0<r<l 


. ' A , 
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aad hence lim Qn(^) ~ 0. Combining these observations with (3.0.19), the 
x-*^ H 


remaining conditions to be used in determining f(x) are easily discerned, 
namely: 


lim 




G„(t) -8 t 

e dt 


e - 0 (3.0.20) 


Since <0, an analysis similar to that of Appendix B will show that (3.0.20) 


will he satisfied if 


-Svt 

- -J 

^ Ai 


(3.0.21) 


Since 


Gjj(t) - G(t)Jj(Xjj)A^ - (Pf' (t) - f'(t))J^(Xjj)Aj^, 


combining (3.0.6)-(3.0.9), and (3.0.21) with two applications of integration by 
parts yields: 






-8 t 

^(^-P) I f(t)e “ dt 


(8j^-P) A(l) + B(l) 


Denoting 


M 


8 (8 -P) 


(P-^)A(D-B(l) 


m m 


the desired properties of the surface function f(x) may be summarized as: 
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f(0) - A(l) and ' f (0) - B(l) 
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lim £(x) exists and Is finite 


and f' (x) and f” (x) 0 as x -« 


(3.0.22) 


■Si- 


i: 


-6 t 

f(t)e ^ dt , M 


To numerically approximate such a surface control function as f(x), let 


NSYS 


(k-l)t 




(3.0.23) 


Then, In light of (3.0.22), set 

=1 + “SlSTS ■ 




B(l) 


and 


NSYS 

2 -4 

k-l -- 


Exp((k-l-6jj)t)dt - Rj, , M-1,2,..'. ,MIERM 


If Che (MTERfH-2) by NSYS matrix L and (MTERM+2) dimension vector b are defined, 
for j - 1,2,..., NSYS, by 


- 1 and bj^ - A(l) 


■" • 


(3.0.24) 


‘’i - V2 


i - 3,4, . . . , >rTERM +2 




t Tne index in the expansion of f(x) starts at k»l instead of k-0 to make 
referencing this section from the accompanying FORTRAN documentation easier 
(Appendix A). The Index limits NSYS and MTE8M noted here will be used in the 
same role in the accompanying FORTRAN codes (Appendix C). 
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then, provided MTERM4-2 > NSYS, the coefficients Cj^ of (3.0.23) may be set to the 
least squares solution of 


(3.0.25) 


- C 


The solution method for Problem Pl-2 is similar to the above and hence 
most of the details are left tc the reader. Using the notation established for 
Problem Pl-2 denote 



G(x) - Pg (x) - g" (x) 

/#(r) - A(r) - A(l) 
and 

arcr) - B(r) - B(l) 

As before, expand /#[r) and^(r) as 

q>_ 

- 

n«l 

Ci 

- ^ J (X r) 

/ ^ o n 

n»l 

Then using the above established notation for S^, a^j , and Sn» the desirable 
properties of an upper solid region surface control function g(x) are summarized 
as: 


and 

^(r) 


g(Q) - A(l) 
g'(Q) - B(l) 

a (Q-t) 

G(c)e dt 


lim g(x) exists and is finite 
x^ 


X J,(X ) f 

4-^ -»a) • J 


(3.0.26) 

(3.0.27) 

(3.0.28) 


(3.0.29) 
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lia g' (x) ■ lim g*fx) • 0 (3.0.30) 

X^ X-H» 


Mimicking the previous analysis, (3.0.28) is simplified by two applications of 
integration by parts. Thereafter, an approximation of g(x) given by: 

MSYS 

g(x) * (3.0.31) 

kSl 


is substi^uted into Equations (3.0.26)-(3.0.28) and the desired coefficients 
determined by a least squares method. 
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4.0 THE HEATING CONTROL FUNCTION 


In five minutes you will say that it 
is all so absurdly simple. 

— Sherlock Holmes, "The Adventure 
of the Dancing Man" 


The ultimate goal of this chapter is the solution of Problem Pl-3. 
Suppose that the temperature distribution T(x,r) for the melt tone is required 
to not only satisfy the state eqxiation 


AT-P^T , 0<x<Q ,0<r<l 


(4.0.1) 


but also must satisfy for 0<r<l the four boundary conditions: 

TC0,r) - C(r) 

T(Q,r) - D(r) 

T^(O.r) - A(r) 

and 

T^(Q,r) - B(r) 

Infortanately, not only is Coo much information supplied for the two end ‘ 
boundaries (:^0 and x“Q; see Figure 4-1), no Information whatsoever is supplied ' 
' r the remaining boundary, r • I (again see Figure 4-1). 


(4.0.2) 

(4.0.3) 

(4.0.4) 

(4.0.5) 


T(Q,r)-D(r) 

Tx(Q,r)-B(r) 


T(0,r)-C(r) 

"jj(0,r)-A(r) 


ST.ATE- EC3UATI0N 




I 


Figure 4-1 


FZ Melt Zone Problem 
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The problem, Cheref ore ,1s Co find some heating concrol function h(x), 
0<x^Q, such Chat Che solution T(x,r) of Che well-posed problem defined by the 
boundary condition T(x,l) - h(x), the boundary conditions (4.0.2) and (4.0.3) 
and Che statj equation (4.0.1) also satisfies (or nearly so) the additional 
conditions (4.0.4) and (4.0.5). 

For simplicity, the functions C(r) and D(r) are both assumed to be 
zero*. The generalization for nonconstant C(r) or D(r) is similar to the 
following analysis and is left to Che interested reader. As in Chapters 2 and 
3, define 


G(x) 

/^(r) 


P jh’ (x) -h” (x) 
A(r) - A(l) - 


//nJ (^„r) 
o n 


and 


ff(r) - B(r) - B(l) - 


00 


a r) 

no n 

n«l 


If T(x,r) is decomposed into 


(4.0.6) 

(4.0.7) 


(4.0.8) 


T(x,r) - 9(x,r) + h(x) 

Chen Equations (4.0.1)-(4.0.3) Imply* 

A0 - P0X + G , 0 < x< Q , 0 < r < 1 


9(0, r) - 9(Q,r) - 0 , 0 < r< 1 


Denoting G (x) ■ G(x)* Jj^ ( X^)/ transforming 
transform X3.0.12), 


r -ps^ -X ^ 

n n n n 


G^ , 0 < X < Q 


9„(0) -9 (Q) - 0 
n u 


(4.0.10) 


by the 


(4.0.9) 


(4.0.10) 


integral 


(4.0.11) 


t For FZ work, both C(r) and D(r) are set to the material melting temperature 
which can itself always be assigned to be zero on some translated temperature 
scale . 

*For convenience, suppress the "1" (liquid) subscript, i.e., P^»P 
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and 


o • (P + SO/2 


8^- (P - S^)/2. 


By a variation of parameters method, the solution of (4.0.11) is 






X 

.^n) “ 


(t)e " dt 


_^n^ 


-K + 
n 


iS "" 


-a t 

(t)e ^ dt 




(4.0.12) 


where 


S „(0 


r a c 

' ^ " 

) e 
n'- 


a (Q-c) 0„(Q-t) 

3 " -e " 


dt 


- e 




I 




For the desired h(x) to be- compatible with Equations (4.0.2)~(4.0.5) (recall 
C(r) and D(r) are set to zero), h(0) ■ 0 » h(Q), h (0) ■ A(l) and h (Q) * B(l). 
In addition, since 0^(0, r) -/#(r) and e^CQ.^) -^(r), 




(4.0.13) I 

. 1 ' 
i 
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Combtnlng Equations (4.0.13) with the derivative of (x) (obtaining by 
differentiating"^ (4.0.12)) yields 




-AU) 


( 

/ 


1 K^(n,t) h(t) dt 
*^0 




and 


(4.0.14) 


n n 1 n 


+ B(l) 



K2(n,t) h(t) dt 


(4.0.15) 


where, if denotes, 

^^n "i (P^-Sq)/C1-E*p(-6^Q)J 

then kernels and K 2 In Equations (4.0.14) and (4.0.15) are defined by 

(4.0.16) 


K. (n,t) - C 
i n 


-o t -(0 t + S Q) 
n n 0 ^ 

3 - e 


and 


■ir 


/ 

K^(n,t) - (l-e " ; 


S/.-Q) 


(4.0.17) 


t The actual process of differentiating (4.0.12) Is routine but laborious and 
is left to the Industrious reader. However, this is not to Inply that great 
care should not be taken; several of the Integrands are the difference of large 
functions (a numerically delicate situation). For the industrious reader 
willing to check these results, the removal of derivatives from Integrands by 
Integration by parts Is necessary. 


{ 

lb- 
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The desired h(x) is numerically approxlmaced by some expansion of the form: 

NSYS 

h(x) « 

k-1 

(for example, let h^(x)»x ). To finish this development, solve for the 

coefficients c. in Equation (4.0.18) by solving (in a least squares sense"^) the 
System (4.0.19)-(4.0.23) given below. 


where 


NSYS 


2 " ° 

k-1 

NSYS 

Z) °’A«> - “ 

k-1 

NSYS 

c^h^CQ) - BCD 

k-1 

NSYS 

k-1 

NSYS 

, n-1,2,---, mTERM, and j - 1,2 


k-1 


jnk 


f'’ 

• I Kj(n.t) h^^(t)dt 
JQ 




(DO. 19) 


(4.0.20) 


(4.0.21) 


(4.0.22) 


(4.0.23) 


and 




t In order to use a least squares method, the index limits NSYS and MTERM 
should be selected such that NSYS/2 < MTERM + 2. 


I 




1 
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And here — ah, now, this really is 
something a little recherche 

— Sherlock Holmes, "The Musgrave Ritual 


5.1 SOLID REGION SURFACE CONTROL FUNCTIONS 

In this section, the methodology developed In Chapter 3 is illustrated 
‘ using material and system data provided by NASA. Recall that the problem is to 
find, after being given the melt zone surface temperature distribution, the 
surface control functions for the solid regions which are compatible with flat 
interfaces. The material and system parameters used are listed in Table 5-1 and 
were provided (and in some cases appropriately modified) by E. Kern (NAiA 
contractor, [llJ ) and E. Cothran (NASA, [4]). Tha material selected was silicon. 


TABLE 5-1 MATERIAL AND STSTEM PARAMETERS 


PARAMETER 

VALUE 

Radius 

i 

0.2413 cm 

Malt Length 

1.1684 cm 1 

Conductivity 


Solid 

7.5 cal/°K tr. sec 

Melt 

16 cai/°K IT. sec 

Density 

, 3 

Solid 

2.28 gm/cm 

Melt 

2.53 gm/cm^ 

Heat Capacity 


Solid 

0.241 cai/°K gm 

Malt 

0.265 cal/°K gm j 

Latent Keat 

431 cal/gm j 

Growth Rate 

2.5 mm/min 

Feclet Number 

i 

Solid 

0.00734 ! 

Melt 

0.00<:*21 1 

Melting Temperature 

16930 R 


-he melt zone surface temperature distribution used was suggested by E. Ke’"i 
lU and is illustrated in Figure 5-1. 


The surface control functions for various combinations of MTERM ar : 
.'JSYS (see Equations (3.0.23) and (3.0.24)) obtained by the methods of Chapter "! 
are illustrated in Figure 5-2. 
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In light of Section 2.4, the variety of control functions depicted in 
Figure 5-2 is expected. In addition, the results of Section 2.4 suggest that 
the tuo lower solid surface control functions that eventually are above the 
melting temperature may be modified as Illustrated in Figure 5-3 without 
substantially changing the thermal gradients at x>0. 


The relative differences between the thermal gradients (in the solid regions) 
required^ at the Interfaces (x"0.0 cm and 1.1684 cm) and those obtained using 
the surface control functions defined in Figures 5-2 and 5-3 listed in 
Table 5-2 


TABLE 5-2 RELATIVE DIFFERENCES BETWEEN THE 

REQUIRED INTERFACE GRADIENTS AND THOSE 
RESULTING FROM THE USE OF THE SOLID 
REGIONS' SURFACE CONTROL FUNCTIONS 




Solid 

Surface Control 

Relative Difference 

MTERM 

MSYS 

Region 

Function Definition 

(in L'^’ norm) 


3 

Upper 

Figure 5-2 

0.0175 



Lower 

Figure 5-2 

0.17 



Upper 

Figure 5-2 

0.00012 

6 

6 

Lower 

Figure 5-2 . 

0.001 



Lower 

Figure 5-3 

0.056 



Upper 

Figure 5-2 

0.000027 

8 

9 

Lower 

Figure 5-2 

0.0013 



Lower 



Figure 5-3 

0.061 


As an aside, in a series of test cases over a range of values for 
MTERM and NSYS, the relative differences between Che required interface 
gradients and chose obtained using surface control functions first decreased and 
Chen increased as NSYS (or MTERM) was increased while holding fixed the value of 
MTERM (or NSYS). Naturally, this is to be expected since an approximate 
solution of an ill-posed problem is attempted by employing an overposed system. 
This, of course, reinforces the old maxim of always examining a computed 
solution for "reasonableness.” In fact, the computer software developed (see 
Appendix A) to determine the solutions of Problems Pl-1 and Pl-2 automatically 
computes the relative-errors between the required interface gradients and chose 
resulting from Che use of the surface control functions. It cannot be 
overstated: always examine .these relative errors before accepting a computed 
solution as reasonable. 


T See the Boundary Conditions (1.2.1) and (1.2.3). 
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The techniques developed in Chapter 4 are Illustrated in this section 
using material and system data as provided by NASA. The problem Is to determine 
a melt zone surface control function compatible with some solid regions' surface 
temperature distributions (provided a priori) such that flat melt-solid 
Interfaces are achieved. The material and system parameters used are listed In 
Table 5-3 and were provided by E. Kern (NASA contractor [11] and E. Cothran (NASA 
l4j). The material selected was silicon. 


TABLE 5-3 MATERIAL AND SYST^ PARAMETERS 


PARAMETER 

VALUE 

Crystal Radius 

0.69 cm 

Melt Length 

1.43 cm 

Conductivity 


Solid 

7.5 cal/°K m sec 

Melt 

16 cal/®K m sec 

Density 

3 

Solid 

2.28 gm/cm 

Melt 

2.53 gm/cm^ 

Heat Capacity 


Solid 

0.241 cal/°K gm 

Melt 

0.265 cal/°K gm 

Latent Heat 

431 cal/gm 

Growth Rate 

2 tnra/min 

Peclet Number 


Solid 

0.01685 

Melt 

0.009638 

Melting Temperature 

16930 K 


The lower and upper solid regions' surface temperature distributions 
used were suggested by E. Kern [1]J and are Illustrated in Figure 5-4. 


The melt zone surface control functions for various combinations of 
MTERM and NSYS (see Equations (4.0. l8)-(4. 0.23)) obtained by the methods of 
Chapter 4 are Illustrated in Figure 5-5. 

Because of the Ill-posed nature of the problem, a variety of surface 
control functions is expected. The relative differences between the required 
melt zone Interface gradients and those obtained using the surface control 
functions defined in Figure 5-5 are listed In Table 5-4. 


t 


See Boundary Condition (1.2.5) 










TABLE 5-4 RELATIVE DIFFERENCES BETWEEN THE REQUIRED 

INTERFACE GRADIENTS AND THOSE RESULTING FROM 
THE USE OF THE MELT ZONE SURFACE CONTROL FUNCTION 




MELT-SOLID 

RELATIVE DIFFERENCE 

MTERM 

NSYS 

INTERFACE 

(in L'^ norm) 

4 

3 

Upper 

.53 



Lower 

.24 

8 

'6 

Upper 

.28 



Lower 

.07 

14 

Q 

Upper 

.11 



Lower 

.13 

20 

12 

Upper 

.05 



Lower 

.04 


In a series of Cest cases over a range of values for MTERM and NSYS, 
Che relative differences between Che required Interface gradients and Chose 
obtained using the surface control functions first decreased and then Increased 
as NSYS (or MTERM) was Increased while fixing the value of MTERM (or NSYS). As 
in the previous section, this Is Co be expected because Che solution technique 
employed uses over-posed systems to approximately solve an ill-posed problem. 
As before, Che computed melt zone surface control function should be checked for 
reasonableness. For example. It Is quite possible chat the surface control 
function could be less Chan the melting temperature on part of Che melt surface 
In which case the control function should be modified or rejected. In addition, 
the relative differences between the required melt zone Interface gradients and 
chose obtained using a candidate melt zone surface control function should be 
examined before accepting the control function as an approximate solution of 
Problem Pl-3. Incidentally, these relative differences are approximated and 
displayed by Che software developed for Problem Pl-3. 
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6.0 FUTURE WORK AND UNRESOLVED ISSUES 

Although the work presented In this report Is, In Itself, rather 
complete, several side Issues remain unresolved and should be Included In any 
continuation of this type of research. In this chapter, some of these Issues 
are addressed. 

6.1 VERIFICATION USING FREE BOUNDARY ALGORITHMS 

The basic Idea of the three methods described In Chapters 3 and 4 was 
to determine the properties a surface control function must have If a flat 
Interface was to be maintained. Unfortunately, both methods Involved many 
numerical approximations and some simplifying assumptions. As an example, for 
the method described in Chapter 4, the thermal distributions in the assumed 
infinitely long solia regions were approximated by numerical methods designed 
for finite length regions. In addition, the interface gradients were 
approximated by finite differences (a second source of error) followed by least 
squares Bessel function fits of these approximate interface gradients (a third 
possible source of error.) Thereafter, the surface control function was 
approximated by solving an overposed system of equations using only a finite 
number of terms in the control function (another source of error). Given these 
several possible sources of error, the actual Interface shapes maintained using 
the computed surface control function should be constructed using some 
multiphase free boundary algorithm (for a survey, see H9l). The results of such 
numerical experiments should hopefully further verify the methods discussed in 
this report and should indicate some future areas to be studied with error 
reduction In mind. 

6.2 MAINTAINING CURVED INTERFACES 

Although thermal stresses, which can generate defects in the crystal, 
are generally minimal for a planar interface [2{j , a slightly curved interface 
shape is also quite desirable in some cases. Specifying the desired shapes, the 
required surface control functions could probably (with sufficient 
investigation) be constructed using methods similar to Chose in Chapters 3 and 4 
after the introduction of transformations similar to those described in [12j. 

6.3 NON-DIRICHLET BOUNDARY CONDITIONS 

Boundary conditions ocher than the Dirichlet type (Equations (FZ6), 
(FZ8), and (FZIO) of Figure 1-2) should be investigated. Fortunately, much of 
Che work for this type of problem will probably be straightforward. For 
example, suppose a question like Problem Pl-3 is to be solved where the 
Dirichlet boundary condition (see Equation (1,2.6)) 

T(x,l) = 0 •' X < Q 

is replaced by a boundary condition of Che type 

T^fx.l; = :<[r7x,l) - S°'(x;] Co. 3.1) 
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where S(x) Is the desired surface control function (for example, If a >4, then 
S(x) might be the temperature of a furnace wall providing radiant heating). To 
solve this problem, first solve Problem Pl'3 as stated In .Section 1.2. Using 
the computed Dlrlchlet type surface control function h(x), next solve Problem 
P2-3 (let Xg ■ 0 and x^ ■ Q) • Then place the resulting temperature 
distribution T(x,r) Into the boundary condition (6.3.1) and solve for the 
desired surface control function S(x). 


6.4 BASIS FUNCTIONS USED TO EXPAND THE CONTROL FUNCTIONS 

The exponentially decaying functions used to expand the solid regions' 
surface control functions (f(x) and g(x) of Equations (3.0.23) and 3.C.3I) 
respectively) were selected because they represented what a typical control 
function would be Intuitively expected to resemble and because they allowed for 
simple integrations In Equations (3.0.22) and (3.0.28). However, from a 
computational point of view, these basis functions are not the best"^ . For 
example, some preliminary experiments Indicate that replacing Equations (3.0.23) 
and (3.0.31) by 


and 


NSYS 

f(x) « c^ + ^ (-<t+k-2)^) 

k-2 


NSYS 

g(x)« EXP (-(c-fc+2+Q)^) 

k-2 


respectively can significantly reduce the least squares residuals of overposed 
systems like Equation (3.0.25). More study Is needed to find basis functions 
that both further reduce the least squares residuals and are not too difficult 
to Integrate In equations like (3.0.22) and (3.0.28). 

6.5 APPLICATIONS OF LINEAR PROGRAMMING 

The linear system of equations (3.0.25) wi^l be underposed If MTEPM+2 
< NSYS. However, the desired coefficient vector C may still be ^ete^mine^ as 
follows. Let T be the residual vector of Equation (3.0.25), l.e., r ■ b - L.C. 


t It Is sometimes dangerous to approximate a function f(x) by a sum: 

N 

k*l 

or most of the functions fjr(x) "resemble" each other, e.g., fj^(x)» 

This, for example, is why Chebyshev polynomials are preferred over the 
standard basis, for polynomial approximation on certain 


where all 
Exp(kx) . 
so-called 
domains . 
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Then solve the linear prograomlng problem 113k pg 15] 


min 


£lr^l 


subject to LC + r • b 


(6.5.1) 


In fact, such a technique might be used tq reduce the chance of say, f(x) of 
(3.0.23), becoming positive*^ (recall that the melting temperature was 
translated to zero) as x approaches negative Infinity. To accomplish this, 
first select a grid, 


partitioning a portion of the lower solid region, and then adjoin to (6,5.1) the 
additional N constraints (see Equation (3.0.23)): 


NSYS 

^ c^ EXP ((k-l)Xj^) < 0 , 1-1, ••• , N 
k-1 


Some preliminary numerical experiments suggest this Idea has sufficient 
potential to warrant further investigation. Although this discussion has 
centered on the lower solid region, these ideas are applicable to either of the 
solid regions or to the melt zone. 

6.6 MODELS AND REALITY 

The FZ process was modeled In this effort as a steady state proce 's on 
an infinitely long boule. Unfortunately for the modeler (but fortunatel / for 
the commercial FZ operator), the boule has finite length*. For finite length 
boules, the problem of finding the proper surface control functions to maintain 
flat Interfaces would now Involve end effects and various time transients. 


t The partial differential equations used to establish the desired surface 
control functions are quite Ignorant of the fact that surface control functions 
for solid regions should always be below the material melting point. In fact, 
In some numerical experiments where MTERM and MSYS were large, the computed 
surface control functions for one of the solid regions became greater than the 
melting temperature. This Is only one of the dangers In trying to solve an 
overposed problem. 

* Fortunately, the assumptions and results of this effort are still quite 
reasonable for long boules with slow growth rates. 
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However, the basic Ideas discussed In this report could probably be extended to 
cover such difficulties. The resulting partial differential equations would 
Involve the additional tern 


1^ T(:c,r,t) 


(where t represents time) and hence would be parabolic lustead of elliptic. The 
boundary conditions would also be time dependent. However, the dual integral 
transform pairs used In this report should still provide enough information 
concerning the surface control functions (required for flat Interfaces) to allow 
for their construction. 

In addition, the fluid dynamics of the melt zone should be 
Incorporated In the computation of the surface control functions. Of the topics 
discussed in this chapter, this Is undoubtedly the most difficult one to model 
and resolve . 
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(or a uaar guide co davelcned software) OF POOR QUALITY 

A. 1 INTRODUCTION 

The data Input procedures and output interpretations for the software 
developed to approximate the solutions of Problems PI-1, Pl-2, PI-3 and P2-1 is the 
subject of this appendix. To begin. Problem P2-I is covered in Appendix A. 2 
and is followed by Problems PI-1 and Pl-2 in Appendix A. 3. To finish, Problem 
Fl-3 is Che subject of Appendix A. 4. 

A. 2 USER CONSIDERATIONS FOR PROBLEM P2-1 SOFTWARE 

The data input procedure and output interpretation for the software 
developed for Problem P2-1 is Che subject of this section. To begin, all Che 
required data are input in the form of punched cards. The definitions and formats 
of this input data are sunnrrlzed in Table A-I. 


TABLE A-1 PROBLEM P2-1 SOFTWARE INPUT 


PROGRAM 

SYMBOL 

VARUBLE 

DEFINITION 

UNITS* 


READ (5, 16) IHFC ,M 
16 FORMATC2HO) 


IHFC 

M 

" 1 if a cubic spline will be used co approximate 
Che boundary function h(x) in Condition (2.2.i). 

■ 0 if the user will supply a functional form of h(x) 

(see Condition (2.2.4)). In this case, the user 
must insert this functional form of h(x) in the sub- 
routine HFC (see the software list in Appendix C.2). 

' 

Number of 'i«nocs used co approximate ;ux) (see Condi cron 

(2.2.4)) by a cubic spline (IHFC-i). If IHFC-0. sec 

vM»0. 



00 32 I - 1,.M 
READ(5,22).'a3(r),YD(I) 
32 CONTINUE 
22 FORMAT (4E20, 10) 

;o(i) 

The dimensionless axial position of Che Ith knot used 
to approximate h(x) (see Condition (2.2.4)). Ignore 
if IHFC-1. XDaX.'CDd+l) . 

x.d 


x,r and R will denote the axial distance, the radial position and the rod radius 
respectively. °T will denote whatever temperature scale the user prefers. 
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TABLE A-1 PROBLEM P2-1 SOFTWARE INPUT (CONT) 


PROGRAM 

SYMBOL 

7ARZABLE 

DEFINITION 

UNITS 

YD(I) 

The surface taaperacure represeatad by the Ich knot 
used CO approxisMiCt b(x) (see Coadlcioo (2.2.4)). 
Ignore If iMFO»0. 



READ (5 , lu) P ,X0 .XN.MSTJM.NGRID , NR 
10 FORMAT(3F10.5,4I10) 


P 

Peclet number 

dimen- 

sionless 

XO 

Dlmensioaless ajr'al position of bottom boundary 
displayed In Figure 2-1. 

Remark: The boundary cemperacure distrlbucion, 

A(r/R) at XO (see Condition (2.2.2)) Is 
user supplied in Subroutine AFC (see Appendix 
C.2). 

x/R 

XN 

Dlmeaslonless axial position of cop boundary 
displayed in Figure 2-1. 

Remark: The boundary cemperacure discribucion. 

B(r/R) at XN (see Condition (2.2.3)) is user > 
supplied in Subroutine BFC (see Appendix C.2). 

x/R 

MSU« 

The first MSDM terms of the expansion in Equations 
(2.2.14) are used to approximates (x,r). MSUM must 
be less chan 21. 


SGRID 

(XN-X0)/NGRID is the grid size employed in System 
(2.2.20). In addition, the final cemperacure dis- 
cribucicn is output for NGRUH-L axial values from XO 
CO XN. NGRID may not exceed SOO. 


NR 

The final cemperacure distribution is 

output for NR+1 radial values' (r/R) from 0 to 1. NR 

may not exceed 100. 



An inpuc saople is illuscraced nexc in Figure A-1. 

^ The output is labeled clearly for ease of use. The inpuc data are 
first viewed followed by the Chennai distribution (the approximate soluCi.on of 
Problem P2-3) given in cable format (see Figure A-2) . Incidentally, the 
axial and radial positions in Figure A-2 are given in dimensionless form 
(x/R and r/R) . The thermal gradients at XO and XN are given last 
in a cable format (again, see Figure A-2). 
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ORIGINAL PAGE IS 
OF POOR quality 

To flAlah, Che user ntusc supply an algorithm co fit (in che least 
squares sense) a linear coablxiacion of functions to a given set of data points 
(see the and of Section 2.2 for a short discussion). This algorithm is 
required in che subroutine COEFS. In addition, an algorithm to evaluate che 
Jq Bessel function (required In che subroutine JO) oust be provided. These 
required algorithms are generally available from the host computer libraryt or 
may be obtained from various sofCx»are packages such as che mSL and RJNPACX. 


• I 

* The user is warned, however, chat many hose computer mathematics libraries . • ■ \ 
(with che general exception of IBM) still contain numerous faux pas chat were ; 
well known years ago and still remain uncorrected. | 


A- 7 


i 

\ 



USER CONSIDERATIONS FOR PROBLDIS Pl-1 and Pl-2 SOFTWARE 


A. 3 


The data Inpuc procedure and output interpretation for the software 
developed for Problems Pl-1 and Pl-2 are the subjects of this section. 
Recall that the general problem is to find the solid regions’ surface control 
functions (f(x) and g(x) in Problems Pl-1 and Pl-2 respectively) which, 
for the sake of flat interfaces, are compatible with the a priori given melt 
zone surface temperature distribution (h(x)). Ac present, all the required 
data are inpuc in the form of punched cards. The definitions and formats 
of the inpuc data are summarized in Table A-2. 
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TABLE A-2 PROBLEMS Pl-1 AND Pl-2 SOFTWARE INPUT 


PROGSAM 

SYMBOL 

VARIABLE 

DEriNTTION 

units' 

READ(5.80)IHFC,M 
ao FORMAT (2110) 

IHTC 

- 1 if a cubic spline will be used co approxlmace che 
melc zone surface camperacure dlscribucion, h(x). 



■ 0 if Che user will supply a functional form of h(x) . 
In chis case, che user muse insert this functional 
form of h(x) in the subroutine HFC (see the software 
list in Appendix C.3) 



Number of knots used co apptoxima.ee h(x) by a cubic spline 
(IHFC-l). If IHFC-0, sec M-0. 



DO 32 I - l.M 
READ(5,22)XD(I),YD(I) 
32 CONTINUE 
22 FORMAT (2E20.10) 

•UjCI) 

The axial position of che Ich knot used co approximate 
h(x). Ignore if IHFC-0. In addition, XD(I)<.'CD(I+1) 
and it is reconmended chat che axial pcslcion be measured 
from che lower melt-solid interface, i.e., ;03(1)-0.0. 

rad 

VD(I) 

The surface temperature represented by che Ich knot used 
co approximate h(x). Ignore if IHFC-0. 

above 

■nelcing 
temp. 

READ ( 5 . 1 0 ) P , -KO , .YN ..MTERM , .MSUM , NGRID . .VR 
10 FORMAXOFIO.5,4110) 

P 1 

aO 

Melt Zone Peclec Number 

Axial position of lower raeic-solid interface; XO-0.0 is 
recommended. 

dimension- 

less 

rad 

CM 

Axial position of upper melt-solid interface 

rad 

••(TERM 

Sec CO Zero 


•(SUM 

The first .NSUM terms of che expansion in Equations 
(2.2.1i) are used co approximate Che melt zone temperacurt 
distribution, T(x,r). MSUM muse be less chan 21. 


MCRID 

(XIJ-XO) /NGRID is che grid size employed in System 
(2.2.20). In addition, che melt zone temperature dis- 
tribution is output for NGRID/ 10+1 axial values from 
XO CO XN . NGRID may not exceed 500 and must be 
divisible by 10. 


.JR 

The melt zone temperature distribution is output for 
NR+l radial values (rad) from 0 co L. NR -la- not 
100. 



t All lengths will be measured in radius 
rod is wide. The temperature scale will 
point. 


(rad) units, e.g., 2 rad is as long as the 
be °K above or below the material melting 
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TABLE A-2 PROBLEMS Pl-1 AND Pl-2 SOFTWARE INPUT (CONT) 


PROGRAM 

SYMBOL 



variable 

DETTNITION 

UNITS 

READ ( 5 . 1 0 ) P . XO , XN , MTERM , MSUM , :JGRID , NR 
10 FORMAT(3F10.5,4II0) 

P 

Lower solid region Feclet number 

dimension- 

less 

XO 

The semi-infinite lower solid region Is, for compucstional 
purposes, truncated to a finite length. XO is the axial 
position of lower end of this truncated region. Review 
Che end of Section ^.3 for details. 

rad 

XN 

Axial position of lower melt-solid interface. 

rad 

MTERM 

Integer parameter determining the size of system used 
CO compute the lower solid region's surface control 
function. See Equation (3.0.24). 


MSUM 

The first MSUM terms of the expansion in Eqiutions 
(2.2.14) are used to approximate the lower solid tempera- 
ture distribution, T(x,r). MSUM must be less chan 21. 


NGRID 

(XN-XO) /NGRID is the grid size employed in System 
(2.2.20). In addition, Che lower solid temperature 
distribution is output for NGRID/ lO^-L axial values from 
XO CO XN. NGRID may not exceed SOO and must be divisible . 
by 10. 

■ 

NR 

The lower solid temperature distribution is output for 
NR+1 radial values (rad) from 0 to 1 . NR may not exceed 
100. 



READ ( 5 . 90 ) Rjes . Ria , RL , NSYS 
10 FORMAT(3E20. 10,110) 

RKS 

Conductivity of material in lower solid region 

cal 


rad sec 

RKL 

Conductivity of material in melt zone 

cal 


‘’K rod sec 

RL 

of Equation PZ4, Figure 1-2. RL is the product 

cai 


of the growth race, the solid materials density and the 
latent heat. 

sec 

NSYS 

Number of terms in expansion of lower solid region's 
surface control function (see Equation (3.0.23)). NSYS 
may not exceed 20. 
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ORIGINAL PAGE 19 
OF POOR QUALITY 


TABLE A-2 PROBLEMS Pl-1 AND Pl-2 SOFTWARE INPUT (CONT) 


PROGRAM 

SYMBOL 

VARIABLE 

DEFINITION 

UNITS 

READ (5,21 )MINTE?M .MAXTERM.DELTERM, MINNSYSM.MAXNSTS .DELUSTS 
21 FORMAT(8I10) 


The sofcvare is designed Co compute Che lower solid 
region's surface concrol function for many combinations 
of MTERM and NSYS. To do this, the user must supply che 
bounds and increments for che cases desired. 


MINTERM 

Lower bound on MTERM. 



MAXTERM 

Upper bound on MTERM. 



DELTERM 

Inceger incremenC for MTERM 


MINMSYS 

Lower bound on NSYS 



MAXNSYS 

Upper bound on NSYS 



j DELNSYS 

Inceger incremenc for NSYS 


READ (5,50) lOPTION , aiP 
50 FORMAT (110, FIO. 5) 


Afcer che lower solid region's surface concrol function, 
f(x), is decermiaed, che user may specify certain modifica- 
tions of che surface concrol function. These options are 
principally of use when che surface concrol function is 
above che material's melting point on portions of che lower 
solid region's surface. If f(x) is interpreted as che '^K 
above or below che melting point, then let A be the lower 
endpoint of che largest subincerval of che form [a,.'CN] on 
which f(x; is not positive. If A < XO, reassign A to be 
XO. Let (.XMIN.FMIN) be che minimum poinc of f(x) on che 
Interval [A,XN J . 


lOPTION 

• 0 if f(x) is not to be 
be assigned any value. 

modified. In this case, CLIP may 
e.g. , zero. 



« 1 if f(x) is to be redefined as; 




f(x) if X > .XMIN 



C 1.3%; 

FMIN otherwise 



■ 2 if f(x) is to be redefined as: 



f (‘f ^ 

f (x) if X > .XMIN 




min |f(x), clip} otherwise 
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ORIGINAL PAGE (9 
OF POOR QUALITY 

TABLE A- 2 PROBLEMS Pl-1 AND PI -2 SOFTWARE INPUT (CONT) 


PROGBAM 

SYMBOL 


VARIABLE 

DEFINITION 


These options are graphically illustrated below 

0 , 



K above 
Che melting 
ceap 


-L 0 


4. FMIN 


UNITS 


10 


READ ( 5 , 1 0 ) P , XO , XN , ITEPM .MTERM ,MSUM , NGRED , .NR 
FORMAT(3F10.5,4I10) 


P 

xo 

XN 

MTERM 

MSUM 

MCRID 

.NR 


Peclet number for upper solid region 
Axial position of upper interface 

The semi-infinite upper solid region is, for computational 
purposes, truncated to a finite length. XN is the axial 
position of the upper end of this truncated region. 

Number of equations (n-1,2, • • • ,MrERM) of the type given in- 
Equation (3.0.28) used in least squares generation of upper 
solid region's surface control function. 

The first MSUM terns of the expansion in Equations (2.2.14) 
are used to approximate the upper solid regions' temperature 
distribution, T(x,r). MSUM must be less than 21. 

(XN-XO) /NGRID Is the grid sire employed in System (2.2.20). 

In addition, the upper solid region’s temperature distribution 
is output for NGRID/10+1 axial values from XO to .XN. NGRID may 
not exceed 300 and must be divisible by 10. 

The upper solid temp. distribution is output for NR+l 
radial values (rad) from 0 to 1. NR may not exceed 100. 


Dimension- 

less 

rad 

rad 



READ ( 5 , 90 ) RKS , RKL , RL , .NSY3 
90 FORMAT (3E20. 10, no) 

1 

RKS 

Conductivity of material in upper solid region 

cal 



°K rad sec 

RXL 

Conductiviev of material in melt zone 

cal 



°K rad sec 

U 

4 of Equation FZ2, Figure 1-2. RL is the product of the 

cal 


growth rate, the solid material's density and the latent 

sec rad^ 


heat 
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ORIGINAL PAGE F3 
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TABLE A-2 PROBLEMS Pl-1 AND Pl-2 SOFTWARE INPUT (CONT) 


PR0C8AM 

SYMBOL 

VARIABLE 

DEFINITION 

UNITS 

NSYS 

Numbar of earns in axponancial axpansion of cha uppar solid 
region's surface control function (sea Equation (3.0. 3D) 


READ (6, 23)MINTESM.MAXTEi!M,DELTESM,MlNNSYS .MAXNSYS .DELNSYS 
21 YORMAT(ailp) 

MINTERM 

MAXTERM 

DELTERM 

MINNSYS 

MAXNSYS 

DELNSYS 


> As previously defined above but applied to 

the upper solid region. 


READ(5,50)IOPTION,aiP 
50 FORMAT(I10,F10.5) 

lOPTION 

Af 

8< 

ei 

pr 

ab 

so 

ar 

11 

0 - 

CLIP - 
GM3-J 

ter the upper solid region's surface control function, 
t) , is datamined the user may specify certain modlflca-. 
ons of this surface control function. The options are 
incipally of use whan tha surface control function is 
ova tha aaterlal's malting point on portions of the upper 
lid region's surface. The definitions of lOPTION and CLIP 
B siailar to their previous definitions above and are 
lustrated below. 

above 

the melting — 

temperature 

xo aaN lOPT-M.''— ^ 


y^IOPTlON-2 

^^^^^idPTION-l 






1 


r 




V 

a 


ORIGINAL PAGE ia 
OF POOR QUALITY 

An input sample is illustrated in Figure A-3. 

The output is clearly labeled for ease of use. The melc .one input 
data is first displayed followed by the melt zone temperature dist.j.xbuclon (given 
in table format) and interface gradients (see Figure A-4) . The lower solid 
region's material and system parameters and the values of MINTERM, ’ * * ,DELNSYS 
are next displayed followed by the required lower solid region interface 
gradient, B(r) of Equation (1.2.1) (again see Figure A-4) . 

For each of various acceptable combinations of MTERM and NSYS (recall 
the definitions of MINTERM, • • • ,OELNSYS) , a lower solid region surface control 
function is computed. For each of these cases, the values of MTERM and NSYS 
are first displayed followed by the coefficients (see Equation (3.0.23)) used 
to determine the surface control function. Using the surface csntrol function, 
the temperature distribution in the lower solid region is next displayed (in 
cable fonn'f') . The lower solid region's interface gradient is then given 
followed last by the relative difference (in the norm) between the required 
lower solid region interface gradient and the interface gradient obtained bv 
use of the surface control function (see Figure A-5) . 

After all the lower solid region cases (various combinations of MTERM • 

and NSYS) are given, the results for the upper solid region are displayed 
in a similar fashion (see Figure A-6) . 

To finish, the user must supply the two algorithms described at the 
end of Appendix A. 2. In addition, an algorithm to solve (in a least squares 
sense) an overposed system of linear equations mu.'.t also be provided for use 
in the subroutines S0LID2 and S0LID3 (see Appendix C.3). 

I 


t t 




' I 



The surface control function values 
column. 


can be read from this cable in the R-1 



i. 


• , 

I 


! 
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2.7 
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1,04027 




It, 442 

0.0 



^00 1 

0,rlO42l 

0.0 4.842 

10 

2f* 

0,40734 

•10,0 0,0 

10 . 

20 

500 1 

0,4180075 

0,038908 

0.238«i 
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QSER CONSIDERATIONS FOR PROBLEM Pl-3 SOFTWARE 


A. 4 


The data input procedure and output Interpretation for the software 
developed for Problem Pl-3 are the subjects of this section. Recall that the 
problem is to find a melt zone surface control function (h(x) in Problem Pl-3) 
which, for the sake of flat interfaces, is compatible with the a priori given 
surface temperature distributions of the two solid regions. At present, all 
the required data are input in the form of punched cards. The definitions and 
formats of the input data are summarized in Table A-3. 
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TABLE A- 3 PROBLEM PI- 3 SOFTWARE INPUT 


PROGRAM 

TOIBOL 

VARIABLE 

DEFINITION 

UNITS 


R£AD(S, 12)F,MSnM,NGRID,KR 
12 FORMAT (E20. 10.6110) 

P 

Peclec a\mb«r of upper solid region 

Dimension- 

less 

MSUM 

The first MSUM terms of the expansion in Equations (2.2.14) 
tre used to approximate the upper solid region's temperature 
distribution, T(x,r). MSUM must be less than 21. 


NGRID 

For computational purposes, the semi-infinite upper solid 
region is truncated to a finite length (SLENGTH) . SLENCTH/NGRID 
is the grid size employed in System (2.2.20) which Is used to 
approximate the temperature distribution in this crunceted 
upper solid region. In addition, the upper solid region's 
temperature distribution is displayed for NCRH)/1041 uniformly 
spaced axial valuas. NGRIB may not exceed 500 and must be 
divisible by 10. 


MR 

The upper solid region's temperature distribution Is output 
for NR-Hl radial valuas (rad) from 0 to 1. NR may not exceed 
100. 




READ(5,22)RKS,Ria,RL,SLENGTH 1 

22 FORMAT (4E20. 10) | 

RKS 

Conductivity of oiaterial in upper solid region 

cal 


°K rad sec 

RKL 

Conductivity of material in melt zone 

cal 


rad sec 

RL 

J. jf Equation F22, Figure 1-2. RL is Che product of the 

cal 


growth cate, the solid material's density and the latent 
heat. 

sec rad*^ 

SLENGTH 

For compucacionel purposes, the semi-infinite upper solid 
reg.i../n is truncated to a finite length, SLENGTH. For 
details, review che end of Section 2.3. 



rad 

READ(5,22)Q 
22 FORMAT (4E20. 10) 

0 

Q is che length of che melt zone 

rad 

R£AD(5,16)IHFC,M 
16 FORMATdOIS) 


All lengths will be measured in radius (rad) units, e.g., 2 rad is as long as the 
rod is wide. The temperature scale will be °K above or below che material' melting 
point. 
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TABLE A-3 PROBLEM Pl-3 SOFTWARE INPUT (GOUT) 


Is. 


PROGRAM 

STMBOL 

7ARIA2LE 

DETTNITIQN 

UNITS 

IHFC 

> 1 if a cubic spline will be used to approxioace the surface 
cainperacura distrlbuclon, g(x), of the upper solid region. 



- 0 if the user will supply a functional form of g(x). In this 
case, the user anise insert this functional form of g(x) 
in the subroutine HFC (see the software list in Appendix 
C.4). 


M 

Number of knots used to approximate g(x) by a cubic spline • 
(IHFC-1). If IHFG-0, sec M-0. 


DO 32 1 - 1,M 
READ(5,22)XD(I),VD(I) 
32 CONTINUE 
22 FORMAT (4E20.10) 

XD(I) _ 

The axial position of the 1^^ knot used to approximate g(x). 
Ignore if IHFC-0. In addition, XD(l)<XD(I+lX XD(1)-Q and 
.'CD(M)-Q+SLENGTH 

rad 

TO(I) 

^ — 

The surface temperature represented by the I^^ knot used 
to approximate g(x). Ignore if IHFC-0. 



°K below 

melting 

temp 


12 


READ ( 5 , 1 2 ) P ,MSUM , NGRID , NR 
FORMAT (E20. 10,6110) 


MSUM 

NGRID 

NR 


Same deflnicions as above buc 
applied CO the lower solid region 


READ ( 5 , 2 2 ) RKS , RRL , RL , SLENGTH 
FORMAT (4E20. 10} 


RKS 


RKL 


RL 


SLENGTH 


Conduccivlcy oc macerlal In lower solid region 
Conductivity of material in melt zone 


of E^iuacion FZ4, Figure 1-2. RL is the product of 
Che growth race, solid material's density, and latent heat 

For conpucaclonal purposes, the semi-infinite lower solid 
region is truncated to a finite length, SLENGTH. For details, 
review the end of Section 2.3. 


cal 


rad sec 
cal 

rad sec{ 
cal 


sec rac 
rad 
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TABLE A-3 


PROBLEM Pl-3 SOFTWARE INPUT (CONT) 


PROGRAM 

SYMBOL 

variable 

DEHNITION 

UNITS 


16 

R£AD(5.16)IBFC,M 
FORMAT (1015) 


IHFC 

Sana daflalciona at praviously given but applied to eha 
surface cafflparacura distribution, £(x), of cha lover solid 
region. 


DO 32 I - 1,M 
READ(5,22)XD(I),YD(I) 
32 CONTINUE 
22 FORMAT (4E20. 10) 

XD(I) 

The axial position of the knot used to approximate f(x). 

Ignore if IHFC-O. In addition, XC(I)<XD(I+1) , XD(1) — SLENGTH, 
and XS(M)-0.0. 

rad 

YD(I) 

The surface ‘teoperacure represented by the knot used to 

approximate f(x). Ignore if IHFC«0. 

°K belov 

melting 

temp. 


12 

READ (5 , 1 2 ) P ,MSUM, NGRID , .'IR 
FORMAT(E20.10,6I10) 


P 

MSUM 

NGRID 


Sama definitions as previously given but applied 
CO the mal'c zone of length Q 


READ (5,16) MAXTERM ,MINTERM .MAXNSYS .MINNSTS , OELTERM, DELNSYS I 

16 FORMATC10I5) { 

1 1 


The sofevare is designed to compute the malt zone surface 
control function for various combinations of MTERM and NSYS. 
To do this, the user must apply the bounds and Increments for 
the cases desired. 


MAXTERM 

Upper bound on MTERM 


MINTERM 

Lower bound on MTERM 


MAXNSYS 

Upper bound on MSYS 


MINN SYS 

Lower bound~^n NSYS 


DELNSYS 

Integer increment for NSYS 
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An input sample is illustrated in Figure A-7. 

The output is clearly labeled for ease of use. The input data for :he 
upper solid region is output first followed by displays of the upper solid 
region's temperature distribution (given in table format) and interface gradient 
(s«e Figure A-8). The lower solid region follows in a similar fashion 
(Figure A-9). Using Equations FZ2 and FZ4 of Figure 1-2, the required 
melt zone interface gradients are computed and then displayed along with 
the melt zone Input data (see Figure A-10). For each of various combinations 
of MTERM and NSYS (recall the definition of MINTERM, •••, DELNSYS), the 
expansion coefficients of the melt zone surface control function (see Equation 
(4.0.18)) are output. Using the computed surface control function, the melt 
zone temperature distribution (given in table* form) and interface gradients 
are displayed next followed last by the relative difference (in the norm) 
between the required melt zone interface gradients and those obtained by use 
of the surface control function (see Figure A-11). 

To finish, the user must supply the two algorithms described at the 
end of Appendix A. 2. Also, an algorithm to solve (in a least squares sense) 
an overposed system of linear equations (for example, see [l7, Chapter 5] ) must 
be provided for use in the subroutine MELTl (see Appendix C.4 for code listing). 
In addition, a numerical integration routine is required for use in the 
subroutines INTEGLl and INTEGL2 (during software verification, Che numerical 
quadrature code CADRE was used and is available in the IMSL package or from 
the open literature [l6. Chapter 7j). 


' The melt zone surface control function can be read from this cable in the 
R*1 column. 
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Figure A-7 Sample Inpui: For Problem Pl-3 S<. .vare 
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Figure A-8 Sample Upper F Id Regioa Output For Problem Pl-3 Software 
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Figure A->< Semple Upper Solid Region Output For Problem Pl-3 Software (Cont) 
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Figure A-9 Sample Lower Solid Region Output For Vroblem Pl-3 ' Software 
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Figure A-9 Sample Lower SoliJ Region Ouiput For Problem Pl-3 Software (Cont) 





Figure A-10 Sample Output Of Required Melt Zone Interface Gradients For Problem Pl-3 Software 
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Figure A-11 Sample Output Of Melt /.onr Sni'fnce Control Functio.n 
Temperature Diuti itiui i.rn, '.n ' Interface Gradients 
For Problem Pl-J Software 
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Figure A-11 Sample Output Of Melt Zone Surface Control Function 
Temperature Distribution, And Interface Gradients 
For Problem Pl-3 Software (Cont) 
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APPENDIX B. CONVERGENCE OF EQUATION (2.3.11) 


Recall, from Section 2.3, Equation (2.3.11) 

“[*Q 4" j • 


where 


a X 

a 


[h-±] V “H* 


V 




a » (P + S )/2 
n n 


(2.3.11) 


r 


H 

i 

A 


and 


a =• (P - S )/2 

n n 

, C- -Sc 

a - \ G e “ < 

n S J n 


A « - B + - 

Q n n 2 


at 




In this appendix, under the assumptions or Section 2.3, it vill be shown that 

lim 0 (x) = 0 
n 


For convenience, suspend the use of the "n" subscript and define 

Since G approaches zero as x proceeds to negative infinity, if given some 
c>0, then there exists some N<J such that ] ^(t) 1 < e if t< N. Then since 
a>0. 


lim 


ax 


G(t;e"“^dt 


< lim e 
.^•^^00 


ox 


f" 

c 0 


1 e-^^dt 

•'x 

■'n 


e/a 


B-l 
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Hence, since cha above e is arbl'-rary and a>0, che first sunanand of (2.3.11) 
converges to zero as x proceeds to negative infinity. Next, for the second 
summand of (2.3.11), 



because 11m G * 0. Hence lia 9(x) » 0. 

-iM 


A clever man understands che need for 
proof . 


— Proverbs 14:15 
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APPENDIX C.O COMPUTER CODE LISTS 


C.l INTRODUCTION 

The computer codes developed to solve Problems Pl-1, Pl~2, Pl-3, and 
?2-l are given in this appendix. The codes themselves contain numerous comments 
correlating portions of the codes with sections and equations in this report. 

The code for Problem P2-1 is listed in Appendix C.2 and is followed by the code 
for Problems Pl-1 and Pl-2 in Appendix C.3. To finish, the code for Problem 
Pl-3 IS listed in Appendix C.4. 

C.2 COMPUTER CODE LIST FOR PROBLEM P2-1 

The computer code for Problem P2-1 is listed in Figure C-1. Before 
using this code, the user should review the remarks given at the end of Appendix 

A. 2. 
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ccrccrccccccccrccrcccttcccccccccccccccrccccccccccccrccrcccrcccccccccccccc 

cc 



CC 

cc 

purpose- 

CC 

cc 

COMPUTE 

steady-state TEHPERaTurf distribution And ThERCAL 

CC 

cc 

gradients for finite length translating fTLINOEP, This PROOLtM 

CC 

cc 

IS OESCRTBEn IN detail In SECTION 2.2 OP PINAL bCPORT (TO f.’ASA) 

cc 

cc 

- The control of float zone InTfRFaCES By THE U 5 C Of SFLEC'EO 

cc 

cc 

BOUNDARY CONDITIONS - BY SCICncf APPLICATIONS, tI<C. 

cc 

cc 

SOURCE- 


cc 

cc 

SCIFNCE 

APPI.ICATIOnS, INC. 

cc 

cc 

HUNTSVILLE, ALABAMA 

cc 

cc 

authors^ 


cc 

cc 

LARRY M 

. foster 

cc 

cc 

JOHN MCINTOSH 

cc 

cc 

REFERENCE 

• 

cc 

cc 

• The control of float zone Intprfaces by the usC of scLcr-EO 

cc 

cc 

boundary conditions - 

cc 

cc 

(FINAL 

REPORT - SaI-63/5034 ♦ HtO 

cc 

cc 

SCIENCE 

APPi ICATiONS, INC 

cc 

cc 

remarks- 


cc 

cc 

- software nevELOPEO and TESTED HN COC 'bOO/NROft ANn 

cc 

cc 

UNIVaC 

noB 

cc 

cc 

- All EOIIATtQNS referenced IU cnoc RELOH are contained VI HE 

cc 

cc 

final report- 

cc 

cc 


- The control of float zone interfaces rY The L’SE r.f 

cc 

cc 


SELFCTED boundary COl.niTIONS 

'cc 

cc 

INPUT Variables and functions- 

cc 

cc 

p 

- PECLET number 

cc 

cc 

MSUM 

- Number of teri;s tn series expansion of 

cc 

cc 


TEHPERATURE distribution (ImE DESIrEO SOLL'TlOf!) 

cc 

cc 

xn 

. axial position OP LOWER EnO OF CyLINoER 

cc 

cc 

XN 

• axial ROSITIOi; op upper end of CyLINoER 

cc 

cc 

NGRID 

. number of DIViSinNS OF C'^LlNOF* aXIS U3C0 Ifl 

cc 

cc 


solution of 0. 0. p. BOUnoaR"' VALUP RROSLCm 

cc 

cc 


resulting proh Transformation of the poe modeling 

cc 

cc 


TMF TEMPCRaTUDE 

cc 

cc 

NR 

. number divisions of cylinder RAotUS used in 

cc 

cc 


OUTPUT OF TEMPfC^TiiPE DISTRIBUTION 

cc 

cc 

IHFC 

- 1 IF A discrete mA^a point form nr The surface 

cc 

cc 


temperature is user provided 

cc 

cc 


- 0 IF A USER OCFIHEO PunoTIONAL Fnf?H of The 

cc 

cc 


surface TEMPCRaTURp is provided 

cc 

cc 

(XO,YO) 

- USER provided data PTS FOR THE AxiAL DISTANCE 

cc 

cc 


(XD) and corresponding surface temperature (TO 

cc 

cc 

M 

- number of data PTS, INPUT IF IMfr » 1 

cc 

cc 


SET TO 0 IF IhFC a 0 

cc 

cc 

MFC 

. USER RROVIOCD flp IHFC a 0) SURFaCE TEmCERaTuRE 

cc 

cc 


distribution 

cc 

cc 

AFC 

• USER provided temperature OISTRiiUTION ON T|)e 

cc 

cc 


lower (ayial dogitton a xo) fno Or the cylimocr 

cc 

cc 

BFC 

- USER provided TEmFERATuRE distribution on the 

cc 

cc 


UPPER (axial POGITtON a XN) END OF TmE CTLINOE" 

cc 

cc 

output VaRIABI.es- 

cc 

cc 

tmold 

« TEMPERATURE OISToISUTION ARRAY Jn THF CYLINDER 

cc 

cc 


(FROM xo TO XN AXIaLLY, hITh NR DIVISIONS OF ThE 

cc 

' cc 


RADIUS) 

cc 

cc 

GRADXO 

« axial thermal gradient array aT »0 

cc 


Figure C-1. Computer Code List for Prob]em P2-1 
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CC 5»*f)XN - *XI*U THERMAL GRADIENT *RR*Y *T jfN CC 
CC USER SUPPLIED m*Tmh 4TIC*L SOr^iARP* CC 
CC - A LEAST SOUtRES ALGORITHU Td aHLVE nvER POSED SYSTEMS QP CC 
CC linear EOUaTIPNS (required Hi SiKJRQUTtNE COCFS.i CC 
CC - AN algorithm to evaluate SCSSFL functions (REoUIRfO IN CC 
CC subroutine ,10) CC 

CC CC 


ccrcccccccccccrccrcccrcccccccccccccccrccccccccccccrccrcccrcccccccccccrcc 

real JIa.iiLam 

CPMHON/CI /R| iM0(20).J! (20),.flLA“(20) 
CnMMON/C 10 /ASCRTP( 20 ), 8 SCRIPC?nA 
COMMON/RpAOt /P»MSUM,xO,XN,(,'GOlO,yR 
COMMON/Cs/Pf J0l).R8I(20, I01),SQ.»i(20) 

CnMMON/C20/A (500), 6(500) »C (500 ),P(5rtO).v (500 ),OfTA(S 05J, GAMMA (SOS) 
COMMCN/C21/THETAB(20,50S),TwCLDYiOj ), T(‘0),GRAOx(i( J01),GRiOx0(loi) 
character*!? R!S,ALPma(6) 

CHARACTER*IA STARS, STAR(6) 

DATA RIS/'R. '/,StARS/'«****************'/ 

DO 210 Lai, A 


AI.PHA(L)sRIS 
STARCL)3STAR3 
PlO CONTINUE 

ccrccccrcccccccccrcccrcccccccccccccccrccccccccccccrccrcccccccccccccccccc 

rr CC 

CC 

RLAMD(M)aROOT OP Jo BESSEL PCN CC 

J 1 (M)aJi (RL aMJ(M)) mmERE Jl Is aCSSEL PCN CC 

JtL*M(M)aJl (H)/RL*M0(M) CC 

CC 


CC 

CC 

CC 

CC 

CC 

CC 


kC CC 

ccrccccrccccccrcLfCccrcccccccccccccccrccccccccccccrccrcccccccccccccccccc 

RLAMO( P-jc' ao««255S7T 
RIAMCC 2;*S, 5200781 103 


,5200781 
RLAMC( 3)a8,653T2TR129 
Rl.AMO( inali ,79I53««3RI 
RL*m 0( 5)ala, 9309177086 
Rl.AMO( 6)al8, 0710639679 
RI.AMO( 7)»2l .2116366299 
RLAM0( 8>a2a,35?oTl53"S 
R|.amO( 9ia2?,a93<iT91T20 
RLAMO( 101 S3n. 6306060684 
RI.AMO(ip« 33, 7756202136 
RLANO(12)»36,9iT0983537 
Rl. AMO (I 3) *40.0554257646 
RL AMO () 4) »45. 19979171 32 
Rt AM0(i5)a46. 301 1883717 
RLAMO( 161*49.4826098974 
RLAMO(l7)»S?.62405l8al I 
RLAM0(18)*55. 7655107550 
RL*mO( 19) *55. 9069839261 
RLAMO(20>*62. 0484691902 
Jl( 1)*0. 5161470973 
Jl( 2)*«0. 3402648065 
J»{ 3)*0,27t4522999 
Jl ( 4)*«0.252459831« 

Jl( 5)10.2065464331 
Jl { 6)»-0. 1577258030 
Jl{ 7)*0. 1732655902 

Jl ( 8)*-O.U17015507 

Jl{ 9)»0, 1521812138 
Jl (tO)*-fl. 1041659777 
Jl (1 l)*0, 1372969434 
<)) (12)*-0. 1313246267 
Jl (13)*0, 1260694971 
J) (|4]*<*fl.l?l39562o0 
JI(15)*0, *172111989 
Jl (l6)*-n.ll3429l9?6 


URIGINAL FAui: 1*:^ 
OF POOR QUALir/ 


Figure C-1. Computer Code List for Problem P2-1 (Cont) 


C-3 




Jl ( !7 )bO. 10«49| UllO 
Jl Cl8)a-o.ln48«TSAA3 
Jt (19)«0,IOV»59S729 
Jl {20)«-0,lrtl?93a989 

DO to Z>t.2fl 

JtLAH(I)«.n tp/«L*HO(I) 
SQJtcnaJi (I)*Jl'I) 
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10 CONTIHUC 
C*l.L INPUT 

ccrccrcrccccccrccrcccrcccccccccccccrcrccccccccccccrccfcccccccccccccccccc 

cc Cf 

CC find COCFS pop PESSCL EXFiHSlOMS Op /1(P)«*(1) Atio BfP)-0Ct) CC 

CC SPE equations (2.2. J7) AND (2.2'.18) OP FfNAL «EpORT CC 

... . CC 

ccrccccrccccccrccrcccrccccccccccccccrccccccccccccccccrcccccccccccccccccc 

CAu AFcn.A.AOFn 
CALL BPCri.n.pOFi) 

on 20 isi.ifti 

«(I)a(T-13 *0.01 
RHOLOaP(I) 

CALL APC(RHOLn.ANS) 

AdlaANS-AQPl 
CALL BPC(RH0L'1.ANS] 

8(Z)aANS-nGPl 
20 C'lNTINUe 

CALL COEFS(P«a. 101*20. ASCntF) 

CALL CO£FS(R, 8,101. 20, 8SCRIP) 

ccpccccrccccccrccrccccccccccccccccccrrcccccrccccccrccrcccrcccccccccccccc 

cc 

cc SOLVE FOP theta 8AB OF EOUATIONP (2.2,19) gy Soi'viNf, THE CC 

cc 7RinlA60N*L STSTeh (2.2.20) - SPE FIN*l REPORT CC 

CC CC 

ccrccccrccccccrccrcccrcccccccccccccccrccccccccccccrccrcccrcccccccccccccc 

OXa(XN.xn)/NCPID 

OX2aOXtOX 

u'Nepic.i 

on 30 Ma«,.H8UH 
00 ao lal.L 

A(I)al .A*0X*P/2.0 
B(r)a«2',o«OX2«RLA«0(N)«PLAMr*(H) 

C(t)at ,n«OX*P/2.0 
Xaxn*I«OX 

CALL GBaRIP.X.ANS) 

0(i)«OX2*AnS 
ao continue 

0(naOn)-(1.8»OX*P/2,0)«Ar.cRTP(H)*SQJi (m). 0,5 

0(L)aO(L)-(l .0-OX*P/2.0)«flScRTr(H)»SOJ| (M)*0.a 
CALL TPlDACrL) 

on so laa.NCPiD 

IIaI-1 

TMCTAB(M.I)av(II) 

50 CONTINUE 

NSTOPaNCPlO*! 

THETa 8(H, |)aASCPTP(“)»Sajl (mi >2,o 

Thetas (M,NST0F)aBScPiP(M)*3oji ('*)/?.<) 

30 continue 

OPal ,0/NP 
NPSTOPaNR^l 
00 60 lal.NPSTOP 
R(I)a(T«naOP 
OO 65 Hal.NSUH 
Fsi(H,r)8p(M,Ptn) 
oS CONTINUE 

60 continue 

cc print TCMPEPATUPES rr 


Figure C-1. Computer Code List for Problem P2-1 (Cont) 


C-4 


ro roRH*T; iMi,a5X,fl9H 


rM 9 C»»iUBE 


nr 5 T»i 0 ijTlo 


!n) 

IFUASaO 

H»IRHT »6 

MlEFTal 

tan cnNT:wuE ^ 

rr(NRSTOP.LF.M«icHT)iFLAc«i original page IS 

H»ir.HT.M;NOfNRSTOP,-»»ICHT) OUALITY 

*<»ITE(6,l90URCJ),J«MLeFT,l!RIr,HT) OF POUR ViUMI-P 1 I 

<90 F 0 RM*T(/////, IM , 17 *, 6 (FJ 2 .«, 5 X>) 

IM 0 l. 0 »H»TCHT-HLEf T» I 

WRITE (6>2a7 U alpha (L)«L ai,IHCLD^ 

?6l FrmHAT(lH+,(7X,AA17) 

WRITE (fa >2661 (STAR(L) >i.aUlHOLO) 

^68 FnRHAT(lH 0 »l 5 X.faAl 7 ) 

00 200 Ial,NST0P 
IIaNSTnP»i-I 
XBXo«(il-n*oy 

on 202 JaMLEFT.MRlCHT 

ccrccficrccccc:t:ccrcccf;cccccccccccccccrccccccccccccr.ccrcccccccccccccccccc 

CC Cil 

CC OETeRHIHF TfmPERaTURE at (X.OfJO rc 

CC SEE equation (2,2,1**) OF FIHAI »£P0RT ' 

0 " . . . c'” 

ccrccccrccccccrccrtc-'cccccccccccccfcrccccccccccccrccrcccccccccccc'^v.ccr.w 

THQUOlJlat.O 

on 2 oa Mai.MsuH 

THni.O(J)aTHflLO{J)*2,0«rsifH,J)»THtTA8fH. li '(") 

20 " CDNTINUF 

CALw -'FC(X,ANS) 

CHOLn vJlaTHntO( J)*ANS 
20? CONTINUF 

wRITF( 6 . 2071 X, (TMOLO(J).JaHi £FT,HRIGwT) 

?07 FnRMAT(3H Xa, Fin, fa, 5 h - ,fa{Fl9.8,2Xn 

20" CONTINIIf 

IF(IFUA 6 ,FQ. 1 )C 0 to 220 
MRIC:pT,hRtG :T** 

HLEFT9HLEFT*fa 
50 TO ISO 

22" continue 

ccrccccrccccccrccrcccccccccccccccccccrccccccccccccrccfcccccccccccccccccc 

CC CC 

CC compute thermal GRaOICN.S at XaFO ANO XM CC 

CC .. , . CC 

ccrccccrccccccrccrcccccccccccccccccccrcccrccccccccrccrcccccccccrccctcccc 

wRI7E(fa,7l) 

71 form*T( iHi,aTx,i 5 H Therm*) ;r.oie-ts) 
WRITE(fa,72)»fl,XN 

72 Ff1RHAT(///,4aX, ihR, 5X, UMCn*0, aT Xa,HO’.'v)aH GR*0, AT Xa,Fm.5 

2 ,//) 

DO 230 lai, ioi 

R(i)p(i«n*",oi I 

on 200 walfUSUH * 

PSr(M,I).FCK.B(I, 5 

?o 0 continue 

230 CONTINUE 

00 a. ox 

on 250 lai, -31 
on 2faO Jai,5 
T(J)bO, 0 
00 270 Hal ,FSUM 

T(,I)aT(J)«2,0»PSI(H,X)"ThETAe(M,j)/SOJl (Hi 
27" CCinTTNUC 

XaX0*(J„n«nx 
CAIL MFr(X,*NS) 


Figure C-1. Computer Code List for IvOi^ien P2-1 (Cone) 
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ORIGINAL PAGE IS 
OF POOR QUALITY 


T(J)«T(.!U*NS 
2bfl CONTINlie 

on 280 JSA, 10 

T(J)a 0 .rt 
JHOLn»NSTPP-tO»J 
00 29n Mal.MJUM 

T(.n»T{J)*2.0«P3IfM,I),THFTA8(H,JHCLO)/Sc;tHH) 

29« COKTTMUF 

XaX0«(,IHfll 0 *-i)*OX 
CALL HrcCX,ANS) 

T(J)aTfJ)*AN3 
28a CONTINUC 

ccrccccrccccccrccrcccrcccccccccccccccrccccccccccccrccpcccccccccccccccccc 
cc , cc 

CC APPPOXIHATE ThEPMAL GPADIEtiTS AT CYLlNOfP E><0S CC 

CC • SEE caiUTrCNS (2.2.21) AliO (2‘.2.22) 0^ PINAL bEPOBT CC 

CC .. . CC 

Lcrccccrcf' crccrcccrcccrcccccccccccrccccccccccccrccfcccrccccccrccccrcc 
CP A" I)a(- 3 *T (6)t 16*T (7)-3ii*Tr8) ♦u8*T(9)-2fi*Tf J 0) ) / CI2*0C) 
CPArA..(I)sC-3*T(5)tl6*T{«)-3o.T>3)*P8*T(2)-2S»Tf 1 ))/( 12*0X5 
W»ITEC6«300)RCI5 .GR*OX0(I)fGn*DTl'!(I) 

AOO FORMATdH ,39X,Pe,A.3X.En.<>,7X,C17.<>) 

<SA CONTINUE 
STOP 

CNO ... 

ccrccccrccccccrccrcccrcccccccccccccrcrccccccccccccrccrcccccccccccccccccc 
cc . cc 

cc THIS SU8ROUTINE APPROXIMATES (Rv FJNItE DIFFF»EnCE5 G CAR OF CC 

cc EOUATION (2‘.2.1«|) 0^ final RCPOpT CC 

cc .... cc 

ccrccccfccccccrccrcccccccccccccccccccrccccccccccccr.ccrcccrcccccccccccccc 
subroutine GBAfifH.X.ANS) 

RFAL Jl,.HLAM 

COMMON/CI/RI AMO (20 5 . jl (20) , J1LA»(20) 

common/reaO' /P.M suM.xo,xN,fiGrtn,MR 

EPSLONaO.Ol 

XIaX-EPSI.ON 

X23X+EP81.0N 

CALL HFC(X,iN3) 

CALL HFC(Xl.ANSl) 

CALL hfC(X2,AnS2) 

G»P*{ANS2-ANSn/(2.R*ePSCQI.) 

CaC-(ANS2*ANSl«2.0*ANS)/(EP3LnN»ePSLCN) 

ANSaC'JtLAMrM) 

return 

ENO 

ccrcccccccccccnccrcccccccccccccccccccrccccccccccccrccrcccccccccccccccccc 
cc cc 

cc This subroutine provides for data input cc 

cc cc 

ccrcccccccccccccccccccccccccccccccccccccccccccccccrccrcccccccccccccccccc 
subroutine input 

C(5MM0N/REA0l /R.MSUH.XO.XN.NGRTO.IjR 

cnMHON/C26/yO(io(n,Tonoo),cifa,ioo),- 

COHH0N/C25/THFC 

OIHENSION C(R.IOO) 

Eoui Valence (Cl n» n»c(i. : 5) 

ccrccccrcccccccccrcccccccccccccccccccrccccccccccccrcci-cccccccccccccccccc 
cc cc 

cc SPLINE INPUT option cc 

cc . cc 

ccrccccrccccccrccrcccccccccccccccccccrccccccccccccrccrcccrcccccccccccccc 
hRITEIB.R) 

5 FORMAT! imi,?7X,20HI NPUT aaT*) 

READCB, 1a)IhFC,m 
IP dMFC.NE.I) GHTObO 


Figure C-1 . 


Computer Code Lise for Problem P2-1 (Cone) 




wRITE(«i,-<0)>* 

JO FnBMAT(//////,<JSH TMf SURFACC TrMPERATURF orSTPrGUTinN IS APP»OXIH 

IaTED 3Y the cubic spline ThPCurh The following, fa, (X, tE*<P) 0 
•aTA points,///, J7H X surface Temp.) 


OP S2 I»1,H 
PFAn(5,2^)xrim,Y0(n 
22 pnRMAT(aE20'.10) 

16 FPRMAT(2H0) 

NRIT£C6,3a)xDfI),Yp(l) 
J« FPRMATCIH ,?E20,10) 

J 2 CONTI'lue 


ORIGINAL PAGE (S 
OF POOR QUALITY 


CAU, COFr.EN 
hO CONTINUE 

RfADCS, injP.xo, xn,m3um,ngRI0,nr 
10 FORMATCJPlO.S.allO) 

i*BITEC6,20)P,X0»XN,MSUH,NCnir<,NO 

20 FPRMATt///,';6H p XO XN mJUM nGRIO 

1 nR,//,ih .3El2.a,I5.2I7) 

RFT'JRN 

£NO 

ccrccccrcrccccrccrcccccccccccccccccccfccccccccccccrccfcccccccccccccccccc 


cc cc 
CC This subroutine supplies LaTCRAI Su»F*cE TE-PEPiTURE cc 
cc - SEE EOI'ATtON (2. 2. a) OF FH-al EEPORt CC 

cc cc 


ccrcccccccccccrccrcccccccccccccccccrcrccccccccccccrccrcccccccccccccccccc 

SUBROUTINE hFC(X.AnS) 

CPMMON/C2S/THFC 


dimension CI7) 

IFdHFC.EO.I ) QOTO60 

ccrcccccccccccrccrcccccccccccccccccmrcccccccccccccccrcccccccccccccccccc 
CC cc 

cc USER SUPPLIFO lateral SURFaCC TfHPeRaTU»E CC 

cc cc 

ccrccccrccccccrccrcccccccccccccrcctccrccccccccccccrccrcccccccccccccccccc 

return 

ccrccrcrcccccccccrccccccccccccccccr.ccrccccccccccccrccrcccccccccccccccccc 
cc cc 

cc LaTFRaL surface temperature rPn«IDED by spline cc 

cc FIT OF USER supplied DATA C*D,Y'') CC 

cc CC 

ccrcccccccccccrccrcccrcccccccccccccccccccccccccccccccccccccccccccccccccc 

60 CALL 3PLJNEfX,AN3) 

return 

END 

ccrcccccccccccrccrcccccccccccccccccccrcccccrccccccrccfcccccccccccccccccc 
cc cc 

cc THIS subroutine suPRlIES RaOIal 'EhpeRaTuRE DIStRIBuTIEN cc 

cc ON lower ENn OF cylinder - SEE eOUaTION (2.2,2) OF FINAL REPORT CC 

cc . . cc 

ccrccccrcrcccccccrcrcrcrcccccccccccrcrccccccccccccrccrcccccccccccccccccc 
subroutine aFC(R,anS) 
common/reaoi /p,MSUH,xo,xN,tjGf(rn,NR 

ccrccccrcccccccccrcccrcccccccccccccccrccccccccccccrccrcccccccccccccccccc 
cc . cc 

cc USER SUPPLiFO lower END TEMPERATURE A(R). CC 

cc . . . cc 

ccrccccccrccccrccrcccrrccccccccccccccrccccccccccccrccrcccccccccccccccrcc 

CALL HFC(XO.ANS) 

rfturn 

ENO 

ccfccccrcccccccccrcccrcccccccccccccccrccccccccccccrccrcccrcccccccccccccc 
cc „ , . cc 

cc this subroutine supplies radial temperature nISTRIBuTirN cc 

cc ON UPPFR END of cylinder • SEf fCUaTIPN (2.2.3) OF FINAL REPORT cC 

cc cc 


Figure C-I. Computer Code List for Problem P2-1 (Cone) 
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ORIGINAL "i 

OF POOR QUALITY 


ccrccccrcrccccrccrcccrcccccccccccccrcrccccccccccccrccrcccccccccccccccccc 

SU0»aUT!NE BFC(I»»*NS) 

COMMOH/PE*Ot/P.HSUH,XO,XN,t;GRir>jKB 

ccrcccccccccccrccrcccrccccccccccccccccccccccccccccrccrc'ccccccccccccccccc 
cc cc 

CC USEP SUPPLIP* UPPER £*<0 TEMPEPatUPE 8fP). CC 

cc . cc 

ccrccccrcccccccccrcccrcccccccccccccrcrcccccrccccccrccrcccccccccccccccccc 

CALL HFC(XN,ANS) 

return 

END 

ccrcccccccccccrccrcccccccccccccccccccccccccccfccccrccccccccccccccccccccc 
cc cc 

CC this subroutine fits BESSEL SeRtES TO DATA BY LEAST SOCARES CC 

cc method - SEE EO'IaTIOnS C2.2.17} , (2.2.18) AND f2.2;23) CC 

cc OE FINAL REPORT CC 

cc .. . cc 

ccrccccrccccccrccrcccrcccccccccccccccrcccccccccccccccrcccrcccccccccccccc 

SUBROUTINE rOFFS(R,V,NR,NCCEr,CPEF) 
integer NR.NCOEF 

REAL F,RnOl),Y(10l),COEF(20),HK(«6 0) 

external f 


ccrccccrccccccrccrcccccccccccccccccccrccccccccccccrccrcccrccccccrrcccccc 


cc ''SEP supplied least squares liETuOD POllOmS here "0 nETERMiriE CC 
cc the COEFFICTENTS of equations (2.2. 17J anO (2.2;i8)'. THE CC 
cc SUBROUTINE TFLSQ BELOM IS THE IhSL LEaST SQUARES FUNCTION CC 
cc FIT routine CC 
cc . CC 


ccrcccccccccccrccccccrcccccccccccccccrcccccccccrccrccrcccccccccccccccccc 

CALL IFLSO(F,R,Y,NR,cnEF,NCOCF,MC* lER) 
IF(lEP.CO.j?q.OP.IER.E0.130)URTTE(6. 10 ) 

«0 EOflHAT(B#,H terminal error IN LEaST SQUaReS meThhC , SuBREUT I’.'E COEFS 
1 ) 

RETURN 

END 

ccrccccccccccccccrcccccccccccccccccccrccccccccccccrccrcccccccccccccccccc 


cc cc 

cc This FUNCTION EVALUATES THE ZERn ORDER BESSEL ElUJCTlON CC 

cc denoted in notation n 2.1 (III) . SEE FINaL PEPOrT C 

CC cc 


ccrccccrccccccrccrcccrcccccccccccccccrcccccccccccccccrcccccccccccccccccc 

PEAL function FfN.R) 

COHHCN/cI/Ri amO( 20)»J1 (20),JiLA“C20) 

X«RLAMD(N)«R _ 

CALL Jo(X.Y) 


F«Y 

return 

END 

ccrccccrcccccccccrcccccccccccccccccccrccccccccccccrccfcccccccccccccccccc 


cc cc 

cc this subroutine computes the jo BESSEL FUNCTION Y«Jo(X) CC 

cc . . cc 

ccrccrcccccccccccrcccccccccccccccccrccccccccccccccrccccccccccccccccccccc 

subroutine ,T0(X.Y) 

ccrccceccccccccccrcccccccccccccccccccrccccccccccccrccrcccccccccccccccccc 
cc cc 

cc user SURPLIFO jn FUNCTION PLACED HERE. In THIS fXamplF, THE CC 

cc IM3L BESSEL FUNCTION mmbsJO IS TLLUSTpaTeO CC 


ccpccccccccccccccrcccccccccccccccccccrccccccccccccrccrcccccccccccccccccc 

REAL MMBSJO 
YbMMBSJOCX.IER) 

return 

end 

ccrcccccccccccrccrcccccccccccccccccccrcccccrcccrccrccrcccpcccccccrcccccc 


Figure C-1. Computer Code List for Problem P2-1 (Cont) 
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4 «V 


cc 

cc 

cc 

cc 

cc 

cc 


SliBROIlTINf FOP SOL^IMG A SrSTFM OF UINf^P SIMULTANEOUS 

eouations Having a tridiaconal occfficiEnt matrtX, oIacohals 
APE STOPFO IN THE arrays 4, B, A.'jO C. T^E C^mPL'iCO 
SOLUTION VErTOR V(U , , . V(l) IS STOREo JN Inr ARRAY V, 


CC 

CC 

cc 

cc 

cc 

cc 


ccrcccccccccccrccrcccccccccccccccccccfccccccccccccrccrcccccccccccccccccc 

SOeROUTINC TRIOAG(L) 

COMM0M/C30/A(S0()),B(Sn0),C(500),C(Sfl0),V(500),eFTA(sOS)»GAfi«A(5ft5) 

BFTA(l)«fl(n 

GAHMAa)«0O )/BFT4(l) 

iFPua 

on t lalFPl.L 

8 ETA<n«afn-A(i)*c(i-n/ 8 CTAi’i-n 

GAMMA(n«fO(n«A(I)«CAHM(»(I»i')/ 9 eTA(I) 

1 CONTINUE 

vfL)sCAMMA(i ) ORfGfWAL F 3 

on 2 ksi,last POOR QUALiTY 


I=L-K 

VinsGAHMACn-CCDAVClAj) /beta ■( I ) 
2 continue 


return 

ENO 

ccrcccccccccrcrccrcccrcccccccccccccrcfcccccccccrccrccrcccrcccccccccccccc 


cc cc 
cc this subroutine ESTIhaTES lateral surface TEMPEoATUrE cc 
cc 8T USE OF cubic spline IF TMC UAER 5UPP<-tES A OtCCReTE SET nF CC 
cc LATFRAL surface temperatures, cc 
cc .. . . cc 


ccrccccccccccccccrccccccccccccccccccrrccccccccccccrccrcccccccccccccccccc 

subroutine SPLINECXInt, TINT) 

CnMMQN/CR6/X0(lOnJ,YoCI0O),Cl(a, J00),'A 
DIMENSION CTa, 100) 

EouiVALENCEfCi ti, n.ccunj 
IF(XINT-X0(I))2»1,2 
I YINTaYOin 
return 


2 Kb 1 

3 IF(XINT-XO(k*I) ) 6*«*5 
a YTNTaYDCKAl5 

return 


5 KaK*t 

IFC(m-k),GT' 0) r.QTOS 
IF( CM-kI.LE.O) KaH-l 

4 TINTsCXOlKTt )»XINT)*fCn AXjACxnTKTJ )«*lMT)A*?»0f3,K) ) 
TTNTaYlHTA(xINT«XO(M) )*(C f 2» K ) • /XlNT-xO t K ) ) ABafr ( 4 , K ) ) 

rfturn 


ENO 

tcrccrccccccccrccrrccrcccccccccccccccrccccccccccccrccrcccccccccccccccccc 


cc .. cc 
cc find The 3P, INE curve fit COCFFtCIEnTS. for USE IN conjunction cc 
cc “TTm subroutine spline. cc 
cc INPUTS • cc 
cc Ha NO, Of OATa pairs cc 
cc xn a array of X (ABCISSA) vaiufS cc 
cc TO a array of V (OROINaTeS) v ai UES CC 
cc outputs - CC 

cc C a 2-0 ARRAY OF SPLINE FIT COCFFICIENTS fu CnCFFlCICNTS CC 
CC PER triplet of OATA POInTr), cc 
cc cc 


ccrccccrccccccrccrcccccccccccccccccccrccccccccccccrccfcccccccccccccccccc 

SUBROUTINE rOFCEN 

COHMON/C?6/X0n«0), Y0< J00),CIf4,i00),M 
DIMENSION CM. 1(10) 

DIMENSION PMOO).ErlO(i),AUOO,3>,BnOn)»?n(in),nClOn) 


Figure C-1. Computer Code List for Problem P2-1 (Cent) 
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E0UlV*uENCEfCl(t#1 J.CCl, D) 


no 2 K»«,M 

n(K)»xofKtO-xno<) OF 

P(K)»0(K)/h. 

ECK)«(Yn(K*l)-VD(K))/0(K) 

no j K»a,M 

nc<>iie(K)-F(K-t ) 

*(li 2 )a>l.- 0 (n/n( 2 ) 

*(l»3)an(n/0(23 

4(2»2)«2.*fPf n*PC23 )-P(n*«(l, 2 ) 

n'«?)>B(2}/A(2.2) 

no a Ki 3 ,H 

*(K,2)*2.»(PrK-n+''('0)-PC'<-1)*A{K-l, J3 

P(»<)a8{K)-P(K-t 

*(K.3)aPtK1/A(K,2) 

P()03B(K)/t (Ki?) 

QaO(M-l)/0fM) 

A(ND,2)a.Q.A(Nn, 1 )«a(P<3) 

B(NO)-afl(M-i ),A(NO, n*a(H) 
Z(NO)aSrN03/A(NO»2) 

no P i«i.Nn «2 

KaNO'I 

Z (K ) a8 (K )«a(K. 3) *ZU»n 

Z(na-Aa,?)«Z(2)-A(l,3)*2(3) 

no 7 Kal,M 

0al,/{6,A0rK)) 

nCl*K)aZ(K1*0 

C(2»K)3Z(K*U»Q 

rC3*K)tyo(K)/o(K)-zco*Po;) 

r(tt»K)aVD(K + l)/0<K)-Z(K*U*n(K> 

MaMAl 

return 

FNn CoFCEN 
ENO 


ORIGINAL FAGi^ 

OF POOR QUALITY 


Figure C-l. Computer Code List for Problem P2-1 (Cont) 
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COMPUTER CODE LIST FOR PROBLEMS Pl-1 and Pl-2 


The computer code for Problems Pl-1 and Pl-2 is listed in Figure C- 
Before using this code, the user should review the remarks made at the end of 
Appendix A. 3. 


C-11 


ccn.i.v,u.Lui.LLu.u.r(.i.LLun.i.ccccccccuccrccccccccccccrccrcccccccccccccccccc 


cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

DO 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 


PURPOSE- fC 

COMPUTE the UPPPR **<0 LOHEn nrCTONS SUR'*CE CONtPOL ruccTirns cc 
SUCH TH*T n *T SOLIO-HELT XHTerfaCES *Rfc iCMlEVrO AS FCRMULaTEO CC 
IN PROBLEMS PI-I ANO Pl-a nr final report {TO MiSA) ThC COf.-TROL CC 


OP float zone interfaces by THE (isE OF SE>-ECTCn BnUNOARY 
CONDITIONS- BY SCIENCE APPLICaTtCNS, INC*. 

SOURCE- 

SCIFNCE aPPI.ICATICNS, INC, 

HUNTSVILLE, ALABAMA 
AUTHORS- 

LARRY M, foster 
JOHN MCINTOSH 
RfFFRfNCI- 

• The CONTROL OF float zone InTfEFaCES 0Y THF Use Of sclectfo 
boundary conditions - 
(FINAL RFPORT • SAI-«S/503a>Uu) 
science applications 
remarks- 

- software OEVELOPEO and TcSTfo PN COC T6O0/fi«0n 
uNivAc no6 

- ALL EQUaTtONS referenced H. code BELQM are contained in The 
final report- 

- the control of float 70 NE Interfaces rt the use of 

SELFCTEO boundary CONDITIONS 

INPUT VaRIaBLFS and functions- 

- peclet number 

- number of TCPMS tN series EXPANSjON oF 
TEMPERATURE OISTRIrUTiON (ThE^OFSIpED sOLL'TIOf!) 

- axial POSITION OF LO“ER EnO OF CyLInoER 

- AXIAL POSITION OF UPPER END OF CYLINDER 

- number OF GRID PnJNTS USEo IN NUmER1C*L 
solution of 0, 0, F. boundary valuf problem 
resulting from transformation of P.iE The modeling 

TEMPERATURE 

- number divisions CF CTLIMoEJ» RADIUS USED IN 

OUTPUT OF temperature DISTRIBUTION 
» SOLID thermal conductivity 

- LIOUIO THERMAL conductivity 

- product OF crystal growth Rate# solid densi'y, 
and latent heat of fusion 

- 1 IF A discrete OAT* ROINT form np ThE SURFACE 
temperature is user provided 
0 IF A USER defined Functional forh of the 
surface TEMPERaTuRf is PROMiDEO 

- USER PROVIDED DATA PTS FOR THE AirlAL DISTANCE 

(XD) and CORRCSPONoINC surface temperature (YO) 

- NUMBER OF data PYS, INPUT IF IHfp a 1 
SET TO 0 IF IMFC a 0 

. USER provided (IF IHFC ■ 0) SURpICE TEMFERATuRE 

function 

- This fpqcram cenfKates thf surface control 

functions for VARinUS COnbInaTIONS CF THE INDEX 


p 

H3UM 

xe 

XN 

NCR ID 


NR 

RK3 

RKL 

RL 

IHFC 


(XO,YD) 


MFC 

CASE limits 


CC 

CC 

CC 

CC 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

Cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

•cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 


Figure C-2. Compucer Code List for Problems Pl-1 and Pl-2 


ORIGINAL PAGE 
OF POOR QUALITY 


C-12 




cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 


m*xTERM 

oelterm 

hInnSYS 

M/tXN3Y$ 

OELTERM 


IPPTION 


tiMiTs mterh «no nsys fsce fo. • (a, a. 23 ) 

OP FINAL REPORT), Tf! define THESE et^MBTNATIONS 
hInTeRm « The mINihUh ali.outc v»Lur 
GF -term 

- the maximum *U,OUrO ViLur 
Of mterm 

- ImCRE«EnT oE mTERm FRON ri'JTCRM 
Tn m|XTER“ 

- The mINtmUm aLLOHfO VALUf OF nSYS 
' The MAXJmUm aLLOWfC ViLUf OF nSTS 

- I"CRE»<C‘'T OF NSYS FROM minnSVS 

To i-AXnSYS 

- 0 IF SOLID SURfafE control FUNCTJCN is TO 

remain unchanged 

- I IF SOLID SURfafC control FUNCTION IS TO BE 
CLIPPED (SEE OEFN, JN aPPENdIS A, 3) AT ITS 
minimum value (hINTHUM Val^^E EOUNq ;n SUBCLUT'nE 

LINSRCH) 

- 2 IF FUNCTION IS TO BE cLiPpEn (sEE nEFf.', If; 
APPENOIs *P3) at SnHE USER SRECTFlpD ViLUf (SfE 
CLIP) 

» USER SUPPLIED VAi UE OF mOdIFIFO jURFaCE CONTROL 
„ function (SEE appendix A. 3 FOR OEFC,) 

OUTPUT VaRTaBI ES- 

• TempeRaTURF OISToiBUTlON aRRay FnR EaCm REGION 

- SEE Subroutine mflt of c°of 

- melt zone lower interface GRaoIEwT 

- melt zone upper interface gradient 

- Thermal CRaoICnt at XO fOr each region 

- Thermal GRAoirNTs at xn fo» each region 
. array of COCFriCiENTS OF SOLIE RfGIOnS SURFACE 
CONTROL. (SEE EOuatIOnS (3.0,23) aCO (3.0.31) ) 

- The L2 relative oIFfeRenCe BETWEFf; The PequIreO 

SOLID REGIONS II.TeoFacE 5»AoIEnT 5 if.O THE 

interface graoientm Resulting from the usf of 
THE solid regions SURFACE CONTROL FUNCTIONS. 

USER SUPPLIED MaTmmaTICAL SOFTNaRf- 

• A LEAST SoUaRES algorithm TP fIT a fuNcTion Tn A lInCAR 

combination of selected functions (REOUlREO IN cUBBnUTINE 
COEFS . ) 

• an algorithm to evaluate flCSSFL functions (REoUIReO in 
subroutine ,10) 

• A numerical integration routine (REQUIRED IN subroutines 


CLIP 


ThOLO 

grade 

GRaD 3 

GRaDXO 

GRadxn 

COEF 

ERRL 2 


cc INTeGLI aNO INTEGL 2 ) 


cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 


cc 

cc 


cc 

ccrccccccccccccccpcccrcccccccccccccccrccccccccccccr.ccrcccccccccccccccccc 

INTeCEP DELTERm.OEI NSTS 
PE*L J1»JILAM»MMBSJ0 

COMHON/c1/R|.aMO( 20). Jl (20 ?, JIlam( 20) 
COMMON/C5/Pri01),PSIC20,l0l),SQ,Tl(20) 
cnMMON/C<»/CnEF(2o),RH 
COMMON/CiO/A3CRIP(20),aSCRIP(2(i> 

COHMON/C20/A (500) ,8(500 ),c (SOO) ,5(500) ,V(50 O) ,BfTa (S 05) ,CAMMA(5n5) 
COMHON/Cpl /THET ab (20,505) ,THOlPM 01 ) , T ( 1 0 ) , CRAOxt; ( 1 n I ) , GRACX 0 () 0 1 ) 
C0Mm0n/C22/TCASE,Tm£lTC3) 

CnMMON/'CP3/r.RA02(lOl),ORAD3(lon 

CnMM0N/C2a/RKS,RKL»PL/N3Y3 

CnHMON/C25/RRAOATO(lOt),GRAOATOnon,RMSi ( 20 ) , RhC ,. ( jo ) , 3 (20 ) , 0 , 

I AHAT) (20,tO),AMAT2(20,10) , AL2(4o, 1 0 ) , RhS ( 4fl ) , wwj; ( 1 5oO ) , 1 1 WK (20 ) 

COMMON/RFAOi/P,MTERM,m8UM,XO,XN,NGRID,nR 

CrMHON/FIXPT/IOPTIGN,XMlN,CMlN,rLlP 

CflMHON/C30/MlNTERH,MAXTERM,oCLTFRM,MiHNSvS,MAXNsYS,OELN8YS 

CrMHON/C3J/IHFC 

COhHON/C32/XO(lOO),TO(100),CIfa,iOO),M 
CPMH0N/C76/TFLAG2. IFLAC3 

ccrccccrccccccrccrcccccccccccccccccccrccccccccccccrccrcccpccci;cccccccccc 
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CC COMPUTE THCRMit OlfTRlBUTION InTEpP»CE BRAOiENTS FOR CC 

CC MELT ZONE CC 

CC CC 

ccrccccrccccccrccrcccrcccccccccccccccrccccccccccccrccrcccccccccccccccccc 

IFLAC2*1 

IFLAC3al 

ICASCil 

wR!TE{6f<J«31 

0B3 rnRMATnHt<*5«X, I8HM E L T Z n N £) 

CALL INPUT ^ 

CALL HFCCXQ.ANS) ORIGINAL PAGE IS 

TPELT (23 aANfl _ _ _ _ _ <M|A||*rV 

CALL HFC(XN,ANS) OF POOR QUALITt 

TMELT(3)»ANS 

call pelt 

ccrccccrccccccrccrcccrcccccccccccccccrccccccccccccrccrcccccccccccccccccc 

CC CC 

CC USING the lower solid RECICN supF*cE control FUmCTIpN, compute CC 

CC the lower solid region thermal distribution and interface CC 

CC GRADIENT, CC 

CC . CC 

ccrccrcrcccccccccrccccccrccccccccccccrccccccccccccrccrcccccccccccccccccc 

IHFCaO 

ICASEaE 

WRITE(6|30) 

30 FORMATCiHlfSEX.ZZHL OWED sOlIO 
CALL INPUT 

DO 2J mTER aMINTERM»MAXTER»‘..OELT'ERM 

MTERM*MAXTERMfMINTERM-MTER 

DO 22 N3T «MINN3YS,MaXNSY3,OCLNSYS 

NSY3sMAXN3YS«MINNSYS-n3Y 

NN»mTERMa2 

IF(NN,lT.NSvS)GO to 22 
CALL 30LID2 

IF(IOPTinN,FG.O) GOTOTO 
CALL LINSRCH(XMTN) 

CALL FUNC CXhIN.GmIN) 

IF(XMIN,r,T,fl,0.ANO,GMlN,LE,O,O) COTOao 
XMINaioonOO'.O 
aO XMIN»»N»XMIN 
70 continue 
CALL melt 

ccrcccccccccccp.ccrccccccccccccccccccccccccccccccccrccrcccccccccccccccccc 

CC CC 

CC determine Relative oifferencc BFTweeN requibeo i ower solid cc 

CC REGION interface GRADIENT anO ThE InTeRFaCE GRADIENT RESULTING cC 

CC FROM USE OF THE LOWER SOLID RECtON 3URFAc€ CONTbOL EUNCTIOfJ CC 

CC CC 

ccrccccrccccccrccrcccccccccccccccccccpccccccccccccrccfcccccccccccccccccc 

CALL ERROR 

WRITE(6.78RT MTERm.NSYS 

tSR FORmaT(//,ih ,Srt*.l2MF0R MTERh « ,12, i2M And N8yS ■ ,12) 

22 CONTINUE 
21 continue 

ccrccpcccccccccccrcccccccccccccccccccrccccccccccc'pccrcccccccccccccccccc 
cc cc 

cc USING the upper solid region surface control FUNC"lrtN, COMPUTE CC 

cc the upper solid region Thermal distribuIion anq interface cc 

CC GRADIENT. CC 

cc cc 

ccpcccccccccccr.ccrcccccccccccccccccccrccccceccccccrccpcccccccccccccccccc 

ICA3ES3 

WRITEC6.10) 

10 FCIRMAT(jhi,B2X,22hu PPER SOLfO) 

CALL INPUT 


Figure C-2. Computer Code List for Problems 
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on 31 MTpP ,HINTcPM,H»XTeRM,riFL^C*’»^ 

MTEPM»MAXTEPH*MtNTERM»MTCR 

OP 32 NSV «mINNSYS,M*XNSYS,OCLNPYS 

ORIGINAL FACE IS 

ir(MN.LT.N8YS5C0 TO 32 OF POOR QUALITY 

CALL S0LI03 

IFdQPTinN.FQ.O) CriTOBO 
CALL LINSPCW(XMTN) 

CA1.L rUNf.CXMlN.RMtN) 
ir(XMiN,r.T.r,.o.4Nn.5HiN,LC.o,o) cotopo 
xHiNsioonoo'.o 

90 XPl»i«XOfXMIN 

80 continue 

CALL «ELT 

ccpcccccccccccrccrcccccccccccccccccccrcccccrccccccrccrcccrcccccccccccccc 
cc cc 

CC OFTpRMINE RfLaTIVF OirrE«CNCE fipTWEEN PEOUlRfO iiPPER SPLIO CC 

cc RFCION INTepFacE GPaoIENT anO TmE INTERFACE fiPAnlENT PCSUL'INC cC 

cc from USE OF The UPPER SOLID nfCtON SURFACE CONTROL FUNCTIOf: cC 

cc . cc 

ccrccrcrccccccrccrcccccccccccccccccccrccccccccccccrccrcccccccccccccccccc 

CALL ERROR 

NPITE (6| 1201MVERM,N3Y3 

«20 FORMaTC//» IH ,50X.12mFOR mTERm , ,12, 12^ »no nsy5 ■ ,12) 

32 CONTINUE 

31 CONTINUE 

STOP 
END 

SD0ROUTIME i'ELT 

c 

C 

real Jl ,.IIuaH,MNBSJO 

COMMQN/Ci/R| ANDr2O),JlC20),JlLAN(20) 

COMmon/ci 0/ASCRTPt20),BSCRIPC2fl' 
coMMON/PFAOi /p,ntern,nsum,xo,xn,:;grio,nr 
COMHON/CS/R fl01),PSl(20, ll31),80,ti;20) 

CnHMON/c20/AC50n),B(S00),CC500) ,Q(500) ,V(50O) «0rTA(?05J ,CAM**A(505) 
C0MMQN/C2l/THET‘Bf 20,505), TMOLnM0l),T(l0),CRAOYi; no t!,GRACxO(lftl) 
CnMP0N/C22/rCA8E,TMCLT(3) 

CnM-40N/C23/r.RA02Cl01),COA03(l01' 

CMARaCTER*17 PIS,ALPmA(6> 

CNARACTER*1A 3TARS»>5TAR(6) 

data RiS/'R» '/.stars/'*****************'/ 

DO 207 iPl.A 
alPna(L)«RIS 
STARfL)*STARS 
?07 continue 

ccrccccrccccccrccrcccccccccccccccccccfccccccccccccrccfcccccccccccccccccc 


CC cc 
cc RLAMO(M)*ROOT of jo BESSEL FCN CC 
cc J1 (M)*J1 (RLaNOC*) ) “here J1 IS *E3SCL FCN cc 
cc J1LAM(M)«J1 (N)/Rla« 0 (H) cc 

cc ... ... cc 


ccrccccrccccccnccrcccrcccccccccccccccrccccccccccccrccrcccccccccccccccrcc 

RI.AMDC 11*210008255577 
RLAMOC 2)«5,520n7E 1 103 
RUANOC 3)*«. 6537279129 
RI.ahoC 0)«i 1 ,7915300391 
RLAMOC 5)«10, 9309177086 
RIAMOC 6)»18,071063R679 
RI.AMDC 71*21 ,2116366299 
RLAMOC 8)*2a, 3520715308 
RLAMOC 9)*27, 0930791320 
RLAMOC 101*30,6306060684 
RLAMOCl 1)*33, 7758202136 
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BLAMOOaiaSA.oiToowSSJ? 
RLAMO(l3^>4n,OSAa2S7646 
RI.AM0( J4)»4T,19'»7'»17J32 
RI.AH0( tS)a4Ai.34UA«371 7 
»LA“0f 16)i4«».4«2609efl74 
Rl AH0n7)aS2.624fl5iaau 
«l.AM0(i4)a59, 76551 (17550 
5LAM0n9)aSA,906983526l 
5LAMO(20)a62. 048464 1902 

J1( naO, 5191074973 
Jl( 2)a«fl. 3402648065 
JU 3)a0. 2714522999 
Jt( 4}a-0. 2324558314 
Jl( S)a0. 2065464331 
Jl( 6)9-0.1877288030 
JlC 7)90.1732658942 
Jl{ 8)a-0. 161701550’ 

Jt( 5)a0. 1521812138 
J1 Cl0)a>0, 1441659777 
Jl (11)30.1372969434 
Jl (I2)a-fl. 1113246267 
Jl (13)30.1260694971 
Jt(ia)3-n.i?i3986248 
Jt (15)30.1172111989 
Jl (16)3-0.1134291926 
Jl(17)a0. 1094911430 
Jl(l8)a-o. 1068478883 
Jl (19)90.1039595729 
Jl (20)3-0,1012934989 

on 555 l9l,?0 

JlLAM(T)»,rl{Il/«LArtOCl) 
SOJHDaJf (l)*Jl(n 
955 cnNTi*Jue 


ORIGINAL PAGE !8 
OF POOR QUALITY 


ccrccccrccccccrccrcccccccccccccccccrcrccccccccccccrccrcccccccccccccccccc 


cc cc 

CC FIND COSFS for RESSEL £XPAti3lCtN8 Of A(R)-A(1) AmO B(R)-B(1J CC 
CC 8FE EQUATI0N3 (2.2,17) ANO (2.2'.iS) OF FiNal REPORT CC 
CC • CC 


ccrccccrcccccccccrcccccccccccccccccccrccccccccccccrccrcccccccccccccccccc 

CACL AFC(1.0<aOFi) 

CALL BFCn.fl,BOFl) 

on 20 131,101 
R(nBa-n*o.oi 
RMOL09RCn 
call AFC(RH0Ln»AN3) 

A(I)sANS«tOFl 
CALL 8FC(RH0Ln,AN3) 

B(I)sANS-AOF| 

20 CONTINUE 

CALL C0eF3(R»A,10t,20,A3CniP) 

CA|.l CoCF3(RfB, 101,20,83CRIP) 

ccrccccrcccccccccrcccccccccccccccccccrccccccccccccrccfcccecccccccccccccc 


CC SOLVE FOR THETi BaR OF EOUATIoNR (2.2, l’) BV SOi VINC Tf-E CC 

CC TRIoIasONaL 3VSTER (2,2.20) - SfE FIN*l REPORT CC 

CC CC 

cccccccccccccccccrcccrcccccccccccccccrcccccccccrccrccfcccccccccccccccrcc 

OX3(XN-XO)/NCRIO 


0X2a0X*0X 

LaNRRXD*! 

DO 556 Mal,MSUH 
on 40 lal.L 

A(T Bl,p40X*P/2,0 
B(I)a«2.0-OX2*RLARO(M)9RLAMn(H) 


C(I)al.O*OX*P/2,0 


XaXO'»I*OX 


Figure C-2. Computer Code List for Problems 
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CALL RBa»(M,X.ANS) 

OC)aOXP*ANS 

ao cnNTiMiJC 

Dt .0*0X*P/2,0)»ACCRTr(H)«SOJl 

0(U»OfL3»( 1 .0-f>X*P'?.0)*82CR?r(M)»1OJl (M) » 0 .^ 

CALL TRtnAOa) 
on so I*2,NCRtO 
IIsI-l 

TH£TAB(M,n«V(in 

50 continue 

N3TOP«nCRTO+J 

THETAB(H, t )iASCPIP<'*)»80J1 (m)/2.0 
THCTAStM.NSTOPJaRSCRIRCIJ^GOJA (H)/2.0 
S56 CrUTINUC 

snTnoRR 
099 CONTINUE 
0»31 ,0/NR 

nbstqpsnr+1 OR!GII\iA> Pit"'' 

on 60 I«l,NRSTOP r\e 1 ’ “ 

R(n = (i-n* 0 R POOR QUALITY 

on 65 M»1,MSUH 

VARaR(I)»RLAHO(H) 

CALL JfUVAR.Y) 

o5 continue 
oO CONTINUE 

c 

ccpccccrccccccrccrcccrcccccccccccccccrccccccccccccrccrcccrcccccccccccccc 
cc cc 

CC PRINT temperatures CC 

cc cc 

ccrccccrcccccccccrcccccccccccccccccccrccccccccccccrccrcccccccccccccccccc 

50 TO (2), 2?, 23) ICASE 

21 WRITEC6.983) 

083 format C JHI ,5flX, tSHl' e L T ZONE) 

C0T02a 

22 hRITE(6.30) 

30 FnRMAT(lHl,52x,22HL ONER nOLIO 

GOTO20 

25 NBITE(6.10) 

10 FORMaTC 1MI»52X,22HU PPER pOliD) 

2« CONTINUE 

nRITE(6,70) 

TO FORmatiIH ,a5x.S0HT ENPEi*AtuRE OIsTrISUTIOn 

• ) 

IFLAG-O 

mRIGHToR 

MLEFToi 

I8n CONTINUE 

IFCNRSTOP.LF.MRICHTlIFLACal 

mricht*minO(nrstdp»nrichT) 

NRITE(6. 190) (R(J).JPPt£FT,«RICHP) 

*90 FORMAT!/////, IH .17X,8(F12.8,5X'; 

nRITE( 6.267)' (ALPMA(L).L»l,MRTr,HT) 

?67 F0RMAT(1H*, |7X,6A17) 

nRITE(6,268) ,STaR(L),L« 1,MRICHTJ 
;>68 FORMATHHO, l5X,fcAl7) 

DO 200 1>1,nSTCR 
lSKIP«l-t 

IHOLOoCISKIPT.OOOOni )/lO,0 
XHOLOadSKlP/lO.Ol'lHOLO 
IF(XHOLD.ST',0.005) COTO20O 
IlaNSTOPti-I 
XoxOT(II-l).OX 

on 202 JaMLEFT.MPiGHT 

ccpccrcrccccccrccpcccccccccccccccccccrcccccrccccccrccrcccrccc?cccccccrcc 
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cc 

cc 

cc 

cc 


oFTEHHiNr tfmpe»*tu»»e *t (x.nf.»<3 
5 FC EQUATTOn ( 2,2, 1«> CF FX*^>^L *>£POlfT 


ORIGINAL PAGE JS 
OF POOR QUALITY 


cc 

cc 

cc 

cc 


20/1 


202 

?to 

20 A 


Zen 


ccrcccccccccccrctrcccrcccccccccccccrcrccccccccccccrccrcccrcccccccccccccc 

THOLn(JliO.i) 

DO 204 Mil,MSUM 

THnLOfJjiTHOLOCJJ^Z.O'PSl'/t;* J)*THkTA8(H» Ii)/3(JJ1 (•*) 

CONTTNUF 
CALL HFC(y,AN3) 

THOLn(J)iTMnuO(JJ*AN8 
CONT!NUF 

white (6.21 0 )X, (TmOLOCJ)»J«M( EFT.HPIGhT) 

FORMATCSM X«.F1«,6.5h • »«>{EIW,8,2X)) 

CONTINUE 

IF(IFLAC,r 0 .nCO TO 220 
MHICHT«MRtCHT *4 
MLEFT»MLEFT +6 
cn TO jBo 

CONTINUE 

ccrccccrcccccccccrcccccccccccccccccccccccccccccfcccccrcccccccccccccccccc 
cc cc 

cc compute thermal gradients at XmVO ano xn cc 

cc .... cc 

ccrccrcrccccccrccfcccccccccccccccccccrccccccccccccrccrcccrcccccccccccccc 

CO TO (3i,3p,331 ICASE 
MRITE( 6 .Ra 3 l 
C0T034 
WBITEC 6 . 30 ) 

G0T034 
w«ITE( 6 , 10 ) 
continue 

wRITE( 6,71) 

FORmathh .flTx.TSM TmeRn*i GRaoIE 
wRITE(6«'2)xO,XN 

F0RMAT(//’/,44X, IhR.SX* I IHGOAO' aT Xa.PlO'.S, 14H 

2 .//) 


3t 

32 

33 
3« 

71 


72 


i; t S) 

GRaD, at X«»F10,S 


on 230 I»l, ioi 

Rf I)»(l-n*rt.0l 

DO 240 Mai,HSUM 
VaR«R(I)*RLamO(M) 
CALL JO(VaR.t) 
PSKM.D.y 

24 P continue 
230 CONTINUE 


ODa«OX 

00 250 lat,;oi 
OO 260 Jat(5 
T (J)aO.rt 
DO 270 Ma],MSUM 

T(.I)aT{J)A2,0«PSI(M,I)aTHpTAB(M,J)/SQJl (M) 

27/1 CONTINUE 

XaxO*(Jal)aOi( 

call HFc(X.AnS) 

T(J)aT(.!)y*NS 
26/1 CONTINUE 

on 280 JaA.lO 
T(J)aO,/» 

JMOLnaNSTOP»lOaJ 

00 26 a Mat,H9UM 

T(J)aT<J)*2,0«P8I(M,I)*THFTA9(H,JMOLO)/SOjj(M) 
26" CONTINUE 

XaXO«(JHOl.O«t)*OX 
CALL HPC(X.ANS) 

T(J)aT(J)«AN3 
28" CONTINUE 
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ccpccrcrccccccrccrcccrccccccccccccLccrccccccccccccrccrcccccccccccccccccc 
cc cc 

CC *PPP0X1MaTE ThEPMal GPAOIEf.TS *t CYLl^iOeP E'^ns cc 

cc - SEE EQU*TtQn5 ( 2 . 2 . 2 a) *(.0 ( 3 *. 2 . 22 ) OE pjnaL pCPCpT CC 

CC CC 

ccrccpcrccccccrccrcccrcccccccccccccccrcccccccccccrrccrcccccccccccccccccc 

C»*nXN(l)j(.l.T (tt)*U«T(7)-36*T>a)»<i«.T(9)-2S«Tf 10))/(12*00) 
GPAnX0m«(-3*T(S)*l6*T(«)-36«T>3)*fl«r»T(2)-2S»Tr 1 n/(l2*0);) 
rsKrP*i-i 

IHOLOadSKIP+.OOOOOt )/10,0 
XMOLOa(ISKIP/ 10 . 0 )-lHOLO 
rPCXHOLO.GT'.O.onSl GQT0251 

«PiTE(6.3ooiP(n,cR*3xo(n,cruf>yN(i) ORIGINAL PAG'?’ 

FOPMATOh ,39X,P«.6.3X<EJ7.«,7X.£17.9) 

coNTiMue OF POOR QUALliif 

iFdCASE.Nc'.ncn to 2?o 

GPAn2(I)aGRA0X0(I) 

GPAn3(I)3GRA0XN(n 
25n CONTINUE 


300 

25 * 


RETURN 

ENO 

ccrccrcr.ccccc->.ccrcccccccccccccccccccrcccccccccccrrccrcccrcccccccccccccc 


cc cc 
cc This subroutine aPPRCXIPaTES fflv finite DIFFERENCE) G EAR Pf CC 
cc equation (a'.a.iA) cf final rcpopt cc 
cc cc 


ccrccccrcccccccccrcccrcccccccccccccccrcccccrcccccrrccrcccccccccccccccccc 
subroutine G8Afi(H,X»ANS) 
real JI.JILaN 

COHnCN/CWRi A“0(20) » J1 (20),J1lah{20) 

cnMMON/OEAOt/P,HTER«.MSUM,Xo,Xn.;iGRID,NR 

EPSUONaO.Ol 

XlaX-EPSl.QN 

X?xX»EPS1.0N 

CALL HFC(X.anS) 


CALL HFCdl,4NSn 
CALL hFC(X2.anS2) 

G«T (ANS2-ANSi)/(2.0*EPSLa!.) 

G«C-(ANS2*ANSl«2.0«*NS)/CEPSLnN*EPSLON) 

ANSaG*Jll.AMfP) 


return 

ccrccccrcccccccccrcccccccccccccccccccrccccccccccccrccrcccccccccccccccccc 


cc cc 
■cc PURPOSE cc 
cc - provide Input data for snFTw*pe cc 
cc - SEE appendix a. 3 FOP DETAILS CC 
cc cc 


ccrccccrccccccrccrcccrcccccccccccccccrccccccccccccrccrcccccccccccccccccc 

SUBROUTINE INPUT 
INTeCER PELTCPn.OELNSTS 
CnHHON/r22/TCA3E.TMELT(3) 

CnHM0N/C?a/PK3.RKL,»L»NST3 

CnMMON/PFAO(/P.NTC»N.NSUH,XO.*N,fICRID,NR 

COMHON/FlXPT/IOPTTON,XMIN,GHlN,rLIP 

C0MMQN/C31/THFC 

COMHON/CS 2 /XO r 1 00 ) » TO(joo),ci(«,ioo),H 

COMMON/CIO/nINTFPm, naXTERM.OClTfKM.hIhnSyS, HAXNsYS.DELf.'SYS 
DIMENSION C(a,t 00 ) 

DIMENSION XmOLO(100),yhOLO{JOC) 

equivalence (CJ(1.1)»C (1,1)) 
pPITE(6,5) 

5 FORMAT!//, IH ,S*X,20mI NPUT OAT a) 

IFdCASE.NE.n GOT060 

ccrccccccccccccccpcccccccccccccccccccrccccccccccccrccrcccccccccccccccccc 
cc . . . cc 


Figure C-2. Compucer Code List for Problens 
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CC INPUT mEI.T 70N6 SUBFiCE TEtl*. (TTSTHIBUrlO** IN A PAT* SCT CC 

CC rnRMAT FOR USE TM A CUBIC 3PI.INF CC 

CC .. . CC 

ccrcccccccccccrccrcccccccccccccccccccrccccccccccccrccfcccccccccccccccccc 

RFAP(5,«ft)lMFC 

80 FORNATCIz) . rvB/'C ■« 

iF(iHFc,Nc.i j COT060 ORIGINAL PAiaC. l» 

read(s,««mi) h ^ pnOR OUAUTY 

««»<» format (IS) Kuun 

MRITE(6.S08m 

•\8 FnRM*T(//////,RSM the SUBFaCC Tf^Pe^ATURE DISTRIBUTION 13 APPR0*|N 

PATEf) BV the cubic spline TmRCUgm The following, T fl,25N (X, TEMPT 0 
•ata points, ///,37h X surface Temp.) 

DO 3.2 lai.M 
REAO(5,??)xn(n»V0(I) 

22 FnRMAT(2F20'.lO) 

MRITE(6.3«))(0(l).TP(n 
3« FORM*T(2E20.10) 

32 continue 
CALL COFGEN 

60 continue 

ccrcccccccccccrccrcccrcccccccccccccccfccccccccccccrccrcccccccccccccccccc 

CC CC 

CC INPUT hEI.T 70NC parameters CC 

CC CC 

ccrccccrccccccrccrcccccccccccccccccccrccccccccrcccrccrcccccccccccccccccc 

RFAO(S, in)P,XO.XW,MTERM,HaUM,NGPlD,NR 
to FORMAT(3Fl0.5,allO) 

«RITE(6.R531 

aS3 format (/, tH . |a*. IMP,20X,2H*0 ,irX.2MXn. T*. 5MMTEpli,6x,aHM3U'’,5X.SMN 
«GRID,6X,2MNR) 

WRITE(6.20)P,X0.XN,HTERM,M5Ut!,Ncf:iD,NR 
20 FORMaTOH ,3E2O.10.«It0) 

IFdCASE.EO'.nRETURN 

ccpccccccccccccccpcccrcccccccccccccccrcccccccccrccrccrcccccccccccccccccc 

CC CC 

CC INPUT material CONOUCTIVITIES CC 

CC . . . CC 

ccrccccfcccccccccrcccrcccccccccccccccrccccccccccccrccrcccccccccccccccccc 

REArXS.OnlRKS.RKL.RL.NSTS 
RO FnRMAT(3E20‘,l0.Il0) 

WB1TE(6,888) 

abb FORMAT!//, IH ,tOX,3MRKS,tTX.3MRFL»jTX,2MpL) 

ccrcccccccccccrccrcccccccccccccccccpcrccccccccccccrccrcccrcccccccccccccc 

CC CC 

CC INPUT material conductivities CC 

CC .. . CC 

ccrccccrcccccccccrcccccccccccccccccccrccccccccccccrccrcccccccccccccccccc 

WRITE(6,aO) RKS,RKL»RL ‘ 

00 FoRMAtUH ,3E20,t0) 

ccrccccccccccccccrccccccccccccccccccccccccccccccccrccrcccccccccccccccccc 

CC CC 

CC INPUT CASE I IMITS CC 

CC 


CC .. . CC 

ccrccccrccccccrccrcccrcccccccccccccccccccccccccccrccccccccccccccccccccc 

REA0(5,21)MTNTe»M,MAXTERM,DELTFor,MlNN3TS,MAXN3Y5,0ELNCT3 

2t FORhaT(BXiO) 

nRITECp.tori 

■»R9 FORMAT!//, JH ,9X,S7MHaXTeRM mtNTeRM OELTFRM m*XN3T3 HINNSY 
-S OELNSYST 

WRITE !6*JB)NAXTERM,MlNTEBM,nCLTrRM,MAXN3yS,MiNuj,y8,oELN3Y3 
18 FORMATdH ,SX.A!5X,IS)) 

REAO!S,SO)lnPTION,CLlP 
30 FORMAT!I10,F10.5) 

return 

ENO 


Figure C-2. Computer Code List for Problems 
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-V r- J.. ' ' -ff- '''■ ir'"-' • ’ju.'tv’.i .-.'Li. . 


> 



cCrCCccrcccccccccfCCcrccccccccccLCCrcrccccccccccccrccrccccCCCCCCCcCCcccc 


cc PURPOSES CC 
CC • PROVIDE USES ENTRY OF FOUCTTOmAL FO»m PF nfCT TOnf SL'RF*CE CC 
CC tfmp. distrtbution cc 
CC • evaluate solid regions SUREacf control FUNCTlnl4S CC 
cc • MnoiEY SOI ID Regions surface rONTROL euncttOns ustng iop'iQn cc 

cc A*'0 CLIP AS OETaILEO in APPCf.ftIv A.J CC 


ccrccrcrcrccccrccrcccrcccctccccccccccrccccccccccccrccfcccrcccccccccccccc 
subroutine hFC(*.*NS) 

CPmnON/CR/CPEF(201,Rh 

cnMMOH/C22/TCASE.TH£LT(3) 

CnNNQN/Cp4/OKS.RKL>RLtN3YS 


CrNMOH/RCiOl /P,NTERN,N8UH,X0,*N,!JG»Tn,NR 
COMMON/EIXPT/TOPTION,XHIN,GHlM,rLIP 
COMmON/cSO/TCKCUT 
COHMON/C?6/rPOLY(20) 

CnMHON/C32/YOnnO), YO( 100 5,Clf4, 1001 ,- 

COHMON/CSl/TMEC 

OIHENSiCN C'a,t^01 

DIMENSION 2f20) 

equivalence f Cl (’.i)*c(i.nj 
rFClHEC.EQ.I ) 6OTO60 

lEdCASE.Ec' 2)COTn20 
IE(IC*SE.EQ.3; COTORO 


ORIGINAL 

OF POOR QUALITY 


ccrccccrcccccccccrcccrcccccccccccccccrcccccrccccccrccrcccccccccccccccccc 


place user SUPPLIED MELT ZONE SiiPEaCE TEmo NpRE 


cc 

cc 


:c 
:c 

^ . cc 

ccrcccccccccccrccrcccccccccccccccccccrccccccccccccrccrcccccccccccccccccc 
return 


RO C*LL SPlINETY. ansi 
return 
20 continue 


ANSsO.O 

OP to K»1,NSY3 

2fR)a(l-«)«rxN«xi 

RMOLOro.d 

IE(7(K),f,T.-2S0.0) RnPLOrEXPCZCkI) 

AN3«ANS*rOEE(Kl*«NOLO 

lO continue 

IECIOPTIPN.fQ.o) G0T0«»5 
IECK.GE.XMIN) cOTOrS 

ccrccrcccccccccccrcccccccccccccccccccrcccccccccccccccfcccccccccccccccccc 

cc 

;c modify loner SOLID “CCICN SUPr*rE cONTROl EUNCTtCN aS CEEII.’eO cc 

:c 8Y Value of ioption cc 

, . cc 

ccrccccrccccccrccrcccccccccccccccccccrccccccccccccrccrcccccccccccccccccc 

AN3«CMIN,(2.ICPTI0N)Aa0PTI0l—l>«AMlNl (*nS,ClIP) 

«5 return 
40 continue 
AN3«0.0 

DC 50 ml.NSYS 
ZfKls. t>Kl«rX»XO] 

Rh0LD«0,0 . 

IE(ZCK).r.T..250.0) RMPLD»EXP(Z(if)) 
anS«anSaCOEE(K1 •RMOLO 
50 continue 

lEdOPTinN.EO'.Ol G0T055 
IEIX.LE.XMJN) COT055 

ccrccccrcccccccccrccccccccccccccccccccccccccccccccrccfcccccccccccccccccc 

cc 

:c MODIFY upper solid region SUREArr control EUNCTtCN aS CErir:EO CC 
:c BY value gr ioption _ cc 


Figure C-2. 
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ORIGINAL PAGE 5S 
r>c POOR QUALITY 


0 - 


cc cc 

ccrccccrccccccrccrcccrcccccccccccccccrccccccccccccrccrcccrcccccccccccccc 

55 RFTIIRN 

end 

ccrcccccccccccrccrcccccccccccccccccccrccccccccccccrccrcccccccccccccccccc 
cc cc 

cc THIS susroutine supplies p*oi*l TeHPc»»TuPE niSTdiBuTirM cc 

cc ON LOWE" ENn OF C»LIN0E« - SEE FGU*TlnM (?.2,2) OF FINAL REPORT CC 

cc ... cc 

ccrccccrccccccrccrcccFcrcccccccccCLCcrcccrcrcCccccrccpcccccccccccccCcccc 

SUBROUTINE aFC(R.anS) 

CnHNQN/RFADl/P.NTFRM,HSUH.X0,XN,U6RI0.NR 

ccrcccccccccccrccrcccrcccccccccccccccrcccccccccccccccrcccccccccccccccccc 
cc cc 

cc user SUPPLIFO lower E* TEUPERa^URe AfR> cc 

cc cc 

ccrccccrccccccrccpccccccccccccccccccccccccccccccccrccccccccccccccccccccc 

CALL HFC(XO.ANS) 

return 

END 

ccrccccrcccccccccrcccccccccccccccccccrccccccccccccrccrcccrcccccccccccccc 
cc cc 

cc this subroutine supplies raoial tehperatore niSTRiBuTirN cc 

cc ON UPPER ENn OF CTLInoeR - SEE fQUaTION (2,2.2) CF FINAL REPORT CC 

cc cc 

ccrccccccccccccccrcccrccccccccccccccccccccccccccccrccfcccccccccccccccccc 

subroutine rFC(R.anS) 

CPMHON/RFAOl/P,HTERH,"SUH,XO,XN,fJCRIO.NR 

ccrccccrccccccrccrccccccccccccccLCCccrcccccccCccccrccccccccccccccccccccc 
cc cc 

cc USER SUPPLIED UPPF» END TEJ'.PERaxURE BfR) CC 

cc cc 

ccrccccccccccccccrcccrcccccccccccccccrccccccccccccrccrcccccccccccccccccc 

CALL HFCfXN.ANS) 
return 

END 

ccrccccrccccccrccrcccrcccccccccccccccrccccccccccccrccrcccccccccccccceccc 
cc cc 

cc THIS subroutine pits RESSCL SFRtCS TO 0 *TA BX least soi.are; cc 

cc nfThQd • SEF equations (2,2.17) , (2.2. ‘(O anO ■f2,2‘.23) CC 

cc OF final report cc 

cc cc 

ccrccccccccccccccrcccrcccccccccccccrcfcccccccccccccccfcccccccccccccccccc 
subroutine rOEF3(R,x,NR,NCCEr,crCP) 

InTcQER nR.ncOEF 

REAL F.ROOt ),mon,COEF(20),M«(a60) 

external f 

ccrcccccccccccrccrcccccccccccccccccccrccccccccccccrccrcccccccccccccccccc 
cc PC 

cc user supplied least squares ueTmOO FOllOwS here to oETCRMIUE cc 

cc The coefficients of EQUaTIOnC (7.2,17) *nO (2,2'ib)’, the cc 

cc subroutine IFLSQ below is the IhSL least sOuaRES function riT cc 

cc routine cc 

cc , . . cc 

ccrccccpccccccrccccccccccccccccccccccrccccccccctrccrccrcccrcccccccccccccc 

CALL IFLSO(F,R,t,nR,cOCF,NCOCF,wK,IER) 

IF (IER.E0.J7R. OR. lER.EO.l 30 )URlTe(R, In) 

•0 FnRMAT(5NH TERHTNAL ERROR IN LEaST SOUaReS "FTmpiO, SUBROUTINE COEFS 
n 

return 

END 

ccrccccrcccccccccrcccrcccccccccccccccrccccrcccccccrccfcccccccccccccccccc 
cc cc 

cc Function F nsco in subroutine coEFs F(n.»)«jo(laMOa,n)*R) cc 

cc cc 


Figure C-2. Computer Code List for Problems 
Pi-1 and Pl-2 (Cont) 
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ccrccccrccccccrccrcccrccccccccccLCcccrCccccccfccccrccccccccccccccccccccc 
REAL ruNrTiriN f(n.R) 

IHTFCER N 
"FAL R,Rl.AHn{20) 

RLAMOC n»2' ao«R2'55'i77 
RlAMOt 2>aS.520«78no5 
RLAM0(3)a8.h53727Rl29 
»LAMO( a) = lt.7<»153««3Rl 

f»LAMO( S'»«ia.i»3nRl77oa6 • 

RIAHO; 61.1A.0710*3R67<» qc pQQ- ..LmLI iV 

RLAMOf 7is2t ,21 163662RR <->r r v 

RLAhO( «)s2a,35?«7t530a 

RLAmOC R^«^7,aR3a7Rl320 

Rl AM0(10)a3a.634«06O68R 

R» AHD(l I)a33.7788202l36 

RL AMO a2)s3<,.R 17(1083537 

RL AMO (iJ) ■4(1,0588257648 

Rt.AM0(141a43, 10070171 32 

RLAHOf 151=46.3411883717 

R|. AMO (1 61 »40. 0826008474 

RLAM0(l7)a52.62405l8al 1 

Rl AMO(i8)=5S.7e55lO’550 

RL AMO (101=58.0060830261 

RL AMO (2 01 *62. 048460 1002 

X*fll.AMO(N)*0 

CALL JQ(X«Y1 
F*Y 

RfTl-RN 

EMO 

ccrccccrccccccrccrccccccccccccccccr.rcrccccccccccccrccrcccccccccccccccccc 
cc , . cc 

CC THIS SUSRCutINE computes the jo 3ESSEL T=Jo(X) CC 

cc . . cc 

ccrccccrccccccrccrcccccccccccccccccccrccccccccccccrccrcccrcccccccccccccc 
SUeRQUTINE .tO(X,Y5 

ccrccccrccccccrccrcccrcccccccccccccccrccccccccccccrccrcccccccccccccccccc 
cc cc 

cc use® SUPPLIED JO FC*<. PLACF.D MfoE. The JmSl routine MHpSJo . cc 
cc IS illustrated 8CL0M CC 

cc . . cc 

ccrccccrccccccrccrcccccccccccccccccccrccccccccccccrccccccccccccccccccccc 
REAL MM0SJO 
Y»MMeSJO(X,TER) 

0''7URN 

-NO 

ccrccccrcccccccccrcccccccccccccccccccrccccccccccccrccfcccccccccccccccccc 

cc . cc 

cc SU6R0UTIHE FOR SOL^INC A SYSTpH OF LINEAR SIMULTANEOUS CC 

cc EOUaTIONS H*VINC a TOIOIaCOnal rOCFFIclENT maTRtX. dIaCOHALS CC 

cc ARE STORfo TN the arrays a, G, iND C. THg cOmPUtED CC 

cc solution vector v(i) , , , v(L) IS stored in The array v, cc 

cc . . . cc 

ccrccccccccccccccrcccccccccccccccccccrcccccccccccccccrcccccccccccccccccc 
SUBROL'VINE TRIOAC(L) 

cnMM3N/C2O/AC5OO),0(5OO),CC5OO),O(5OO),V(5on).erTA (5053 , gamma (505) 

ccrccrc. : rccccrccrcccccccccccccccccccrccccccccccccrccfcccccccccccccccccc 
cc cc 

cc compute intermediate arrays OETa and gamma cc 

c .. . cc 

career crccccccccccccccccccccccccccccccccccccccccccrccrcccccccccccccccccc 

aET4(n«fl(i) 

CAMMA(1}bO(1)/BETa(1) 

IFPl«2 

00 1 lalFPUL 

BCTAdiism-Am-cci-n/ecTAi-i-i) 
GAMM4(n«fO{U-*(I)«CAMMtt(I.p}/BET4(I ) 




Figure C-2. Computer Code Lise for Problems 
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1 continue 

ccrccccrccccccrccrcccrcccccccccccccccrccccccccccccrccrcccccccccccccccccc 
cc cc 

CC COMPUTE FIN*|. S'iLN, VBECTOR V CC 

cc CC 

ccrccrccccccccrccrcccccccccccccccccrcrcccccccccccccccrcccccccccccccccrcc 

VrL)«CAHHi(| ] 

L*3T*t-l 
on 2 K«t,LAST 

v(i)!lcAMMA{n-c(i)*T(i»i)/nFTicr) ORIGINAL PAGE fS 

2 CONTINUE OF POOR QUALITY 

t.ETURN 

END 

ccrcccccccccccrccrcccccccccccccccccccccccccccccccccccrcccccccccccccccccc 
cc cc 

cc THIS SUBBOUTIne OETEBnINES T[if lO“EP SCltfO BfClnUS SU»r*CE cc 

cc CONTROL FUNCTION AS OUTLINCO IN CMAPTFR 3 OF FImAL REPORT, CC 

cc cc 

ccrccccccccccccccrcccccccccccccccccccrccccccccccccrccfcccccccccccccccccc 
subroutine sol 102 

real J1».MLaM,MN3SJ0 

COHH0N/Ci/fl| AM0(20).JI (20),JILA“(20) 

COMMON/CS/Rnon»RSI(20,101),3fl.»K20) 

COMMON /CQ/COEF (20 )»Rh 
COMMON/C22/TCASE,TmElT(3) 

CnMMON/C23/ORAO?not },CRAO3d0p 
COMMON/C2a/RKS.'!KL.RL.N3Y3 

common/reaoi/p*ntern»msum,xo,)(n.:,’grio,nR 

C0MM0N/CS3/F(«) 

CnMM0N/C76/TFLACa, IFLAG3 
COMMON/CT7/AHAT{20»20),RM3(20) 

oi“ension Ynon»c(«).iwK(20),Ni^(95o) 

OIHENSION A|PMA(20) 

IFCIFLAC2.EO.O) COTO120 
MRITE(6.aA4I 

flAA format(ihi,3ix.aohl oner solio Thermal cra 
• 0 I E N T S) 

NRITE(6,333I 

^33 FORMaT(///, 1M ,52X, IhR.27X,4MCR*0) 

00 30 JJaifiOl 
J»102-JJ 

ccrccccrccccccrccrcccrcccccccccccccccrcccccccccccccccfcccrcccccccccccccc 
cc cc 

cc OFTeRMINF loner solid RECICN interface '•RAOIENT see equation cc 

cc (F24)» FTCURE 1-2 OE FINAL REPORT CC 

« . ... cc 

ccrccccccccccccccrccccccccccccccccccccccccccccccccrccccccrcccccccccccccc 

GRAn2(J)»{flKL«CRA02(J)-RL)/RKS 

RRITE(6»555’R(J)»GBA02(J) 

R5S F0RmaT( 1H ,3RX.E20,»0,lOX,E20',lii) 

Y(J)«CRAO2(j)-CRAD2(l01) 

so continue 

ccrccccrcccccccccrccccccccccccccccccccccccccccccccnccrcccccccccccccccccc 
cc cc 

cc determine coefficients in OESSEI expansion (3.0". J4) cc 

cc . . cc 

ccrccrcccccccccccrcccccccccccccccccccrccccccccccccrccrcccccccccccccccccc 

CALL COEFS(R,Y,tOI»20,COEF) 

ccrccccccccccccccrcccccccccccccccccccrccccccccccccrccfcccrccficcccccccccc 
cc cc 

cc determine matrix and vector Clements as oefineo in cc 

cc equation (3',0.24) cc 

cc cc 

ccrcccccccccccrccrcccccccccccccccccrcrcccccccccccccccrccccccccccrccccccc 

DO AO M«i,MTERN 


Figure C-2. Computer Code List for Problems 
Pl-1 and Pl-2 (Cont) 
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L»H*2 

rhs(i.)«rl*hp(m)*ji (H).cor.F(M)/2,o 
l>0 CDNTIMUC 

anFiaCRin2( I on 

*nri»THCl.TCTCiSC) 
on 103 '^■i.mteb*' 

ALPHA (M)»P*P*a.n«BL*MO(H)«Rl.AHP>K) 

ALPMA(H)»(P-S0BTt4LPHA(H)))/2,0 

L»H*2 

RHS(L)»nM8(l ) + {ALPHA(M).p),AOFl*60n 
RHS(L)aRHS(i )/(4LPH*(”)»(‘*-*I.PH* (M) ) ) 
*03 CRNTIMUE 

on 106 h»1,mTeRh 
on 107 K»i,NSYS 

RKHnLDaK-l 

L«Ma 2 

AHAT(L»lO»r.O/(RKMOLO-ALPHA(M) ) 

*07 CONTINUE 
*06 cnNTINUE 

DO 108 K=i,nSVS 

amat(i,k) = i‘.o 

*08 CONTINUE 
RHSmaAOPl 
on lOR «bI,nSY 3 
AHAT(2,<5a-fK*n 
*09 CONTINUE 

RHS(2)a-B0Fl 
*20 CONTINUE 
IFLAC2aO 
NaMTERHT2 


ORIGINAL PAGE iS 
OF POOR QUALITY 


on 111 lal,a 
EtOaO.O 
«n CONTINUE 

ccfcccccccccccrccrcccrcccccccccccccccrccccccccccccrccrcccccccccccccccrcc 

cc 

CC SOLVE THE Over POSFO LINEAR SYSteh OF eOuATTONJ £3,0.2':) It: CZ 

CC the least SOUaRFS S€nSE, The IMRL ROUTIHe LLBOF is ILLUSTRATEO cc 

cc bfloh (See remarks at the ent, of arrenoIk a, 3) .cc 

cc cc 

ccrccccrccccccrccrcccrcccccccccccccrcrccccccccccccr.ccrcccrcccccccccccccc 
CALL LLBOF(AMAT,20,M,NSYS,RH5,2A,l,0,E,CafF,?0,TWK.wF.:ER) 
RMOLOaO.n 
on 113 Kb2,NSYS 

RH0L0aRHnL0*Cfl6P (K) 


1*3 continue 

COEFUIbaOEi-RmOlo 

ccrcccccccccccrccrcccrcccccccccccccccrccccccccccccrccrcccccccccccccccccc 

cc . cc 

cc 0T3PLAY COEFFICIENTS USED IN THE EXPANSION OF ThE LOHCF CC 

CC SOLID RcriiON surface control function (See equation ( 3.0.23) ) cc 

cc . cc 

ccrccccccccccccccrcccccccccccccccccccccccccccccccrrccrcccccccccccccccccc 

MPITEf6.1«0) 

lart formaT(ihi) 

HRITE(*,,7771 , 

^77 FORHATfiMl, 19X.K3HL ONER .KOLID SURFACE CON 
•7R0L COEFFICIENTS) 

HRITE(6.90> 

90 F0RHAT(///, *H ,a9X, lHK,22X,anc(K) ) 
on 778 lal.NSYS 
HRITE(6.686) I. COEF(I) 

^86 F0RHAT(/,1H , aSX. 12* 10X.E20. 10) 

778 continue 
RETURN 
ENO 

ccfccccccccccccccrcccccccccccccccccccrccccccccccccrccrcccccccccccccccccc 
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OF POOR QUALITY 


thermal 


C R A 


CC CC 

CC this subroutine OETERHINES T(,r iiPRe" sot-ro REOInUS SURf*CC CC 

CC CONTROL function iS OUTLINCO IN CHARTER 1 OF FINAL REPORT, CC 

CC CC 

ccrccrcrccccccrccrcccrccccccccccccLccrcccccccccccrrccrcccrcccccccccccccc 

SUBROUTINE SOLins 

RFAL JIjJILaM.HHbSJO 

COMMON/CI /R| AMO<aO) » jt (20) , JlLAH(ao) 

COMHON/CS/RaoD.RSKRO, ioi),sa.u( 20 ) 

COMMON/CR/CnEF(20)#RH 

common/C22/tcasC,thElT(3) ORIGINAL PAGE FS 

COMM0N/C23/RRA0an0l),CRAD3(l0p - 

COMMON/caa/RKS.RFLfPL.NSYS 
COMH0 N/Rea01/R.hTERh.hSUM,X0,xn,NCRI0,nR 
C0MM0N/C53/F(a) 

CPMM0N/C7S/TFLAG2, 1FLAC3 
COMM0N/C77/iMAT(20,20),RH3(2flI 
OIMFNSION YnOl),C(«) , IWK (20) (950) 

OTMENSION A|.PHA(20) 

IFnFLAG3.EO.0) GOTO120 
NRITE(6,uain 

naa FORMAT(lHt,31X,60HU PP£R SOLID 
• R I C N T s) 

RR1TE(6.333) 

SS3 FORMAT!///, 1H ,52X, tHR,27X,aHGRA0) 

OO 30 Jj3l,t0i 

ccrccccccccccccccrcccccccccccccccccccrcccccccfccccrccrcccrccccccccccccr: 

CC CC 

CC OFTeRminf upper solid region TNTCRF4CF <^RiOlFNT SEE EOLATION cC 

CC (F22), FIGURE 1-2 OE FINAL REPORT CC 

CC .. . CC 

ccrccccrccccccrccrcccrcccccccccccccccrccccccccccccrccccccccccccccccccccc 

GflA03(J)«(RKL«CRAn3(J)-RL)/RKS 

h»1TE(6»5S51R(J’.6RA03U) 

955 FORhaTOH , 397 ,£ 20 ,tO,tOX,r. 2 n,lrt) 

Y(J)«0RAn3(.n«GRAO3(10l) 

SO continue 

CALL C0EF3(R,Y, t01>20«COEF) 

ccrccccccccccccccrcccr.cccccccccccccccrccccccccccccrccrcccccccccccccccccc 

CC CC 

CC OFTERHINE MaTRI* and vector CLE“GHT3 as defined in CQUATIOr.’S CC 

CC (3.0,26) • (3.0.28) CC 

CC CC 

ccrccccccccccccccrcccccccccccccccccccrcccccccccccccccrcccccccccccccccccc 

DO 80 M»i,MTERM 


L“H*2 

Rh3(l)«RLaH0(h)*J‘ (h)*cOEF<m)/(-2,0) 

so continue 

BOFia6RAn3(i01) 

AOFi«TMCLT(TCa8E) 

OO 103 ‘'■i.mtERH 

*!.PHA(M)RP«P*a,<l«RL*MO(M)«RLAHO>H) 

Al Pma(m)«(P*30RT(alPma(m) ))/{ 2’,0) 
L«Mt2 

RM8rL)RRH3(|,)A(P-ALPHA(M) ),Ar.Fl-eOFI 
Rh3(L)rRhSCU/( alpha (H)*(P-alph*(M))) 

<03 continue 

00 106 HrI.mTERh 
on 107 Kal.NSYS 

LRM*2 

AMAT(L*K)«l'.0/(-l ,0*KPALPHA(t;)) 

<07 continue 
<06 continue 

OO 108 K«1,nSy8 
AMAT(l,K)«r.0 
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original PAGP fcr 
OOOR 


loa continue 

RHjn)«AflPl 
OP 109 <«1 ,nSYS 
*NAT(2,K)al-i< 

»09 CONTINUE 

RHS(2)»B0F1 

*20 continue 

IFLAGlaO 
M«MTERma 2 
00 111 I«l*a 

Ecn«o,o 

*ll continue 

ccrcccccccccccrccrcccccccccccccccccccfcccccccccrccrccrcccccccccccccccccc 
cc cc 

cc SOLVE the over poSEO LINEAH System of EQUATIONS {3,0.26) - CC 

cc (3,0,28) FOR the COEFFICIENTS Tn BE USED iN (!.n,31). THE CC 

CC IMSl RQUTinf LLRQE Is illustr*tfo below cc 

cc .. . cc 

ccrccrcrccccccrccrcccccccccccccccccccrccccccccccccrccrcccccccccccccccccc 
CALL LLBOF(AHAT,20|M,NSYS,nM5,2»*, 1,0,E,C0 Ef.? 0, iWK.WK, lER) 
RHQLOao.n 
DO 113 Ka2,NSYS 
RHQLOaRHOLDACCEF(K) 

1*3 continue 

COEFdlaAOFi-RHnLO 

ccrcccccccccccrccrcccccccccccccccccccrccccccccccccr.ccrcccccccccccccccccc 


cc 

cc 

cc 

cc 

cc 

lOrt 


OTSPLAY rOEFFICIENTS USED IN TMF EXPANSION OF TwE UPPfO SOLID 
REGION surface CONTROL FUNCTION {SEE EDIJaTION (t,0.31) ) 


cc 
cc 

cc 

Lcccccccccccccrcccrcccccccccccccccrccccccccccccrccrcccccccccccccccccc 
wRITE( 6, taOl 
FORMAT( I •*! J 
HRlTE(6*V77l 

777 FORMATClMl/ iRXtRSHU PPER SOLID SURFACE CON 
•TRQl COEFFICIENTS) 


HfllTE(6»90) 

90 fORmaTC///, IH ,U9X« lMK»22X,aj,c(i') ) 

00 778 lal.NSYS 

W»ITC(6*n 86) I* COEFCn 

a89 FORMAT(/,iH ,a6X,I2*10X,C2a.tO) 

778 continue 
3OTO180 

I8n Continue 
return 

ENO 

ccrccccrcccccccccrcccccccccccccccccccrccccccccccccrccrcccccccccccccccccc 


cc cc 
cc PURPOSE cc 
cc - perform line search to DCTERMtHE min. pT, ON the surface cc 
cc CONTROL function, USeO IF lOPTinN al OR 2. CC 


'•crccccrccccccrccrcccrrccccccccccccccrccccccccccccrccccccccccccccccccccc 

SUBROUTINE LIN3RCH(XmIN) 

COMMON/RF A0 1 /P • MTERM , nSUM , X 0 , XN.MCRTO , nR 
DIMENSION FTBdOS) 


AaO.O 

aaXN«XO 
STOREa-l.O 
0ELBa8/2O0,n 
00 90 I»l,2n0 
XalaOELB 
CALL FUNC{X,Y) 
RMQLD«Y*8T0pC 
IF(BHOL0,LC,0’.0)CO TO 100 
90 CONTINUE 
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10ft 8aX 

ALPHiao.m 
OFLsB*i 
FTBnJal.O 
FT0(2)«2.0 
5 CflNTnuC 

BRat.O/Al.PHi 

!P(nB-2.fl)lA.IU.11 

10 Cn TO 14 

11 CONTINUE 
J.Ia2 

12 J.taJJti 

FIBfJJ)aFIBfJJ-t)*f IB(JJ-2) 
CraFIB(J.>) 

IP(CC-B6)l3«l5<tS 

13 SO TO 12 

14 WBITE(6,2) 

2 PORMATC///, lOX.SPMMUST ChANSE 

15 laO 

KKS,IJ<2 

IKaJJ-2 

BLaR-A 

Ai.LaFia(rK)*aL/Fin(Jj) 

NaA^ALL 

VaB-AUL 

CALL FUNCCW.T) 

C*LL PUNrtV.U) 

JKat 
IKalK-l 
J.taJJ.i 
00 TO lat.KK 
IP(U-T)2ft,2ft,22 
20 AaA«ALL 
BLaB-A 
WaV 

CALL FUNr(H,T) 
iLLari8<IK).8L/FI8(Jj) 

VaB-ALL 

CALL rUNC(V,U) 

ITaltl 

IXRIK-l 

J.IaJJ-l 

IF(IK-J)28»?Pi2P 

28 IKai 

29 CONTINUE 
cn TO 70 

22 BaB'ALL 

BLan>A 
VaH 

CALL FUNr(V,U) 

ALL aPia(IK)aBL/F!BCJJ} 
NaA«ALL 

CALL FUNrCW.T) 

ITaI*i 

IKaiK»l 

JJaJJ-1 

IFCIK-D30»^1»3I 

30 IKai 

31 CONTINUE 

70 continue 

EP3aO.OOlaH 

Ol.»MT£PS 

CALL FUNc(DI.,VL) 

IF(TL-T) 80.80*81 
80 CALL PUNr(B,8P) 

XNlNa(H»R)/»,o 


ORIGINAL PAGE 19 
OF POOR QUALITY 


Al PHa In SgBflnuTiwC LTNSrCM) 




oriToaT 

«l CitL FUNC(A.Ar) 
XM1Ni(W^A)/?, 0 
1 F10,a,2x,2HXa,E10.«) 
87 4CC»(W-*)/{nEl ) 

OP9 continue 


ORIGINAL PAGE fS 
OF POOR QUALITY 


ENO 

ccrccccrcccccccccrcccrcccccccccccccccrccccccccccccrccrcccccccccccccccccc 
cc cc 

CC PURPOSE cc 

cc • evaluate basis riiNcTiONS usrn in cqua'ion a.n.sj) cc 

cc . . cc 

ccrcccccccccccrccrcccccccccccccccccccrccccccccccccrccrcccccccccccccccccc 

SU8P0UTINE FUNCfX.V) 

COMMON/C9/CnEFCPO>.PH 

CnMMQN/C?a/BKS,»KL#'’L.NSYS 


T»0.0 

DO 10 Ksl ,NSTS 
Zc(l.K)*X 

IF(Z.LE.-2Sn.O)r.O TO 10 
YaY*COEF(K)*£xP(Z) 

«o continue 

RFTURN 


END 

ccrccrccccccccrccrcccrcccccccccccccccrCCccccccccccrccFcccrcccccccccccccc 


cc cc 
cc PURPOSE cc 
cc - Evaluate 1.2 diffenence between The rb“uired Solid peoions cc 
cc interface gradients and those Obtained by use of the selid cc 
cc REGIONS surface cOnTrol FUNCTIOmS, CC 
cc . . cc 


ccrccccrccccccrccrcccccccccccccccccccrccccccccccccrccf cccccccccccccccccc 

subroutine FRROR 

C0MHaN/C22/TCASE,TH£LTC3) 

COMHON/C23/fiRAOP(lonf CRAOSaop, 

CnHMON/CEl/THeTAB(20*BOS),TNrLOYlOl),T(>0).CRAOYNC'ni5,CRAOxO(lft) 

IP(ICA3E,EQ'.3) 50T050 

OXNL2«0.rt 

CXNLINPan.O 

DO 10 J«! , mi 

GxNL2«CXNL2*CRAOxN(J)*CRAOXNCJ) 

XHAGlaABSICRAOXNCJ} ) 

G*NLINF«AMaX 1 (XMACUGXNLINf) 

10 continue 

GXNL2*S0RT(r.XNL2) 

E2NUM«om 
EINFNUMafl.O 
DO 20 Kai, mi 

E2NUH8<gRAOXN(K)>CR*D2(K] )«*2.i)*C2nuh 
XHAC3aABS<GRA0AN(K)«GRA02(K}} 

EINENUHbAMAXI (XMAG 3>EINFNUfl) 

20 continue 

E?NI/HaSQRT{F2NUN) 

ERRL2«E2NUM/GXNL2 

ERRLINEbEINPNUH/CXNLINF 

GOTnao 

50 GX0L2a0m 
GXOLINEan.O 
DO 80 JJal,iOi 

6XOL2aCXnL2*GRADXO(JJ)*GRAOXOYJ.T} 

XHA61aABS(CRAOXO(Jjn 
GXOLlNEPAMAXl (XMA61 «GX0LlNn 
60 CONTINUE 

GXOL2aSORT(r.xa2) 

EPNUMao.n 

EINENUHan.o 
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OP TO KKmI, 101 

E;*N(tM>(GRAOXO(KK)-GRA03(KK)}ii*2‘,0*e2NiiM 
XMA63BA0»(CRADXn(KK)-GnAO3()(i:}) 

ETNFNUHaAHAXt (XMA63«ElNrNUtO 
70 cpnt:«ue 

E2NUHaSGRT(P2NUM) 

e*>RL2aE2NUM/GX0l.2 

CRRLIMfaEtNfNuM/GXOLINP 

MRITC(6*AA91 

AfcO FPRmaK///, (M ,a6X,37MPEI.ATIVe filPrEReNCpS BPTupCN REOflREP) 
«"ITC;i,A70l 

a70 FPRHATIim ,93X,23« *nO OBT^IliPp CRaPIE^Ts) 
mRITE(6,AA6) 

A6« FPRHAT(//ih ,98X,9HL»a ERROR) 

SO MRITE(6*30)ERRL2 
30 rnRHATUH ,aSX,2(9X,PtO,5)) 

RETURN 

END 

ccrcccccccccccrccpccccccccrccccccccrcrcccccccccrccrccfcccccccccccccccccc 
cc cc 

cc Tnis Subroutine estim4Tes lateral surface tehpepature cc 

cc 0T USE OF CUBIC spline IF THE USER SUPPLIES A OtSCRFTE S£T oF CC 

cc lateral surface temperatures, cc 


CC 


cc 


ccrcccccccccccrccrccccccrccccccccccccrccccccccccccrccrccc:cccccccccccccc 
subroutine SPLINE(XINT,Y1NT) 
cnMMON/C32/xo(ioo),Torioo),cj(fl,ioo),M 
dimension Cffl.iofl) 

EouivALENCErc) n»n*c(i*n) 

IF(XINT»X0{i)) 2» lf2 

1 YINTaTOd) 
return 

2 Kal 

3 IF(XlHT-xO(KAmt.,4»5 
« YiNTaYOtKAl) 

return 


5 KaKAt 

IF((M-k),OT'. 0) G0T03 
IF((M.K),LE',0) «aM-l 

6 YTNTa(XO(KAl )-XINT)A(Cn,K)«{xn'fKAn-XlNT)**2AC/3»Kl ) 
YINTaYINT*(XINT-X0C<))»(C(2»‘:)*>XlNT-X0<K))*«2»cC«.K)) 

RETURN 

END 

subroutine cofcen 

CCrCCCCCCCCCCCCCCrCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCrCCCCCCCCCCCCCCCCCC 


cc 



cc 

cc 

FINf) The SPt.INE CURVE FIT COCFFICIEnTS, FOR 

USB IN conjunction 

CC 

cc 

MTTh subroutine spline. 


cc 

cc 

INPUTS • 


cc 

cc 

M a NO, Of OaTa pairs 


cc 

cc 

XD a array op X (aRCISSA) VALUES 


cc 

cc 

YD a array of T (ORDINATBS) v*|UE3 


cc 

cc 

OUTPUTS . 


cc 

cc 

C a 2«n array of spline fit coefficients 

(a CnEFFiCICNTS 

cc 

cc 

PER triplet of data POINTS), 


cc 

cc 



cc 


ccrccccrcccccccccccccrcccccccccccccccccccccccccccccccccccccccccccccccccc 
cnMMON/C32/xo(ioo)»Tonoo),ci(«,ioo),H 
DIMENSION CM, 100) 

DIMENSION Pn0g),E(>0O),Aao4.3T,0n0O)'Z(lon),n(100) 
eauivALeNceTctn,t)ACU,n) 
c 

NpaM 

MaM«t 

no 2 Kai,M 
n(K)aXOCKAi)«xncK) 


Figure C-2. Computer Code List for Problems 
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o 


P(K)»OfK)/h, 

r (K)»<vn(K*i)-roo<) )/0fK) 

no 3 K«3*M I 

B(K)aE{K)-F(K-l) 

*Cl»2)»-l,«0n)/0(2) ( 

*(l»3)»naWD(2) 

* (2i 3) » rp (?)-Pt I )•*( 1.3) )/*(?,?) 

PC2)«B(2)/*(2#2) 

no <t k*3,m 

A(K.2)»2.«fP(K-l)*<’('<) )«P(k-1).a(k-1,3) 

P(K)«B(K)-P{K-n»6(K-l) 

*(K,3)»P(K)/*(K,2) 

B(K)«B(K)/i(K,2) 

OanjM-D/OrM) 

*(MO. l)al.*Q»AfM-l»3) 

A(N0.2)>-0«A(Nn, n*A(F4.3) 

P(NO)s0fM-i )»A( n0, 1 )*D(M) 

7(NO)a0fNDWA(NO.2) 

no 6 I>l.Nn«2 

KaNO'I 

Z(K)«B(K)-A(K,3)«7(KA1) 

Zn)«-An,?)*Z(2)*A(l.3)*2(3) 

no 7 K«1 ,M 

Q«1 ,/(6.*'OfK) ) 

r Cl .K)«Z(Ki»o 

n(?.K)«7CK*n<''J 

C(3.K)«vOO<)yO(M-Z(K)*PtK) 

n(O»K)«y0CKAl)/0(A5-ZCK»l)»PfK' 

return 


ORIGINAL PAGE l« 
OF POOR QUALITY 


END CQPGEN 
END 


figure C-2. Computer Code List for Problems 


Pl-1 and Pl-2 (Cont) 




C.4 


COMPUTER CODE LIST FOR PROBLEM PI -3 


The comptuer code for Problem Pl-3 is listed in Figure C-3. Before 
using this code, the user should review the remarks made at the end of Appendix 
A. 4. 



origiijal [S 

OF POOR Q JALITY 


*-• wtcu.i.i.'.wuwcw^vi 

cc 

PRncRAH PURPOSE* 

COMPUTf Thc'mflT ZO^C 3URF4CC cnNT«Ow FUNfTinN pCOUIREC fop 
SOLIO-^ELT interfaces aS OFXCRTBEO IN PROOlJH PJ-: ANC 
SOLVED IN C‘*APTFR « OF THE fjnai REPORT (To N*si) Entitled the 
CONTROL OF rLO*T ZONE INTEOFaPEF BY ThE USE OF sCLECTEC 
BOUNOARY conditions- BY SCIENCE APPLICaT iOnS, InC, 

SOURCE- 

SCIENCE APPI ICATIONS, INC, 

MUNTSVILI.E, ALABAMA 
AUT jRS- 

LA' •• M. foster 
JOHN MCINTOSH 

RFFERENCE- « 

. THE CONTROL CF float ZONE InTfRFaCES By THE UsE OF SPLEC . FO 

boundary conditions - 
(FINAL RFPORT - SAI-83/503«»HU) 

SCIENCE APPI ICaTIOnS 
NEMaRkS- 

- software developed and TES'Fn ON CCC TaOo/aROo 

-"^ALL^EOi'aTtOnS referenced IN CODE BElQN *«E COnTaInED in "mE 
FINAL report- 

- TwE control of float fPKE interfaces f«Y the LSe PF 
SELFCTED boundary COI.OItIONS 
INPUT variables and FUnCTIONS- 


cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 


p 

MSUH 

BASIS 


XA 

XN 

NP.RID 


nR 

RKS 

RKL 

RL 

SI ENCTh 
Q 

IHFC 


(X0,TD) 


PECLET NUMflC" 

- number of terms tN series expansion of 

TpMPCRATUHE CISTRlpUTlON (ThE DFSIrED SOLUTIOI.') 
-USER PROVIDED FunoTIOnS tO EXPanq THE MELT ZONE 
surface control function. SjE EO, (4.0. IBJ OF 
FINAL REPORT ANO SuBROUTINt BASIS of THIS CODE. 

. axIaL position Of L0“€R End OF CvLlNoER 

- axial position Of upper End of Cylinder 
. number of grid points USfO IN numerical 

SOLUTION OF 0. 0. F. BOUnOART VALUF PROBLEM 

resulting from Tr*m3F0Rm»T10M of PnE The modeling 
temperature 

- number divisions of CYlIHoFR BaOtUS msec in 
OUTPUT OF temperature OISTRiBUTTON 

- solid Thermal conductivity 

- Lifluio Thermal conouctivitt 

- product of cRVSTaL cROhtm Rate, sCLId dcnsI'y, 
and latent meat of fusion 

. length of solid PfClONS To BE COkSIOFREP. 

• length of helT zone 

- I IF A discrete data point FORM oF ThE SURFACE 
temperature is use* provided 

0 IF A user OEFlNCf. functional FORm OF THE 
surface temperature is provided 

• USER provided data PTS FOR THE ayIAL DISTANCE 
(XD) and CORRESPONoINC surface TEmpCRaTURC (TO) 

. number of OaTa RtS, input if THfr ■ 1 


Wb 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 
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cc 


SET TO 0 IF IllFC s 0 

CC 

cc 

MFC 

. USER provided (IF IHFC « 0) SURpiCE TCMrCRA^ijRC 

cc 

cc 


FUNCTION 

CC 

cc 

CASE LI“!TS 

. This program genfRates the suoPAfc contpol 

cc 

cc 


functions for VARinUS COMBINATIONS OF THE IllOP* 

cc 

cc 


LIMITS mTERM anO M5YS (SEF F». (•.".18) - (0.0. 2J) 

cc 

cc 


OF final REPORT). TO oFFINE ThESE fOMBTNATIONS 

cc 

cc 


mInTERM - tmC mIntmUn aLI.ONfO value 

cc 

cc 


Cf nTERm 

cc 

cc 


haxTERm - TmE haxtnUn alI.OUfO valuP 

cc 

cc 


Of nTERm 

cc 

cc 


OFLTERM - Increment of mTEOm from riUTERM 

cc 

cc 


To mAXTERm 

cc 

cc 


mInNSTS - TmC mTntmUh aLLOUfO vaLUC OF nSTS 

cc 

cc 


maxnsts - The maximum aLI.Ol'fO v»luP of nSts 

cc 

cc 


DELTERm . IncRe"EnT of NSTS PROM MiNfjsvs 

cc 

cc 


TO maxnsts 

cc 

cc 

OUTPUT VaRIaBi 

LES- 

cc 

cc 

THOLO 

. TEMPERATURE OISToiBUTION a»RAT pnR EaCM REGION 

cc 

cc 


- SEE Subroutine mflt of cOof 

cc 

cc 

GRAOXO 

- AXIaL THERMaL GRjOIENT aT Xq FOR RECTOnC. 

cc 

cc 


(XO IS SET TO 0 FOo UPPER SpLID REPION, aKD TP 

cc 

cc 


negative SLEUCTh for LOHfR SOLID RpCION) 

cc 

cc 

GRAOXN 

- axial thermal gradient «t xn for regtonc. 

cc 

cc 


(XN IS SET TO SLENrTH ♦ t) FqR UPPER SOLID REGION, 

cc 

cc 


ANO TO 0 FOB LOUeR SOLID REgIQn) 

cc 

cc 

GRAOATq 

. axial Thermal gradient at bottom of melt zonf 

cc 

cc 


(SEE main of COOF) 

cc 

cc 

cparato 

. AXIaL thermal gradient aT TQP of melt ZPNE 

cc 

cc 


(SFE main of CODE) 

cc 

cc 

CPOLT 

- aRBaT of COCFFFCvCNTS used TO EXPiNO The MELT 

cc 

cc 


surface CONTROL FUNCTION (SeE EO, f«.0'.18J OF 

cc 

cc 


final report *ng Subroutines mci.ti and hfc of 

cc 

cc 


cone) 

cc 

cc 

ERRL20 

- THE Relative L2 "IffeRpnce betmepu the 

cc 

cc 


OESIRFO gradient a^ X«X0 aNq THE CpAOIENT 

cc 

cc 


obtained BT using the surface CONTROL FUnCTIOI.’, 

cc 

cc 

ERRL20 

- The Relative l? oIffcrence betmefu the 

cc 

cc 


desired gradient av X«»N aNi) the CB4DIFNT 

cc 

cc 


obtained aT USII.G VHE 3URF*ce CONTROL funCTIOI:. 

cc 

cc 

user supplied 

haThmaTIcaL SOFTliARF- 

cc 

cc 

- A least squares aLC!)R!TMH to r!T A fuMcTion Tn A iINfAR 

cc 

cc 

COMnINATtON 

OF SeLFCTEO FUMCTIOnS 

cc 

cc 

- AN algorithm to evaluate bessfl functions (Reouiped in 

cc 

cc 

subroutine .10) 

cc 

cc 

- A numerical integration routine (REOUIrFO in sUBRnuTiNES 

cc 

cc 

XNTegU and 

INTEGL2) 

cc 

cc 



cc 


ccrccccrcccccccccrcccccccccccccccccccrccccccccccccrccrcccccccccccccccccc 

PRnSRAM MATN{INPUT, OUTPUT, TARE^»lNPUT,T*PE6«0tJTP'JT) 


CnMHON/BE*0’/P**^TERP,MSUM,X0,)tN,:;CBI0,NP 

crtHHON/cai /THfT*B (2 «. 505) » thclom 0 n * t( ) c ) ,6B*o»n (ion, craoxo ( mi) 

CnMMON/C22/TCASF,THeuT(3) 

cnHH0N/C9</RKS,RKL,RL,NSYS 

CnM**ON/CP5/KR*0*T0(lOl),6RAOATO>J0n.RH3) (20) ,Rii52f ?0 ) , 3 (20 ) , 0# 
1*P*T1(20, J0),*«*T2(?0. lO),AL2fa«,2o),PH3(«*),t»Mif (moO),II**K(20) 
CnHHON/COO/RLENCTH 

coMHOR/c32/xo(mo),To(ioo),cira, ioo),** 

OIMeNSlOH.RflOl) 

ccrccrcccccccccccrccccccccrccccccccccrcccccccccrccrccfcccccccccccccccccc 

cc 

CC OFTERMINf TFHPERaTURC OISTRIOUTtON aNO «R*OIENTs FOR UPPER CC 

CC SOLID RERION rr 

« CC 

ccrccccrcrccccrccrcccccccccccccccccccpccccccccccccrccfcccccccccccccccccc 


Figure C-3. Computer Code List For 
Problem Pl-3 (Cont) 
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ORIGINAL PAGE 19 
OF POOR QUALITY 




• fT.'. 


tj 


r • 


ORIGINAL PAG!: E9 
OF POOR QUALITY 


ICASE «3 

10 FO»s*T( jMi , ja*.??*"' PPER SOLfP) 

C*UL INPliT 
XK» 0 *SLCNGTW 
xnsQ 

c*u. ielt 

Lcrccccrccccccrccrcccccccccccccccccrcrcccrccccccccrccrcccccccccccccccccc 
cc cc 

cc compute GRAnlFNT IN mCLT ZCNC a» UP»E" InTePF*CF CSeC COUA"iOn cC 

cc FTC. 1 - 2 .' AND STOBE RESULT IN RPAOaTO CC 

CC CC 

ccrccccrccccccrccrcccccccccccccccccccrcccccccccfccrccrcccccccccccccccccc 

00 20 isi.im 

C»ArATQ(l)sfRKS«C»*OxO(n*nU/RifL 
20 continue 

ccrccr.crcrccccrccrcccccccccccccccccccrccccccccccccrccrcccccccccccccccccc 
cc cc 

cc determine temperature orsTniriuTTpN ano <*p*urENT tor loker cc 

cc SOLID RECIC*: CC 

CC cc 

ccrccccrccccccrccrcccccccccccccccccrcrccccccccrcccrccrcccrcccccccccccccc 
ICASE =2 
mRITEio.^o) 

JO FOBMAT(IHl,aax,22NL Ok^EC «OltC) 

CALL INPUT 
XNSO 

xPs-SLENGTh 
CALL melt 

ccrccrcrccccccrccrcccrcccccccccccccccrccccccccccccrccrcccrcccccccccccccc 
CC cc 

■cc COMPUTE gradient in melt ZCNC at LOtER InTERFacf ( 3 fE COUA-lON cC 

CC FZa, nc, l-rj) »N 0 STORE result in GRaOATO cc 

CC . . cc 

Ccrccrcrccccccrccrcccrcccccccccccccccrccccccccccccrccrcccccccccccccccccc 
Of ao Ia»,;ot 

CRAOATotI)sfR»(S»6RAOXN(I)*RL3/RMt. 

UO continue 
c 
c 

NRITEfO.aoOT 

nOO FORMATClMli xax.SpMM ELT Zdnf THERMAL GRAOI 
« E N T S) 

rtOl EORmaT(/,ih ,50X,2aHA T IntE»eACE) 

NRITEC6,«0n 

■PITEf 6 . 72 ) 

T2 EriRMAT(///,aaX,lHR,7X,15HCnA0. *T X» 0 . 0 , I 2 X , I JmCR AO . AT X« 0) 
NRITE{6«a02T 
n02 EORmaT(/,ih ) 
on pa? KKai.lCt 
RfKtc;»(KK«n«o.'>i 

••RITEC6. JOOiR(ER).G'’AOXO(KKJ,GR»OXN(KO 
XOO eORMatC^H , 39 X,E«.fc, 3 X.ei 7 .'?, 7 »,El 7 .p) 

oaT continue 

ccrccccrccccccrccrcccccccccccccccccccrccccccccccccrccrcccccccccccccccccc 
cc cc 

cc determine melt zone control ruNpTION (EOuatION a. 0 . 18 ) ano cc 

cc resulting TfMPERaTURE DISTniRUTTCN cc 

cc cc 

ccrccccrccccccrccpcccrcccccccccccccrcrccccccccccccrccpcccrcccccccccccccc 

ICASE-Il 

NRITEC6.R83) 

o83 EORMATtlMlfROX, 18 mm E L T Z n N E) 

call input 

CALL MELT! 

STOP 
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CNO 

ccrccccccccccccccrcccccccccccccccccccfcccccccccrccrccrcccccccccccccccccc 


cc cc 
CC PIIRPOSC cc 
cc - 0CTCR“TNC zn*ic TEHPr»ATilPC OISTPIRUTION cc 
cc • COMPuTF Rft*TTVC FPROPS neTwffK gRaoIInTs REOPlREn THC CC 
cc HELT ZriNE InTFREaCES AND TliOSr OBTaImEo USIpC “HE sucfacc cc 
CC C^lNTROl. rilHfTTP** cc 
cc cc 


ccrccccrccccccrccrcccrcccccccccccccccrccccccccccccrccrcccecccccccccccccc 

StiBROUTTNE foster 

CPHMON/RFACt /P<xTFPH.MSUM,XO,XN.t: 6 RtO.HR 

cnHHON/C?l/THETAB{ 20 ,SOS),THCLnfion.T(‘ 0 ), 6 RAO»iUlM),CRAC«Otini) 
CnHHON/C 22 /IC ase , TmElT ( J ) 

CnMMQN/C2C/RF8,RKL."t*NSV3 

CnMHON/C2S/RRAOATO(lon»6»AOATnfioi),»MSi{8n),fJi;G2(?0),S(20)»0, 
lAHATt C20» lOTf AHAT2(20« 10]«AL2(4A,2O)>RHS(a4)>HUH(lSOO)« IIH|t(20) 
CnHMQN/CaO/SLENfiTH ■ 
xnao 


C 


XNsQ 

CALL melt 


50 


60 


670 

666 

667 

660 


GATQL2=0,0 

CaTQLINsA.O 

GATnL2s0,0 

GaToLINsa.o 

OP so j«» , mi 

GATnL2«6ATQ| a*GRAOATa(J)*Gr;ACATn(J) 
GATnL2aGAT0l 2 *GPAOATO C J ) aGnAGATn ( I) 
XHAGlSABStGRAOATOCjn 
XHAG2saBS(CRA0aTQ(J)} 

gatolinsamaxi(xmagi#gatolii:i 

GaTQLINsAMAXI (X>«AG2f GATOLIM 

continue 

CAT4L2»SORTrGATOL2) 

CATQL?aS0flTrGAT0C2) 

C?N(IMOs0,O 
E?NIIMQ*0,0 
ENFNOaO.O 
EHFNQao.n 
00 60 

E2NllMOaCr.RAnATO(K)«CR*OXO(K)3*«?.OrE2*IUMO 
EPNIIMOa (ORAP A70 ( K 1 -CBAOXN (K ) ) •*?. 04E2HUMQ 
XHAg3>ABS(GRADATO(K}*6RAOXO(K)) 
XMAG4aA8S(CRAOATO(K)-6RAOXr.(K)) 
ENFNOaAMAXl fXMAG3.ENFN0J 


ORIGINAL PAGE 13 
OF POOR QUALITY 


ENFHOaAMAXj f XHAG4,ENFNQ) 

continue 

C?NUH0aSQRTrE2NUM0) 

C2NIIMOaSORT rE2NUN0) 

ERLINFOaFNFNO/GATOLiN 

ERLINFQbFNFnQ/GATQLIN 

CRRL20aE7NUM0/6AT0L2 

ERRL20aE2NUM0/GATQL2 

h»ITE(6.660T 

FORMATf///, IH .46X.37HRELATIVF pIFFE»FnCES SFTufCN BEOUIREP) 
hRITE(6»670) 

FORhATOh >537. 23h anO OBTaILEP CRaOIFnTs) 

HRITE(6.666) 

FORnaT(//,im .5RX.OHL-2 ERROR) 
h»ITE(6.667) ERRC20 
FPRmaT(IH .42X.4MAT 0.2(5*.FI(j,f;)) 
mRITE(6.66B)EPRL20 

FORMATOM .42X.4HAT 0.2{9X,FI0.5J) 

return 

ENO 
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siiaf-ouTiNP melt 
«E*L JJ ,.tiLA*‘,K*'BSjO 
COMMON /Cl/Pl A“0(?8) .Jl (20) , JiLAM(20) 
cntiMQH/ci 0/iSCRTPf2O),aSCPIP(2n' 
cnM«ON/priO*/p.MTepM,**suM,xo,»N,i;CRio,N» 
CnMMON/cS/Rnon,PSI(20.10t )»SO.ri(2n) 

CnrtMON/C20/*(50 0).H(5<'0),CC5«tfl),C(5'’0),V{S0O1,BFTACS0S),C*HMA(5n5) 
cnMMON/CPl/TMFT*B(20,S05),TMCLr)rioi),Ttlo),GR*D»'.(lOl),(!n*PxO(inn 
CnMMON/C?2/TCASF,TMELT(3) 
cnM»,0N/C?3/f:fli0?(101 ) .CR*03( lot < 

CHAR4CTFR*t7 RIS,*LPH*(6) 

CMARACTER*iR STARS. STARC6) 

OaT* flIS/'R« '/.stars/'*****************'/ 

on 207 LSI,/. 
aLPha(L)*RI.S 
STaR{L)sSTaRS 
207 CONTInuC 

ccrccrcrccrrccrccrcccrcccLcccccccccf'-rcccccccccCccrccrcccccccccccccccccc 

cc CC 

CC AMOCM)-.RCnT OF JO BESSEL FCn Cr 

CC J1{M)=JURLaMD(")) "MERE J) :s rESSEL FCn cl 

cc J1LAH(M)=JI fM)/RlAMD(M) Cj; 

CC 

ccrccrcrcccccccccrccccccccccccccccf.ccrCcccccccccccrccrcccccccccccccccccc 

RLAMOt n*2' ao«R255577 
RLAmO( 2>sS.520O78t 103 
RlAMOt 3lsfl‘.653727912R 
RLAMOC oi»ii ,7Rt53a«3Rl 
RIAMOC Slsla.R30R177oA6 
RIAmOC 6)alA.071 063**67R 
RLAMOt 7l«2t .2H63662RR 
"I AMOt 8is2a. 3520715308 
RlAMO( Rls27,aR3a7*»l320 
31 AM0(10)s3ft,63O60hO6Aa 
RLAmoCI I)s33,7758202l36 
Rl AMC(l21a3A.RlT0R«35T7 
RLAMO Cl 31 *40, 0580257606 
RLAM0(iaia43,lRR7Rl7l32 
RLAM0(lS)*46.301 1883717 
Rl AM0(I61 s 4O, 0826098470 
Rl AH0Cl7)«52. 6200518011 
31 AmD( 181*55, 7655107550 
R| AMO{l9l*58,906983R26l 
RL AMO (20 1*62,008069 1902 
J){ 11*0,5191470973 
J)( 2)3.0.3402648065 
JtC 3)30,2714522999 
jlC a)*-o, 2820598310 
Jl( 3)10,2065460331 
JlC 6)»*o. 1877288030 
JlC 71*0,1732658902 
Jl ( 8)3-0.1617015507 
JlC R1=0, 1521812138 
Jl ClO)s-o. 1401659777 
Jl (111*0,1372969030 
Jl (1213-0.1313206267 
Jl (131*0.126069097! 

Jl (14)*-n. 1213986208 
Jl (151.30,1 1721 11989 
Jl (16)3-0. 11342919E6 
Jl (171*0,1099911030 
Jl (l8)*-0. 1068078883 
Jl (19)*0, 1039595729 
Jl (20)*-n, 1012930989 
on 566 r*l,20 

J1LAM(I1*.M CIl/RLAMOm 
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snji (Dcji H) 
continue 

ccrccccrccccccrccrcccrcccccccccccccccrccccccccccccrccrcccrcccccccccccccc 
r.c cc 

CC FTNO COEFS rc» HfSSEL E*F*NSI0M« of AmO cc 

cc SFE equations (?.2.1T) and (2.2*.i«) OF EINaL »EpOBT cc 

cc . . cc 

ccrccrcrccccccrccrcccrcccccccccccccccrcccccrccccccrccrcccrcccccccccccccc 
CALL ArCn.n.AGFl) 

CALL BFCM.n.BCFi) 

00 20 

B(I)«(T-n«o,ni 

RHQLO«R(n ^ 

CALL AFCCRHOLQ.anS) nOlfilNAL PAGE t3 

A(n»ANS-AOFl rtliai ITY 

CALL BFC(»NOLn,ANS) OF POOR QUALIIt 

BCDsANS^nOFl 

20 continue 

C*uL COEFS(R,A,tO!.20.A3CRIP) 

CALL COEFStR, 0.101, 20, BSCRIP) 

ccrcccccccccccrccrcccccccccccccccccccrccccccccrcccrccrcccccccccccccccccc 
cc cc 

cc SOLVE FOR ThETa 0aR OF EQUATIONS (2.2.t'») 0T SOiVInr THE CC 

cc TRiniACONAL SVSTep (2.2.20) - SfE FINAL REPORT CC 

cc . CC 

cc.-ccccrccccccrccrccccccccccccccccccccccccccccccccrccrcccccccccccccccccc 
oxs(XN.vn)/NGRin 
0X2>0X*0X 


■ L«NRRID-» 

00 S6S SatfMSUH 
DO 00 I«l.L 

ACI)al.n»0X*P/2.0 

8(n«-2’.0-OX2*RLAMO(Mj«RLAMf'{M) 

ctllat .0-OX«P/2.0 
X»X04I«OX 

CALL GBaR(h,X,ans) 

0(na0XR«ANS 
00 CONTINUE 

0(l)»0n)-(l .OsOX*P/2,0)«*CCPTp(N)«3QJl (M)«0,5 
0(L)aOa)-( J .'’-0**‘’/2.0)«BSCRTr(H)«S0Jj (mJaO.s 

CALL TRlnAC(L) 
on so ia2,NCBI0 

IlaI-1 

THETAB(H»naV(II) 

50 CONTINUE 

NSTOPaNGRrOsl 

THETABfH,! )jaSCRIP(P)*3QJ1 (*‘)>2.0 

TMETABfM,NSTOP)»B8CRlP(M)«50j; C”)'2.0 

S65 continue 

ORal ,0/NR 

NR3T0PaNR*l 

on hO lat.NRSTQP 

R(i)«(i-n«oR 
on 65 Hal.HSUN 
psitM,naF(M,Brn 


o5 CONTINUE 

oO CONTINUE 

ccrcccccccccccrccrcccccccccccccccccccrccccccccccccrccrcccccccccccccccccc 


cc 


CC 

cc 

PRINT TEMPERATUREB 

cc 

cc 

. . • 

cc 


CCrCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCrCCCCCCCCCCCCCCCfCCCCCCCCCCCCCCCCCC 
IFCICASE.CQ.3) C0T06TB 
»RITE(6.30) 


30 F0RmaT(IH1,52X,22HL ONER sOliD) 
•' G0T067R 
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^7e cpnti'iue 
"BITE f6, 10) 

10 FOSM*TnHl,«;2X,22HU PPER ^OUlO) 

hV> continue 

«BrTE(6.T0) 

70 FOPmaTCIM ,u5X,a«»M TEhpfP/iTupE niSTPiauTlO 
In) 


IFLAGaO 
HPir.HT*6 
Ml EFT = 1 

I9n CniJTiMUE 

IF(hRSTOP.i,F.mPTGmT)TFI.AGsI 

HOiC.HTsMTNOfNRSTOFfMBIGHT) 

»»ITEt6» '90UP( I).J»»*l.EFT,f;P:CHT3 
190 FORMATC/////, iH , l7X,fc(F12.8,SX' 3 
«»ITE (6,2671 ( ALPHA (U,L*t,'*ritGHT) 

?67 F0RHAT( IH*, 1 7X.6AI 7) 

-BITE (6, 2681 CSTAR(L) ,Lst,tf»lGH-r) 

?68 FORMATCIHO, i5X,^A)7) 
on 200 laljNSTOP 
ISKIP=I-1 

IHOLD=(I5iK:P*0.10000 0'n)/13.0 
XhOLC=CISKIP/10,0)-IhOLO 
IF(XHOLO.GT'.O.On5) GOTO200 
IIaNSTOP>i-I 
X»XO+(II-l)«OX 

DO 202 JsmlEFT, “RIGHT 

ccrccrccccccccrccrcccrcccccccccccccccrcccccrccrcccrccrcccccccccccccccccc 


cc cc 
CC DFTFBhINF TrHPEBATUBE AT (X.ni.I') cc 
cc SEE equation (2.2.1“) OF FINAL oE?CBT CC 
CC CC 


ccrccrcrccccccrccrcccccccccccccccccccrccccccccccccrccccccccccccccccccccc 
THCLn(jiso.o 
DO 2<la Mel, “SUM 

THntD(j)sTHnLO( J)»2.Q*FsrrH. J)*TwtTAB(M, lT)/SOJl (“) 

20/1 CONTtnup 

CALL hfc(X,anS) 

THQLn(JlsTHOLD(J)+ANS 
20? COnTtnuf 

mRitf (6,210 IX, (ThOLD(J3,J«MI CFT.moiGhT) 

?10 FORmaT(5H X»,Fi0.6.5h • ,6(EIR.8,2X)) 

20 n continue 

IF(IFLAC.FQ.1)G0 to 220 
MRIGhTbMRtGHT*6 
HLEFT»MLEFT+6 
GO TQ 180 

22« continue 

ccrccrcfcccccccccrcccccccccccccccccccfcccccccccccccccrcccrcccccccccccccc 
cc cc 

cc compute thermal G»a 0 iEnt 3 at X3»0 anO xn CC 

cc cc 

ccrccccrccccccrccrcccccccccccccccccccrccccccccccccrccrcccccccccccccccccc 

39«: continue 

IFdCAse.EQ'.I) GCT0i8t 
HRITEC6.10) 

G0T0682 

n81 continue 

hBITE(6. 10) 

fc82 continue 

-BITE(6,7I) 

71 FORMATdH ,a7X,35H T H E R h a ( 

«RITE(6, 723X0, XN 

72 FORMAT(///,aax, 1HR,5 x, llMGRAn, aT 

2 ,//) 

DO 230 I«l,i01 


CRAOTefiTS) 

X»,F10'.5, 14H GRaO, AT X»,Flo.5 


ORIGINAL ^ 

OF POOR QUAUTY 


Figure C-3. Computer Code List For 
Problem Pl-3 (Cent) 
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Rf i)«ci-n*n.ei 

00 ?.H0 Msi^mSuh 

C*LL jn(V*R.Y5 PAGE IS 

PRHM.rj.v original PAVat w 

24A contimuc -c poor quality 

asrt continue ^ 

cn«-ox 

oo aso i«i , 1 01 

on 260 J3i ,s 

TCj)«o,n 
00 270 Ms),HSIlM 

T(.n«T(J)*2,0*RSI(H,t)«TNFT*efR»J>/S8J> (H) 

27ft CONTINUE 

x«xo*(j«i)*nx 

CALC HPfCX.ANS) 

T{JJ»T(.r)«AN3 
26ft continue 

on 2 ao jaA.to 

TtJ)«0.ft 

JHOLflsNSTCiP-jO + J 
00 aon m=i,>isUm 

T(.I)aT(J)+a,0»PSl(M,i)«THFTAB(H,JN0L0)/SCjl(N) 

29n CONTINUE 

x»xo4(.iHOi.o-n*ox 

CALL hfcCX,AN3) 

T(J)sT(J)*AN 3 
28n CONTINUE 

ccrccccrccccccrccrcccrcccccccccccccccccccccccccrccrccfcccccccccccccccccc 
cc cc 

CC approximate TMEPmaL gradients at CyLInOFR ENOS cc 

cc • SEE EOHaTtCNS (2.2.21) ANO (a'.2.22) Of FINAL rCPOrT CC 

CC . cc 

ccrccccrccccccrccrcccrc' ccccccccccccrrcccccccccccrrccrcccccccccccccccccc 
C»AnXMcI)«(-3*T(6)T«6*T(7)-36«Tfe)*a8.T(9)-25*Trl0))/(12*O0) 
GRAnX0CI)»(-J*T(5)Tl6»TC«)-5«i«Tr3)^«8*T(a)-25»Tf n)/(l2*0X) 
ISKIP«I«» 

IHOLO«nSKIP*0,ftOOOOoni)/10,Q 
XMOLOa(ISKip/'JO,0)-lHnuO 
IFCXfinLD.CT.O.OOS) GOT0250 
NRITEC6»300)Ra).GRAOX0(n,CRAOyN(p 
TOO FORMATcJN ,T9X,Fe.6.3X,£17,9,7X,E17,9) 

25ft continue 
RFTURN 
END 

ccrccccccrccccrccrcccccccccccccccccccrcccccccccccccccrcccccccccccccccccc 
cc . cc 

cc TNI3 SUBROUTINE APPROXIMATES (flv FINITE DIFFERENCE) G BAR OF CC 

cc EOUaTION ( 2 ’. 2 . 16 ) OF FINAL "CPORT CC 

CC .. . , CC 

ccrccccrccccccrccrccccccccccccccccLccrccccccccccccrccrcccccccccccccccccc 
subroutine G8aR(H>X.anS) 
real J1..ULAM 

COMMQn/CI/Rl AMO <20 )»J1 C20),J1lam(20) 

CnMHQN/RFAOi/P,MTERM»MSUM,XO,XN,HGflID,NR 

EPSLQN«o,01 

XJaX-EPSl.ON 

xaiXTEPSI.ON 

CALL HFC(X,AN3) 

CALL HFC(Xl,AN3t) 

CALL HFC(X2.ANS2) 

G»P«(ANSa«AN3l)/(2,0«EPSLOt;) 

G«C-(AnS2tAn3j»2,0*an3)/CEP3LON«E'’3LOn) 

AN3«C*Jll.AMfM) 

return 

ENO 
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ucrccccrccccccrccrcccccccccccccccccccrcccccccccccrrccccctccccccccccccccc 
cc cc 

cc PURPOSE . cc 

cc SUPPLY input o*ta. See *pPE>*r,rx a. a rnn details cc 

cc cc 

ccrccccrccccccrccrcccrcccccccccccccccrcccccrcccccrrccfcccccccccccccccccc 

SUBROUTINE ynPUT 
INTrCER OELTERm.OELNSYS 
CnMPON/PfAOl/P,»‘TE»*',NSUM,XO,JrN,NGPID,N» 

CnMMON/CP2/TCASE.THELT(3) 

COMMON/Cpa/RKS.RKL/RLf N3Y3 

CnMMON/C25/r,flAO*TOnon ,CRhDATCMOl )»RMSi ( 20 1. nwC2 f?0), 3 (20 )/Q» 

t AH ATI (24, 101,AMAT2(2O.1 O}«aL2(4a.2O)>RHS(A4).h»;;(1SoO],IZN(;( 2O) 
COMM0N/C31/THPC 
cPMMON/C32/xonoo)»''onoo),ci(a, ioo),H 
COMMON/CaO/SLENGTH 

CDMMON/CSl /MAX TERM INTERM, HAxfji»YS,-JlNN3y 3, DFLTpnH,nELf»'SY3 
NRITEf6,S) 

5 FnRMAT(//,lH ,Snx, 2 ftMl NPUT DATA) 

GO TO(1O,20.2O),ICaSe 

ccrccrccccccccrccrcccrcccccccccccccccrccccccccccccrccrcccccccccccccccccc 
cc cc 

cc mflt Parameters cc 

cc .. cc 

ccrcccccccccccrccrcccccccccccccccccccrccccccccccccrccccccccccccccccccccc 

10 READ(S, 12)P,MSUh,NGRI0,NR 

12 FnRHAT(E20.(0,atl0) 

NRITE(6.25) 

25 FORMAThm ,a8H P ncRjd nr "SUM) 

HRITE(6#ta) P, ncRI0,nr,m3Uh 
14 FORMATCF^O. (0,3114) 

REAO(S, 1A)MAXTERm.hINTCRM,MAXN8vs,mINN3TS,DE1.TEbI;, 0ELNCY8 
16 FORMATCinrS) 
hRITe(6.T9R1 

T99 FnPMAT(//,IH ,9X,57HMAXTERM HTf|TERM OElTF»H h*XNSTS MINNSY 
•S OELNSYSI 

NRITE(6, 18)haXTERm,hinTERM,0CLTfRm,haxnSyS,hinnsYS, 0ELN3)'S 
la PHRMaTHh ,«!X,6(5X,I5)) 

RFTtlRN 

ccrccccrccccccrccrcccrcccrcccccccccccrccccccccccccrccfcccccccccccccccccc 
cc cc 

cc loner SOl.rO/UPPER SOLID PaRAMFTfRS cc 

cc . cc 

ccrccccrccccccrccrcccccccccccccccccccrccccccccccccrccrcccccccccccccccccc 
20 continue 

RFAn(5, 12)P,M3 UM,nCRID,NR 
NRITEC6.25) 

h»ITE( 6, ia)P,NCRlOiNR,MSUM 
RFAn(5.22)RKS,RKL»RC,3LeNGTM 
22 FORMATcarao" la) 

IFCICASE.NcCs) SQTQ26 

RE*0C5,22)Q 

NRlrE(6,aaa) 

P88 format {//, iH , 14X,3MRKS,17X,3 hRi,L, 1TX,2HRL ,20X,7’H3LENGTM,8X, 1 IMMEL 

-T LENGTH) 

NRITE(6,20)RK3»RKL/PL«SLENCTi,,Q 

24 FORMAT(5E2o'.io 5 
2b CONTINUE 

IFCICASE.EQ'.S) G0T048 
hRITE(6#«88) 

nRITEC6,20)rk3,RfS,Rl»SLENCTN,(j 

a8 CONTINUE 

ccrccccccccccccccccccccccccccccccccccrccccccccccccrccrcccccccccccccccccc 
cc cc 

cc SPLINE INPUT option cc 

cc cc 
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tcrcTcrccccccrccrcccrccccccccccccuccrccccccc'cccrrccrcccccccccccccccccc 

Br*n(5, 

'i»ITEC6..'fO}M 

JO FnRMAT(//////,9SH TMe SUPfUCC TFl'RfRATUPt 0 ISTffTnUT TON IS APPPOXIM 
2*TEn BY the cubic spline THPnur.M the FO.LOwiNG,Ta,2jH (X, TENPl 0 
• *T* POINTS,///, J7H X SURFirE TehP',) 

IFCTHPC.EQ.n) RETURN 
00 32 I»1,M 
RFAn(S,22TXn(n.YDtn 
NRIT£(6»Ja)xo(I’.Tn(i) 

3« FPRNaT(2E20'.1O) 

32 continue 
CALL COFr.EN 
RETURN 
CNO 

ccrcccccccccccrccrcccrcccccccccccccccrccccccccccccrccrcccrcccccccccccccc 
cc cc 

cc This subroutine prqvioes for oata input cc 

cc this subroutine supplies LATfRAI SURFaCE TEHPERirURE CC 

CC • SEE EOIIaTtON (2, 2, a) OF FINAL REPORT CC 

CC ... ..... CC 

ccrcccccccccccrccrcccccccccccccccccccrccccccccccccrccrcccccccccccccccccc 

SUBROUTINE HFC(X.*NS) 

COHH0N/C2R/BKS,Nt<L.«L.NSYS 

C0MM0N/C3t/THFC 

COMMQN/C26/rPOLT(20) 

COMH0N/C22/TCASE,TmElT(3) 

CO TO(10,20.3O]«ICaSE 

ccrccccrccccccrccrcccccccccccccccccecrccccccccccccrccrcccrcccccccccccrcc 
cc cc 

cc melt SUPFaCF CONTROL TEMPER aTijRc, SEE EOuation S.O.lS CC 

CC • .. . CC 

ccrcccccccccccrccrcccccccccccccccccccrccccccccccccrccrcccccccccccccccccc 

to ANSaO.O 

00 12 K«1,NSYS 
call BAStsCk.X.anSI ) 

ANS«AN8TrP0l.Y(K)*ANSl 

t2 CONTINUE 

return 

ccrccccrcccccccccrcccrcccccccccccccccrcccccccccfccrccrcccccci:cc:cccccccc 
cc cc 

cc loner SOLID SURFACE TEMPERATURE COMPUTED NgxT CC 

“ .. . CC 

ccrccccrccccccrccrcccrcrccccccccccr.ccrccccccccccccr.ccrcccccccccccccccccc 

IFUMFC'fiQ.IJ COT022 

ccrccccrccccccrccrcccccccccccccccccccrccccccccccccrccrcccrcccccccccccccc 
cc cc 

cc USER SURPLIrO LOmER SOLID SURfafE TEMP, OISTRTBiiTIOn CC 

CC placed mere cc 

, CC 

ccrccccrccccccrcccccccccccccccccccccccccccccccccccrccccccccctccccccccccc 

PPTttPM 

22 CALL SPLINEIX.aNS) 
return 

ccrccccrccccccrccrcccccccccccccccccccrccccccccccccrccrcccccccccccccccccc 
cc cc 

cc UPPER SOLID SURFACE EmPERaTURF COMPUTED N£XT CC 

CC CC 

ccrcccccccccccrccrcccccccccccccccccccfcccccccccccccccrcccccccccccccccccc 

30 IFCIHFC.EC..,) COT022 

cc USER SUPPLIED UPPER SOLID SURFACE TEMP FlACED HeRC IE IN CC 

CC functional form CC 

RETURN 

End 

ccrccrcrccccccrccrcci ccccccccccccccccrccccccccccccrccrcccccccccccccccccc 
cc .... cc 
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cc 

cc 

cc 

cc 


this SUBROLTIHE supplies RAOIai. T£mPE(|*TuPE nlOTHIBUTirN 
ON lower ENn OF C'^UNOER • SEE fOUaTIUN f?.2,2j OF FINAL 
REPORT 


cc 

CC 

cc 

CC 


ccrccccrccccccrccrcccrcccccccccccccrcrccccccccccccrccrcccrcccccccccccrcc 
subroutine AFC(Rt*NS) 

COMMON/RpiO' /P*“TFRP» hSUh,XO,XN,NCBIO*nR ORiniMAI 

CALL HFCfXJ.ANSl 

return of poor 

ENO 

ccrccrcrccccccrccrcccccccccccccccccccrcccccccccccccccrcccrcccccccccccccc 
cc , . cc 

cc This subroutine supplies radial temperature nlSTRlBUTIfN cc 

cc ON upper end of cTLInDER - SEE equation (2.2.2) CF FINAL REPORT CC 

cc .. . cc 

ccrccccrcccccccccfcccccccccccccccccrcpcccccccccccccccrcccrcccccccccccccc 
subroutine pFC<R.*N8) 

COMHON/RFAOl /P.HTERM,NSUM,XO,*N,KGRtO,NR 

call HFCtXN.ANSl 
return 

end . . . . , 

ccrccccrccccccrccrcccccccccccccccccrcrccccccccccccr.ccrcccccccccccccccccc 
cc cc 

cc this subroutine fits BESSEL SFRtCS TO 0 *TA by LfAST squares cc 

cc method • SEE equations (2,2,17) * (2,2. Ifl) AND r2.2’,23) CC 

cc OF FINAL REPORT CC 

cc . . cc 

ccrccccrccccccrccrcccrcccccccccccccccrccccccccccccrccrcccccccccccccccccc 
subroutine rOEFS(fl,Y.NR,NCCFr,COEF) 

INTfGCR nr.ncoef 

real F.RMO) ).Y(101),COEFC20),WpC«60) 
external f 

ccrcccccccccccrccrcccrcccccccccccccccfccccccccccccpccrcccccccccccccccccc 


phGE ry 

QUALITY 


cc 

cc 

cc 

cc 

cc 

cc 


USER SUPPLIED LEaST SQUARES METhOO FQllOwS hFRE TO DETERMINE 
THE COEFFICIENTS OF FQUATICnO (?.2,17) *nO 2,2.76), THf 
subroutine TFL3Q BELOW IS THC Im5L LEaST SQUARES FUNCTION 
FTT routine 


cc 

cc 

cc 

cc 

cc 

cc 


ccrccccrccccccrccrcccccccccccccccccccrccccccccccccrccfcccccccccccccccccc 
CALL IFLSO(F»R»Y,NR,cOEF»NCCCF,wK»IER) 
IF(IER,En,i?R,OR.IER,EO.no)URiTE(6,lo) 
lO F0 RHaT(5AH terminal error in least squares NETHoO*3uBRrUT:r;E COEFS 

n 

RFTtIRN 

ENO 

ccrcccccccccccrccrcccccccccccccccccccfccccccccccccrccrcccrcccccccccccccc 
cc cc 

cc This FUNCTION evaluates THE CfRO order BeSSEI. FllliCTlON cc 

cc DENOTED IN notation n 2»1 (IH)) - SEE FInAl REPORT CC 

cc . . . cc 

ccrccrcrccccccr.ccrcccrcccccccccccccrcrccccccccccccrccrcccccccccccccccccc 

REAL FUNCTION F(N,R) 

COMMON/CJ/R|.AMO(20)*JI C20),J1LA“{20) 

X«RLAMD(N)*R 
CALL JOlX.V) 

F.Y 

return 

ENO 

ccrccccrccccccrccpcccrcccccccccccccccrcccccrccccccrccccccccccccccccccccc 
cc cc 

cc this SUBROUTINE COMPUTES THE JO BESSEL FUNCTION Y«Je(X) CC 

cc . . . . . cc 

ccrccccrccccccrccrcccccccccccccccccccrccccccccccccrccccccccccccccccccccc 
SUBROUTINE ,I0(X,Y) 

cccccrcrccccccrccrcccrcccccccccccccccrccccccccccccrccccccpcccccccccccccc 
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cc cc 

cc use# SU^PLlFO Jfl FUNCTION PUACEP Nf*C, In THIS rXANpLC# THE CC 

CC IMSL flCSSCL FUNCTION mm08J0 IS »LLUSTBaTEO CC 

CC . - . CC 

tcrccccrccccccrccrcccrcccccccccccccccrccccccccccccr.ccrcccrcccccccccccccc 

RF*U nmBSJO 
TtMMeSJOfX# TER) 

RETURN 

ENO ... 

ccrccrcrccccccrccrcccrcccccccccccccccrcccccccccrccrccrcccccccccccccccccc 


cc 

cc 

cc 

cc 

cc 

cc 


SUBROUTINE FOR SOLVlNfi A SYSTEM OF LINF*R SINULtANECUS 
equations HaVINC a TRIOIA60NAI eOEfFIcIEnT naTRtX, oIaCONALS 
ARE STORf.0 TN THE ARRAYS A, 0, aNO C. THe CONRUfEO 
SOLUTION vector V(1) . , . V(t) IS stored !N The array V, 


CC 

cc 

cc 

cc 

cc 

cc 


ccrccccrccccccrccrcccccccccccccccccccrcccccccccccccccrcccrcccccccccccccc 

subroutine TRIOAC(L) 

COMMON/C30/A(50('),Bf590)»C(500),0(SOO),Y(Soo)»BETA(-«05),CAf'HA(Sn5) 

ccfcccccccccccrccrcccrcccccccccccccccrccccccccccccrccrcccccccccccccccccc 
cc cc 

cc compute INT^-RHEO’aTE arrays QETa and gamma cc 

cc . • cc 

ccrccccrccccccrccrcccrcccccccccccccccrccccccccccccrccrcccccccccccccccccc 

BCTA(n«fl(l) 

GAMMAIDaOn )/8ETA(1) 

IFPt*2 

on 1 I«IEP1,L 

BCTA(i)«Ba)-A(i)*ca-n/BCTA>i-n 

CAMHA(T)«TO(H-A(I)«CAMMA{l-n)/BETA(I) 

I CONTINUE 


c COMPUTE etnal soun, vECfcR V ORIGINAL PAGE 19 

vtl)»cahma(,) of poor quality 

00 ?. k«i,last 
I»L«M 

v(i)«CAMHA(n-c(i)*vaTt)/ 0 £TACi) 

2 CONTINUE 
RETURN 
ENO 

ccrccccrcccccccccrcccrccccccccccccccccccccccccccccrccccccccccccccccccccc 


cc cc 
cc PURPOSE - cc 
cc I, generate melt zone surface control function cc 
cc 2 , generate thermal OISTRIOUTiOn in melI IOnE cc 
cc . . . cc 


ccrccccrcccccccccrcccccccccccccccccccrccccccccccccrccccccccccccccccccccc 
subroutine MELT) 
integer nEL'ERM»OELNSTS 
real JImMLaM.MMBSJO 
COHMON/Ct/Rl AM0(Z9)«J1 (20),JILAM'20) 

COMMON/C5/Rflon .P9I (20# 101), 30.M 120) 

COMMON/CR/cnEE(20)»RM 

COMMCN/C 10 /ASCRIP{ 20 ).eSCRIP<? 0 > 

COMMON/C20/A <500 ). 8 (SOO), C (500), Q (MOO), veSflO), beta (R05), gamma (SOB) 
COMMON/CRl/TMeTA8(20,50M),THCLnMOj)#T(Io),GR*Oxi;()nn#GRAOxO(tni) 
C0MM0N/C?2 /ICaSE,TmELT V J) 

CPMMOW/C23/r.RA02MOn#CRAOJ<lOl> 

C0KH0N/CP4/FKS,RML»RL#N3TS 

COMMON/REAOi/P#mtERM,m8UM,XO,J(N,MCRIO,nR 

COMMON/FIXRT/IOPTION,XMlN,CMlN,r|,IP 

COHMON/C25/eRAOATO(t01)#6RAOATQMOl)#RNSt (2 O)#RhS2(20),S(20)#Q# 
1AMAT1(20,101#AMAT2(20# )0),AL2(«p,20),RhS (oa),HNi;cU00),IlHf:(20) 
COMMON/C2 */cPOLY(20) 

COMMON/C27/MERNeL#NMeRNEL#KKGRNE4 
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0 ‘ 


CnMMON/CSO/F(«)>AL6(a«,20),Mli6(fl4) 

CnM‘*0N/CSl/M*xTFBM,»'r*'TCilM,«AXN<Y8,‘<lNNSyS,f>FtTrrM,nELf.'3YS 

ccrccccrcccccccccrcccccccccccccccccccrccccccccrcccrccrcccccccccccccccccc 
cc cc 

CC GCMFRiTr BESSEL E*P*nSION COCFF »C IE^'TS OF EOUATtCNS («.07) iNO CC 
cc or FINAL »C»OBT cc 

cc cc 

ccrccrcrcrccccrccrcccccccccccccccccccrcccccrccccccrccrcccccccccccccccccc 
on 20 

AfDmCRAOATnCD-CCAOATOClOI) 

Ba)t6RAnATi}(n*GBAOiTO(10I) 

20 cnNTIHUC 

'cltt oR'g'nal page ra 

mterh.maxtcpm of poor quality 

NSYS«MAXnSYS 
on 50 lal.NTERM 
ITaI*a 

ccrccrcccccc.i.rccrcc'Tcccccccccccccrccccccccccccccrccrcccccccccccccccccc 
cc , cc 

cc gfnfrate right hand sioes nF fqiution (“.o.e.tj nr final report cc 
cc cc 

ccrccccrccccccrccrcccccccccccccccccccrccccccccccccrccfcccrcccccccccccccc 

RM3fII)«-0,R*A3CRIP(I)*Jl (I)*RLAHO(IJ-GRAOATO(ion 
ITalltMTFRM 

RH3(II)» 0,RaBSCRTP(I)*J 1 (I)«RLAHO(I3*6"ADATO(lflI) 

8(I).iP*P*8.(1«RL*H0(I)»RLAM0(1J 

Sa)«S0RT(8fI)) 

30 CnNTiNUE 

ccrccccrccccccrccrccccccrccccccccccccrccccccccccccrccrcccccccccccccccccc 
cc cc 

cc generate right hand sioes CF EOmATIONS fo.O.lR) " (4,0.22) Of CC 
cc FINAL REPORT cc 

CC . . CC 

ccrccccrcccccccccrcccccccccccccccccccpccccccccccccrccrcccccccccccccccccc 

RHSimo.o 

RMS(2)»0,0 

RHS(3)aGRAOATO(tOt) 

RH3(4)cGRA0AT0(t01) 
on 40 Nal,HTEPM 

on so Kai.NSYS 

NNaN «4 

CALL INTEGL’ (N»K»*R3) 

ccrcccccccccccrccrcccccccccccccccccrccccccccccccccrccfcccccccccccccccccc 
cc cc 

cc generate coefficients in LEPT Hand side of equation (4'.0,23) cc 

cc OF final report cc 

cc ... .. .... cc 

ccrcccccccccccrccrcccrccccccccccccccccccccccccccrcrccrcccccccccccccccccc 

AL2(NN,K)bANS 

CALL INTfGL2(N.K,aNS) 

NNaNN^HTERH 

AL2(NN,K3aANS 

so coNTnuE 
40 continue 
or 60 Kbi.nsYS 

ccrccrcrccccccrccpccccccccccccccccccccccccccccccccrccfcccccccccccccccccc 
cc cc 

cc generate coefficients in lEPT Hand side op cc 

cc EQUATIONS (4.0. 1«) - (4.0,22) CC 

, .. . . cc 

ccpcccccccccccrccrcccccccccccccccccccrccccccccccccrccrcccccccccccccccccc 

call 8aSIS(K,0.ANS) 

A( 2n.l()BANS 
call BASTSCK.O.ANS) 

AI.2(2,K)bANS 
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CALL 08A»ISrK,0>AMS) 

A(.2(3fK)aAN.1 

CALL OBAiitSrK,a<ANS) 

AL2r«.K)aANS 
60 CnNTIMUf 

on 810 Mj(YS«HtNN3YS>MAXNSVS<Ci;t.'‘'3YS 
on 500 MTERM«HINTCRf*,M*XTeR*t,OKl term 
NNaat2*HTeRM 


if(mn.l£,njvs) '".ornsoo 
on 540 ral«A4 
RM4(naRHS(T) 
on 550 Jal«R0 
A(0 <I,J)bALR(!, 1} 

580 cnNTiNue 
540 cnNTiNue 

on 570 Ial,MTER5 

I25MAXTERMt4-»l 

I65MTERM54tT 


ORIGINAL PAGE (S 
OF POOR QUALITY 


on 5«0 Jat ,?0 
AL6(l6,J)sAI 2(12, J) 

56C cnNTINUE 

RH6(16)sRHSM2) 

570 CnuTlNUE 

Efnao.O 

E(2)a0.0 

C(3)»0,0 

E(4)a0.0 

ccrccrcrccccccnccrcccccccccccccccccccfccccccccccccrccrcccpcccccccccccccc 


CC SnLYC FOR LFA3T 30UARE3 30l.UTTnn OF EQUATIONS (J,3,tR) • CC 

CC - (4,0.23), The IM9t "C'TifiE lLfOF 13 ILlUSTRATfO meJ’E CC 

CC . CC 

ccrccccrccccccrccrcccccccccccccccccrcrcccccccccrccrccrcccccccccccccccc'“c 

CALL LLROF(iL6.a4,NN,MSY3,nH6,4fl,i,o,E.CPOLY,20,H55,Yi;K.ICR) 

CC OISPLay twE crEFFirlENTS OF CQU^TION (4,0,18)< the melt ZOr.’E CC 

CC surface control temperature « SfE FTNal report CC 


wRITECO.PO) 

30 FnRMAT(lHl,?lX,7RMH ELT 7nNE 3URFACE CQHTR 
• nt COFFFiCiCfiTSi 
p»ITE(fc,789) hTERH,N3Y3 

t 80 FnH5AT(//,iM ,50X,12HF0R MTERh ■ ,I2,i2H and msvS t ,12) 
m"ITE(6*90) 

99 FORmATC///, tH »«9X, 1HK,22X,4HCCK)) 
on 85 rai.NSYS 
M»ITE(6#»ie6) I» CPOL»(I) 
a 86 FnRMAT(/,lH , I2» 1 OX, E20 , 1 OJ 
85 CONTINUE 
CALL foster 

508 continue 
518 continue 
return 
end 

ccrcccccccccccrccrcccccccccccccccccccccrccccccccccrccfcccccccccccccccccc 

CC CC 

CC PURPOSE • CC 

ec - COMPUTE Matrix tat elemeuts c\w equation (o'.o.js) cc 

CC . cc 

ccrccccrccccccrccrcccrcccccccccLcccccrccccccccccccrccrCccccccccccccccccc 
subroutine INTEQL1(m,K,AN3) 
external c 

CnN-ON/C28/r.RA0AT0n0l),6RAOATQrt0l),RMS,(28),nHG2(?O),S(20),Q, 

lAMAri(20,101«AMAT2(20«10),AL2(a4,2o),RN8(«4),MUi((15oO),IZUt;(20) 

CnHMON/C27/5ERNEL.N5ERNEL,KKCRNPL 

KERNEL*! 

MKERNELaN 


kkeRNEL*K 
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ifR»« 0,0 
RfRP«t ,OF-lrt 
RFR»«t ,OF-l? 

*NS* 0 C* 0 RC(t;, 0 ,Q,AEBR,RCRR,Er,»no, lER) 
return 


ORfGfNAL PAC^ fS 
OF POOR QUALITY 


EMO 

ccrccccrccccccrccrcccccrcccrcccccccccrccccccccccccrccfcccrcccccccccccccc 


cc cc 
cc PURPOSE . cc 
cc • COHPUTE MATPIX ?A? CLEME!;TS fir EOUATIPn CC 

cc cc 


ccrccrccccccccrccpcccccccccccccccccccccccccccccccccccrcccccccccccccccccc 

SU8ROUTINC TNTE6L2(*'.'<.*NS) 

COMM0M/C25/CRiC*T0f lOI)»CRiOATQriOj),PHSi f 20 1 .RuS2 t? 2 ) , 3 (20 ) »0 f 
< AMArt( 20 ,l 01 «iMAT 2 ( 20 « 10 ).AL 2 rPP« 20 )>PHS(aa),MMK(iSoO)>IlMr( 20 ) 
CnHMON/C^T/KEPNrL.N^EPNEL.KKEPNrU 

external c 

KFRNEL=2 
NKERNELaN 
KKERNelsK 
AFRRso.O 
rerrs 1 , of * 1 0 
RFRPai ,OF-l? 

ANSaOCAORE((:, 0 ,(J,AeRR,RERR,ERPO«».lER) 

RETl/RfJ 

END 

ccrccrcrccccccrccrcccrcccccccccccccccrccccccccccccrccccccccccccccccccccc 


cc cc 
cc PURPOSE . cc 
cc • Evaluate the kernal FUncTICn nCFivEn 8 v ERUATjfiN (o.C.joj CC 
cc cc 


ccfccrccccccccrccrcccecccccccccccccccrccccccccccccrccrcccccccccccccccccc 

subroutine kePncl1(n,t,*n3) 

CPMHOU/C25/r.RAOATO( lOU ,CRAOATQ>jot ),PhSi (20),Rm 32(20),S(23) »0. 
1AMATH20, tfl1,AMAT2(2n,10),AL2C4ii,20)»»HS(O«),»«wiCt5e0),IIH|;(20) 
CnMHON/REAOl/P,“TE»N,"SUH,XO,XH.KCRIC,sR 
Z»-«»S(N) 

TFRMa(p*P«3(N)»S(N) )/«.0 
RHOLOaO.n 

IF(Z,CT.-25o.O)PMOLOaEXP{2) 

TrRHsTERH*(t .O-R hOLO) 

Za.(P*S(NI ),T/ 2,0 

Zla-S(N)«Q*fS(N)-P2*Ty2,0 

Rlan ,0 

R2»ff,0 

irCZ,CT.- 25 O. 0 )Rlae*P(Z) 

irCZj.CT., 2 so.O>P 2 »CxPCZO 

AN3«TeRM,(R,.P2) 

rfTurn 


ENO 

ccrccrccccccc-Tccrcccccccccccccccccccrccccccccccccrccrcccrcccccccccccccc 


cc cc 
cc PURPOSE - cc 
cc • evaluate the VePNAl function nEFINEO 8 r EOUATyfiN (O.OtJZ) CC 
cc CC 


ccrccrcrccccccrccrcccrcccccccccccccccrccccccccccccrccrcccrcccccccccccccc 

SUBROUTINE kEPnEL 2 (‘<»t,an 3 ) 

CnMHON/C2S/r,RAOATO(lOt),CPAOAT(3TlOn.RMSi (20 ) , (?i02 (20 ) , 3 (20 ) » 0, 
lAH 4 TU 20 , 102 »AH 4 T 2 ( 2 o.l 0 )»AL 2 (an» 2 O 2 ,»M 3 (a(,).owi(iSoO),Hv»;( 20 ) 
CPMHON/PFAO;/Pf NTFPH,HSUH.XO,*N,;iGRI0,NP 
T)»(P«P.S(ni* 3 (n) I/R.O 
Za-sA(M)*fl 
RHOLDaO.rt 

IF{Z,CT.-25n,0)PMOL0«EXP(Z) 

TlaTi/(t,o»PHOLO) 

Za(P-3(N))«n,5«(0-T) 
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ORIGINAL PAGE IS 
OF POOR QUALITY 


T7S0.0 

rF{Z,i:T.- 25 n. 0 )T 2 « 6 *P( 2 ) 

Z»S(N)*T 
T5«l,0 

ir(Z.CT.-<Sn.0)T3«l.o-CXPt2) 

*N8«TI»T2*T1 

upturn 

CNO 

ccpccrctccccccr.ccrcccrcccccccccccccrcrccccccccccccrccccccccccccccccccccc 
cc cc 

cc PURPQse * cc 

cc - cnHPuTp iMTfcRANo useo Tc compote h*tpi» ?*? of cc 

cc EOUaTIOM cc 

cc .. . cc 

tcrccccrccccccrcccccccccccccccccccccccccccccccccccrccccccccccccccccccccc 

«»r*c FUNCTION c(X) 

COMMON /CaT/KeBNFL.N«E"NeL.KKCBNFL 
RF*U X 

IF(FE«N6|..E0,I) CALt KEBNEL1 (NKFKNEL»*.*nS) 

IF(KeRNEl..Eo,Z) call KERNELa<NKFRNEL.X.*NS) 

CALL HASI5(KRERNEL.*,*NSn 
CaANS*AN51 


return 

end ..... 

ccrccrcccrccccr.ccrcccrcccccccccccccccrcccccccccccccccrcccccccccccccccccc 


cc cc 
cc purpose • cc 
cc - provide user supplied set Cf functions used In EXPANCION cc 
cc OF nclt zonf surface control function. See ebuation (a.o.jai cc 
cc . . , cc 


CCPCCCCrCCCCCCCCCrCCCCCCCCCCCCCCCCCCCrCCCCCCCCcCCCCCCrCCCrCCCCCCCCCCCCCC 
SUBROUTINE RASIS(K,T,anS) 

CnHMOM/c2S/ORAOATO(IOI)fCR«0»Ta'rlOnrRNS| (2 o).RhS 2(2O),S(Z0)>Q. 
lAMATl(20,tOnANAT2(20f l0),AL2taA,20)»BMStaa),Nw5(15o0),IINKf20) 
iPCt.Nc.nco TO JO 
ANSal.O 
RFTIIRN 

10 ANS»(T«0/2.0)»«(k-1 ) 

return 
END 

subroutine nBAStS(B»T.AKSJ 

OFLTaTiO.OOI 

XpT*OELTAT 

call 8ASTS(B»X,ANSt) 

XpT-oELTaT 

'■ALL BaSIS(b.X,ANS2) 

aNS»{ANSi»AN32)/(2.0*OELTAT) 

rfturn 

end ... 

ccrccccrccccccrccrcccrcccccccccccccccrcccccccccccccccrcccccccccccccccccc 


cc TnIS SUBPOUTINC estimates LaTeRaL surface temperature cc 
cc BV USE OF C«ialc SPLINE IF The UREP SUPPLIES A DISCRETE SET OF CC 
cc LaTePal surface TEmPEPATURCS. CC 
cc cc 


ccpccccpccccccrccrcccccccccccccccccccrccccccccccccrccrcccccccccccccccccc 
subroutine SPLINE(»^vT,YINT) 
cOMNON/c^z/xonooj.Toi joo},c(o, «oo),M 
IF(XINT-X 0 (in 2 »l *2 

1 YINTarod) 
return 

2 v*l 

3 If (X lNToXD(KRl)10fa,S 
« YlNTaY0(«*n 

return 
5 K8K«1 
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0 - 


■ir((M-K),cT',o) RtiToJ 

ir 

6 rtNT»(X0fK>O«xrwT)*(C(l,K)«()(nfK^l)»»l'*Tl»*?*f ^3.*)) 
TTNT«yiVT*(xINT»xnC») )«fC(2»K)*i-XlNT-X'){K) )«*2trC«.tO ) 

RETUffN 

END 

ccrccccrccccccrccrcccrcccccccccccccccrcccccccccccccccrcccrcccccccccccccc 
cc cc 

CC The 3Pt INE PIT COrErTCIC^TS. for use in CONJU'.'C'ION cc 


cc nITh SUBROUTINE SALINE. CC 
cc INPUTS • • cc 
CC N S no, GF DAT* P*l»s cc 
cc xo » *PR*Y OF X f*8ClSS*) V.'.LUfS CC 
cc TO s *RR*Y OF Y (oROINATES) V*| UES CC 
cc OUTPUTS - cc 
cc C « 2 -n *PR*Y OF SPLINE FIT C^CFFIC TE'^tS (a Cf'CFFiCIt NTS CC 
CC PER triplet of DAT* POINTS), CC 
cc . . cc 


ccrccrccccccccrccrcccccccccccccccccccrcccccccccccrrccrcccrcccccccccccrcc 
subroutine rOFCEN 
c 

coMHaN/C32/xo(ino).Toaoo),ci(fl, ioo).'< 

OIHENSION era, too) 

DIMENSION Pri00),En0O),*(l00,3'.3{i0o)»Zfino),r>Ci0n) 

EOUIv*lence re t ( t * 1 ) »c r I , I ) ) 

c 

NOSM 

M»M-l 

no 2 K«i,M 

nCF)»XDrK*\ )-XD(X) 

P(K)*0(K)/n. 

2 r CK)s(vf^(K*n-YO(t<) )/o(K) 

DO 3 K*2«H 

3 B{K)aEtK)-F(x-«) 

A(l. 2 )«-l,- 0 (U/ 0 ( 2 ) 

An,S)«n(n/o( 2 ) 

A'2,2)*2.*rPtn*P{2))-P(n»H(i,2) 

AC,-. :;*fP{p)«p(i)*A(t,3))/*(2,?j 
BC 2 )a 8 ( 2 )/A(?, 2 ) 

DO 4 KX3tH 

A(K,2)«2.*fPtK-n»PCO )»PC«-n«A(K-t,3) 

B(K)»8tK)-P(lc-l )«B(K-l) 

A(K,3)*P{(0/A(K,2) 

4 P{K)* 8 (K)/*(K, 2 ) 

O«D(M«l)/0fM) 

A(nD» l)st.* 0 «A(M»lf 3 ) 

A(N 0 . 2 )>- 0 «A(Nn.l)*A(H, 3 ) 

P(N0)«8(N-» )»*(NO» n*B(M) 

Z(NO)aB(NO VACND.2) 

no 6 I> 1 ,ND «2 

KsND*I 

0 ?(K)aa(K)**(K. 3 )*Z(KTl) 

Z(l)*»An,?)«2C2)-*n*3)*2f3) 

DO 7 K«t,H 

Dsl./( 6 .*OrK)) 

C(l,K)a7(K)«0 

C{2,K)aZ{K«n*Q 

CC 3 »K)avO(K)/OfK)-ZC< 5 *PfK) 

7 C(a,K)avO(x*n/0(F)-ZCXtl )*P(Ki 
MaM«t 
•'TURN 
END 


Figure C-3. Computer Code List For 
Problem Pl-3 (Cone) 
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